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ABSTRACT 


The  potential  of  degenerate-four  wave  mixing  (DFWM)  spectroscopy  for  the 
analysis  of  molecular  species  in  high-pressure,  high-temperature  environments  has  been 
explored  theoretically  and  experimentally.  In  addition,  the  reliability  of  DFWM  has  been 
assessed  by  comparing  to  laser  induced  fluorescence  (LIF)  measurements  under  the  same 
experimental  conditions.  The  results  of  these  studies  suggest  that  DFWM  can  provide 
quantitative  information  regarding  trace  molecular  species  in  environments  in  which 
collisions  and  source  emission  hinder  analysis  by  other  means. 

DFWM  has  been  used  to  detect  the  CH  radical  in  an  atmospheric-pressure 
oxyacetylene  flame  via  the  (0,0)  and  (1,1)  bands  of  the  A^A-X^II  system.  The  CH  radical 
is  a  minor  flame  species  and  is  important  for  understanding  the  primary  reaction  zone  of 
many  combustion  environments.  We  have  observed  CH  radicals  with  comparable 
sensitivity  by  both  DFWM  and  LIF.  From  these  measurements,  we  estimate  a  DFWM 
detection  limit  of  4  x  lOH  molecules/cm^  (4  x  10^  molecules/cm^  per  quantum-state)  for 
CH  at  atmospheric  pressure.  Vibrational  temperatures  and  concentration  profiles  of  CH 
obtained  by  both  techniques  are  in  good  agreement. 

Diagrammatic  perturbation  theory  combined  with  a  spherical  tensor  treatment  is 
used  to  decompose  the  DFWM  signal  resulting  from  an  isotropic  molecular  sample  into  a 
sum  of  three  multipole  moments  in  the  weak-field  (no  saturation)  limit.  The  zeroth 
moment  gives  the  relative  internal-state  population  contribution,  the  first  moment  the 
orientation  contribution,  and  the  second  moment  the  alignment  contribution  to  the  DFWM 
spectra.  This  treatment  makes  explicit  how  the  magnitude  of  the  DFWM  signal  depends 
on  the  polarizations  of  the  other  three  beams  and  the  collisional  relaxation  caused  by  the 
environment.  A  general  expression  is  derived  for  the  DFWM  signal  for  an  arbitrary 
geometric  configuration  of  the  beams  (arbitrary  phase  matching  geometry).  Under  the 
assumption  that  the  rates  of  collisional  relaxation  of  the  population,  the  orientation,  and 
the  alignment  are  the  same,  simple  analytic  expressions  are  found  for  the  most  commonly 
used  experimental  configurations,  which  should  facilitate  the  practical  analysis  of  DFWM 
spectra. 

Effects  of  field  polarization  and  field  intensity  in  DFWM  have  been  investigated, 
and  rotational  distributions  (temperatures)  have  been  obtained  for  the  CH  radical  using 
saturated  DFWM.  As  the  field  intensity  increases  optical  pumping  rates  become  larger 
than  the  population  relaxation  and  collisional  dephasing  rates,  and  the  DFWM  signal 
saturates  becoming  relatively  insensitive  to  collisions.  The  saturation  behavior  predicted 


by  nondegenerate  two-level  models  is  in  close  agreement  with  the  observed  power 
dependence  of  the  (0,0)  band  transitions  of  the  A^A-X^IT  system.  The  DFWM  saturation 
intensity  does  not  change  significantly  as  the  polarization  states  of  the  excitation  fields  are 
changed.  However  large  differences  in  the  DFWM  signal  intensities  are  observed  as  a 
function  of  input  field  polarization  and  rotational  branch.  These  differences  are  nearly 
independent  laser  intensity.  The  DFWM  signal  intensity  differences  are  rationalized  using 
diagrammatic  perturbation  theory.  The  weak-field  expressions  reproduce  the  data  taken  at 
saturating  intensities  with  an  accuracy  of  10-30%.  The  important  aspects  of  the  reduction 
of  saturated  DFWM  signal  intensities  to  relative  internal-state  distributions  in 
environments  where  population  relaxation  and  coherence  dephasing  events  are  dominated 
by  collisions  are  outlined.  Rotational  temperatures  obtained  by  this  method  are  estimated 
to  be  accurate  to  5%. 
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CHAPTER  1 
Overview 


-  Into  the  storm  - 


In  between  the  extremes  of  ultrahigh  vacuum  and  the  condensed  phase  lies  the 
regime  of  high-pressure  (several  hundred  Torr  and  above)  and  high-temperature  (several 
hundred  Kelvin  and  above)  gases  of  which  flames  and  plasmas  are  the  most  common 
examples.  The  dynamics  of  these  environments  are  complicated  compared  to  low- 
pressure  studies  because  of  the  high  collisional  rates.  Of  course,  many  "collisions"  occur 
in  the  condensed  phase,  but  the  situation  can  be  reduced  by  only  considering  nearest 
neighbor  interactions;  for  example  the  problem  can  be  simplified  to  a  solute  in  a  solvent 
shell  or  a  dopant  in  a  crystal  lattice.  In  high-temperature,  high-pressure  gases  the 
collisional  rates  approach  those  of  the  condensed  phase,  but  the  mobility  of  interacting 
particles  is  much  higher.  Therefore  a  given  molecular  species  can  interact  with  a  variety  of 
other  molecules  and  surroundings.  In  addition  these  environments  are  "hostile"  in  the 
sense  that  they  are  usually  associated  with  high  levels  of  emitted  radiation,  limited  access, 
nonequilibrium  molecular  concentrations  and  energy  distributions,  and  large  gradients  in 
temperature  and  concentration.  Most  of  these  environments  pose  a  severe  challenge  to 
those  who  wish  to  probe  their  chemical  composition,  kinetics,  and  dynamics  as  a  function 
of  spatial  distribution  and  temporal  evolution.  Even  though  these  environments  are 
technically  challenging,  their  importance  warrants  the  experimental  effort  required  to 
investigate  them. 

The  goal  of  my  thesis  is  to  develop  a  laser-based  diagnostic  tool  for  the 
investigation  of  trace  molecular  species  and  reaction  processes  in  collisionally  dominated 
environments.  Specifically,  the  development  of  a  technique  that  can  be  used  for  the 
investigation  of  the  plasma  chemistry  existing  in  the  flow  of  an  atmospheric-pressure  rf 
inductively-coupled  plasma  that  is  part  of  a  materials  processing  facility  in  the  High 
Temperature  Gasdynamics  Laboratory  (HTGL)  at  Stanford.  This  goat  is  accomplished  by 
exploiting  a  relatively  new  nonlinear  spectroscopy  called  degenerate  four-wave  mixing 
(DFWM).  The  results  are  complemented  by  companion  studies  using  laser  induced 
fluorescence  (LIF). 

Reactive  plasmas  are  harsh  environments  important  to  energy  science  and 
technology.  The  nonequilibrium  behavior  of  thermal  plasmas  is  of  particular  concern  in 
the  area  of  plasma  chemistry  and  plasma  processing.  Previous  experiments  have  probed 
the  nonequilibrium  plasma  state,  excitation  temperatures,  electron  densities,  radiation 
escape,  recombination  rates,  and  electronic  quenching  effects.  Emission  and  LIF 
measurements  have  met  with  limited  success  because  of  the  brightness  of  the  source  and 
the  high  temperatures  and  pressures  within  the  plasma.  A  need  exists  for  in  situ 
nonintrusive  techniques  that  can  probe  nonequilibrium  behavior,  characterize  related 
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chemical  properties,  and  monitor  trace  molecular  species  free  from  source  interference  and 
collisional  effects.  The  experimental  and  theoretical  efforts  presented  in  this  thesis  have 
been  and  are  currently  being  applied  at  HTGL  to  understand  the  gas-phase  and  gas- 
surface  interactions  encountered  in  the  boundary  layer  of  their  reactive  plasma. 

Operated  on  a  variety  of  gas  mixtures,  especially  CH4/H2/Ar,  the  HTGL  plasma 
reactor  has  been  used  in  the  past  few  years  to  gain  experience  in  studying  chemical  vapor 
deposition  of  diamond  thin  films.  A  small  scale  "reactor"  was  set  up  in  our  laboratory  that 
consisted  of  an  oxyacetylene  flame  that  was  operated  under  conditions  known  to  be 
capable  of  growing  diamond  thin  films.  The  temperatures,  collisional  rates,  and  CH 
radical  concentrations  are  very  similar  for  the  flame  and  the  HTGL  plasma,  so  experiments 
that  emphasized  the  investigation  of  DFWM  as  a  diagnostic  were  conducted  in  the 
relatively  controlled  environment  provided  by  the  flame. 

DFWM  was  chosen  because  it  offers  a  means  of  measuring  spatially  resolved 
temperature  and  concentration  information  with  a  very  high  degree  of  spectral,  spatial,  and 
temporal  resolution.  In  addition,  the  problem  of  plasma  luminescence  is  largely  overcome 
with  DFWM  because  it  provides  a  coherent,  directional  beam  that  can  be  detected 
remotely  without  loss  in  intensity.  In  Chapter  2,  DFWM  is  evaluated  with  regard  to  its 
potential  to  provide  sensitive  detection  of  molecular  species  and  quantitative 
determinations  of  CH  vibrational  temperatures  and  relative  concentrations  in  the  primary 
reaction  zone  of  an  atmospheric-pressure  oxyacetylene  flame.  DFWM  is  found  to  be  as 
sensitive  and  as  accurate  LIF  under  our  experimental  conditions. 

Vibrational  temperatures,  however,  may  not  be  equilibrated  with  the  translational 
temperature.  Therefore  rotational  temperatures  are  needed  because  they  are  more  likely 
to  represent  the  translational  temperature  of  the  environment.  The  ability  to  reduce 
DFWM  rotational  spectra  to  relative  populations  is  key  to  the  ultimate  success  of  the 
technique.  Chapter  3  shows  how  the  magnitude  of  DFWM  signal  resulting  from  an 
isotropic  molecular  sample  depends  on  the  degeneracies  of  the  levels  involved,  the 
polarizations  of  the  excitation  fields,  the  phase  matching  geometry,  and  the  collisional 
relaxation  and  dephasing  caused  by  the  environment.  Simple  analytic  expressions  are 
provided  which  should  facilitate  the  practical  analysis  of  DFWM  spectra.  However 
obtaining  an  accurate  and  complete  set  of  relaxation  rates  for  a  given  molecule  in  a 
specific  environment  is  a  formidable  task.  This  difficulty  results  because  these  rates  are 
typically  dependent  on  temperature,  collision  partner,  and  quantum  state,  and  therefore, 
are  not  known  for  a  wide  range  of  conditions. 


3 


An  appealing  aspect  of  DFWM  is  that  the  signal  intensity  becomes  relatively 
insensitive  to  the  specific  value  of  the  relaxation  rates  in  the  saturation  regime.  In  this 
regime,  the  problem  of  extracting  relative  population  distributions  essentially  reduces  to 
knowing  the  absorption  coefficients.  This  aspect  is  discussed  in  detail  in  Chapter  4  where 
DFWM  signals  are  presented  as  a  function  of  field  polarization  and  field  intensity.  Also  in 
Chapter  4  the  important  aspects  of  the  reduction  of  saturated  DFWM  signal  intensities  to 
relative  internal-state  distributions  are  outlined  for  environments  where  population 
relaxation  and  dephasing  events  are  dominated  by  collisions,  and  a  rotational  temperature 
analysis  is  presented  of  the  CH  radical  in  an  atmospheric-pressure  oxyacetylene  flame. 
Rotational  temperatures  determined  using  the  procedures  discussed  in  Chapter  4  are 
estimated  to  be  accurate  to  5%. 


CHAPTER  2 


Evaluation  of  DFWM  as  a  Spectroscopic  Molecular  Probe 


A  person  with  one  thermometer  knows  the  temperature. 
A  person  with  two  is  not  sure. 


Kermit  C.  Smyth 


L  INTRODUCTION 


Highly  luminous  sources,  such  as  arcs,  sparks,  flames,  explosions,  plasmas,  and 
discharges,  pose  a  severe  challenge  to  those  who  wish  to  probe  their  chemical 
composition,  kinetics,  and  dynamics  as  a  function  of  spatial  distribution  and  temporal 
evolution.  Nonintrusive  optical  methods,  ^  like  LIF,  are  hindered  by  interference  from 
background  emission,  distortions  caused  by  source  temperature  and  density 
inhomogeneities,  and  collisional  energy  transfer  and  quenching  processes  that  readily 
occur  at  high  pressures.  One  method  to  overcome  these  difficulties  is  a  four-wave  mixing 
technique  called  coherent  anti-Stokes  Raman  scattering  (CARS).  1.2  xhis  can  provide  a 
coherent  signal  propagating  in  a  known  direction  so  that  its  detection  can  be  accomplished 
with  excellent  rejection  of  source  emission.  Another  nonlinear  technique  that  shares  the 
same  characteristics  is  DFWM^  with  the  extra  advantages  that  (1)  DFWM  involves  one- 
photon  resonances  and  is  therefore  more  sensitive,  and  (2)  DFWM  employs  a  single 
frequency  provided  by  one  laser  source  and  therefore  is  easier  to  implement. 

DFWM  is  a  third-order  nonlinear  optical  process  in  which  three  fields  of  a  single 
frequency  ©  overlap  in  a  medium  to  produce  a  coherent  fourth  field  which  also  has 
frequency  ©.  This  process  may  also  be  understood  in  terms  of  light  scattering  from  a 
grating.  We  employ  the  collinear  phase-conjugate  geometry  which  consists  of  two  coaxial 
and  counterpropagating  pump  fields,  denoted  (forward)  and  (backward),  which  are 
crossed  at  a  small  angle  0  by  a  probe  field,  E^,.  The  grating  picture  of  DFWM  involves 
the  interference  of  E^,  E^,  and  E^.  When  the  frequency  is  tuned  to  an  atomic  or 
molecular  resonance,  the  interference  leads  to  a  spatial  modulation  in  the  (complex) 
refractive  index,  and  hence  forms  a  grating.  Two  important  modulation  patterns  result,  as 
are  shown  in  Fig.  1 :  E^  and  E^,  form  a  grating  that  causes  E^  to  be  scattered  into  E^ 

(the  DFWM  signal),  which  is  the  phase  conjugate  of  E^  (£,  oc  £"*);  and  E^  and  E^  form 
a  finer  grating  with  spacing  and  orientation  to  scatter  E^  into  E^ .  A  third  grating,  that 
arises  from  the  interference  of  Ej  and  does  not  contribute  to  this  process  because  it  is 
not  phase  matched  for  the  scattering  of  Ep  into  E^.  The  depth  of  modulation  of  the 
gratings  is  minute,  and  their  scattering  efficiency  is  poor  (about  1:10^  under  the  present 
conditions),  but  Ej^  and  E^  have  so  many  photons  that  E^  is  readily  detected. 

Although  DFWM  is  in  its  infancy  as  a  nonlinear  laser  spectroscopy,  it  has  been  used  to 
detect  atomic^"  12  and  molecular  species  1^  “21  in  a  variety  of  conditions.  In  combustion 
systems,  sodium^*^.^  and  lithium  1 1  atoms  have  been  detected  in  atmospheric-pressure 


FIG.  2.1.  Grating  picture  of  DFWM.  (a)  The  interference  pattern  between  forward  pump 
Ej  and  probe  establish  a  large-spaced  grating  relative  to  (b)  the  small-spaced  grating 
established  between  backward  pump  and  probe  E^,.  The  phase-conjugate  field  E^  is 
the  generated  DFWM  signal. 
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Grating  Spacing  Grating  Spacing 

D  =  A,/[2sin(e/2)]  D  =  X/ [2cos  (0/2)] 


flames  with  high  sensitivity.  Ewart  and  O’Leary  were  the  first  to  detect  a  molecular 
species,  OH  radical,  by  DFWM  in  an  atmospheric-pressure  flame,  and  Rakestraw, 
Farrow  and  Dreier  have  quantified  OH^^  in  an  atmospheric  flame  and  NH^^  in  a  low- 
pressure  flame. 

We  report  here  the  first  DFWM  study  of  the  CH  radical,  found  as  a  trace  species 
(30  ppm),  in  an  atmospheric-pressure  oxyacetylene  flame.  From  these  measurements,  we 
estimate  a  detection  limit  (signal-to-noise  ratio  of  1)  of  4x10  D  molecules/cm^  (4x10^ 
molecules/cm^  per  quantum  state)  for  CH.  At  high  pressures,  rotational  distributions  from 
LIF  measurements  become  difficult  to  interpret  because  of  the  rotational  level  dependence 
of  the  fluorescence  quantum  yield,  which  results  from  collisions.  ^2, 23  Therefore  to  make 
a  direct  comparison  between  DFWM  and  LIF,  we  compare  the  (0,0)  and  (1,1)  intensities 
of  a  single  rotational  level  of  the  CH  A^A-X^H  system  to  obtain  CH  concentration  profiles 
and  vibrational  temperatures.  Saturation  properties  of  DFWM  and  LIF  signal  intensities 
and  line  shapes  are  also  presented. 

n.  EXPERIMENTAL  SECTION 
A.  Experimental  Apparatus 

Figure  2  shows  the  experimental  configuration  in  which  three  fields  of  the  same 
frequency  overlap  in  the  reaction  zone  of  an  atmospheric-pressure  oxyacetylene  flame  to 
generate  a  fourth  field  (the  DFWM  signal).  The  two  coaxial  and  counterpropagating 
pump  fields,  denoted  (forward)  and  (backward),  are  crossed  at  a  small  angle,  0=2°, 
by  a  probe  field,  E^,.  The  DFWM  signal  field,  E^,  is  coherently  generated  and  propagates 
6  m  to  a  photomultiplier  where  it  is  detected.  At  the  same  time,  the  LIF  signal  is  obtained 
by  imaging  the  reaction  zone  of  the  flame  on  the  entrance  slit  of  a  0.75-m  monochromator 
and  is  detected  with  a  photomultiplier.  The  specific  detuls  of  the  experiment  are  given 
below. 

Excitation  is  provided  by  a  KrF  excimer-pumped  dye  laser  system  (Lambda  Physik 
EMG  102  MSC  and  FL  2002)  operated  over  the  wavelength  range  of 425-450  nm  with  a 
bandwidth  of  0. 18  ±  0.02  cm*^  and  a  pulse  width  of  16.3  ±  1  ns  (fwhm).  The  bandwidth 
was  measured  with  an  etalon,  and  the  pulse  width  was  measured  with  a  fast  photodiode 
and  a  digitizing  oscilloscope  (HP54510A).  The  output  beam  of  the  dye  laser  is  spatially 
filtered  to  improve  the  beam  profile.  The  beam  profile  was  measured  by  the 


FIG  2.2.  Experimental  configuration  for  the  simultaneous  measurement  of  DFWM  and 
LIF  signals.  Spectra  are  recorded  with  unfocused  beams  (diameter  ~  0.5  mm)  having 
energies  of  12  pJ  for  the  forward  and  backward  pump  beams  and  3  pJ  for  the  probe  beam. 
The  flame  consists  of  a  standard  welding  torch  fitted  with  a  0.94  mm  diameter  nozzle. 

The  oxygen  to  acetylene  volumetric  flow  ratio,  R  =  O2/C2H2,  is  0.947.  The  abbreviations 
used  in  the  figure  are  defined  as  follows:  M:  mirror,  P:  polarizer,  L;  lens,  PH:  diamond 
pinhole,  PEM:  photoelastic  modulator,  BS:  beamsplitter,  EM:  energy  meter, 

BPF:  bandpass  filter,  PMT:  photomultiplier  tube. 
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procedure  outlined  in  reference  24  and  was  found  to  be  well  described  by  a  cylindrical 
Gaussian  beam  profile  with  a  fwhm  of  0.47  mm.  The  laser  energy  is  controlled  using  a 
photoelastic  modulator  (Hinds  International  PEM-CF4)  between  two  crossed  polarizers 
which  allows  the  laser  energy  to  be  adjusted  between  1  and  800  pJ  with  no  significant 
beam  walk  over  distances  of  several  meters.  The  energy  of  the  probe  beam  is  one  quarter 
of  the  forward  pump  beam,  and  the  energy  of  the  backward  pump  beam  is  adjusted  to  the 
value  of  the  forward  pump  beam  using  a  half-wave  plate  and  a  polarization  beamsplitter 
which  reflects  S-polarized  and  transmits  P-polarized  light.  We  measure  the  forward  pump 
beam  energy  by  placing  an  energy  meter  (Molectron  J4-09)  after  the  polarization 
beamsplitter,  which  allows  the  energy  of  the  forward  pump  to  be  measured  directly  and 
accurately  to  a  few  microjoules.  In  calculating  the  laser  spectral  intensity  for  these 
experiments  from  the  measured  beam  energies,  we  use  Av  =  0.18  ±  0.02  cm'^.  Ax  = 

16.3  ±  1  ns,  and  A=%{dlTp-  where  dll  is  the  1/e  intensity  radius  of  a  cylindrical  Gaussian 
beam  ("top  hat"  beam  diameter).25  For  example,  the  beam  energies  employed  in  these 
experiments  are  12  pJ  for  the  forward  and  backward  pump  beams  and  3  pJ  for  the  probe 
beam,  which  correspond  to  1.65  x  10^  W/cm^-cm’^  and  0.41  x  10^  W/cm^-cm’l, 
respectively. 

The  experiment  is  configured  with  the  forward  pump  and  probe  (Ej  and  E^) 
having  the  same  linear  polarization  (P  plane)  and  backward  pump  {E^)  having  a  linear 
polarization  rotated  by  90°  with  respect  to  the  forward  pump  and  probe  (S  plane)  .  The 
DFWM  signal  {E^)  has  the  same  linear  polarization  as  the  backward  pump  (S  plane)  and 
propagates  along  the  same  axis  but  in  the  opposite  direction  of  the  probe.  Since  the  signal 
(£,)  is  produced  in  a  known  direction,  {k,  =  -k^),  alignment  of  the  detection  axis  can 
easily  be  achieved  by  retro-reflecting  the  probe  beam  through  a  retardation  plate.  The 
DFWM  signal  is  collected  with  a  40%  beamsplitter  and  passes  through  a  series  of 
apertures,  a  polarizer  and  a  bandpass  filter  before  being  detected  with  a  photomultiplier 
(Hamamatsu  R2393P). 

The  LIF  signal  is  collected  with  frl  optics  at  90°  with  respect  to  the  pump  beam 
axis  and  focused  on  the  entrance  slit  of  a  0.75  m  single  monochromator  (Spex  1702)  with 
a  1800  gr/mm  grating.  The  entrance  slit  of  the  monochromator  is  set  at  1  mm  x  2  mm 
(width  X  height),  and  the  exit  slit  width  is  set  at  2.5  mm,  which  results  in  a  trapezoidal¬ 
shaped  spectral  bandpass  function  with  a  1.4  nm  bandpass  at  the  top  and  a  2.2  nm 
bandpass  at  the  base.  When  the  grating  is  tilted,  the  monochromator  functions  as  a 
narrow  variable  bandpass  filter  though  which  the  LIF  light  passes  before  being  detected 
with  a  photomultiplier  (Thom  EMI  9558).  The  grating  of  the  monochromator  is 


positioned  to  give  a  flat  bandpass  response  from  430.4  to  43 1 .8  nm  (X/vac)  which  passes 
the  Q(2-12)  transitions  of  both  the  (0,0)  and  (1,1)  bands  while  discriminating  against  the  P 
and  R  branches. 

Signals  fi'om  DFWM  and  LIF  are  amplified  (Lecroy  WIOOBTB)  before  being 
processed  by  separate  boxcar  averagers  (SRS  SR250).  The  laser  energy  is  processed 
directly  by  a  boxcar  averager  (SRS  SR250).  The  three  channels  are  read  at  each 
wavelength,  digitized  and  averaged  over  30  laser  shots.  The  data  is  then  stored  on  a 
microcomputer  (IBM  XT)  for  further  analysis. 

B.  Flame 

The  atmospheric-pressure  oxyacetylene  flame  consists  of  a  standard  welding  torch 
fitted  with  a  0.94  mm  diameter  nozzle  operated  in  the  open  air.  The  flows  of  acetylene 
and  oxygen  are  maintained  with  mass  flow  controllers  (MKS)  at  380  cc/min  and  360 
cc/min  respectively.  This  results  in  an  oxygen  to  acetylene  volumetric  flow  ratio, 

R  =  O2/C2H2,  of 0.947  and  an  equivalence  ratio,  O,  of  2.63  (fuel  rich)  with  respect  to  the 
combustion  reaction 

C2H2  +  5/2  O2  =  2  CO2  +  H2O. 

Joklik26  has  extensively  studied  the  CH  radical  in  an  oxyacetylene  flame  at  low 
pressure.  At  a  pressure  of  40  Torr  the  CH  concentration  decreased  from  23 .6  ppm  at 
0=1  to  1 1.0  ppm  at  0=1.4,  as  measured  by  absorption.  Similarly,  the  CH*  (A^Il) 
concentration  measured  by  emission  decreased  from  0.018  ppm  at  0=1  to  0.009  ppm  at 
0=1 .4.  Furthermore,  the  CH  and  CH*  concentrations  in  ppm  were  found  to  be 
independent  of  pressure  from  20  to  100  Torr  and  to  increase  with  temperature.  These 
low-pressure  studies  give  much  insight  into  the  production  of  CH  in  oxyacetylene  flames, 
but  quantitative  extrapolation  of  these  findings  to  our  experimental  conditions  at 
atmospheric-pressure  is  unreliable.  However,  lessen  and  Gaydon^^  measured  the  CH 
concentration  in  a  premixed  atmospheric-pressure  oxyacetylene  torch  by  absorption 
spectroscopy  and  estimated  the  concentration  to  be  tens  of  ppm.  Bonczyk  and  Shirley^S 
measured  the  CH  concentration  to  be  23  ppm  by  saturated  LEF  and  57  ppm  by  absorption 
in  a  premixed  atmospheric-pressure  flame  with  3)=1.  Finally,  Matsui,  Yuuki,  and 
Sahara^^  modeled  a  premixed  atmospheric-pressure  flame  with  C>  >  2.5  and  measured  the 
relative  C2  concentration  profiles  in  an  oxyacetylene  circular  torch  burner  at  atmospheric 
pressure.  The  agreement  between  the  calculated  C2  concentration  and  their  measured 


values  was  good,  and  the  code  predicted  C2  concentrations  in  agreement  with  absolute 
measurements  made  by  other  workers.  The  equilibrium  CH  concentrations  estimated  from 
their  code  (see  Fig.  10  of  reference  29)  is  20-30  ppm  for  O  ~  2.5.  Therefore,  we  estimate 
a  CH  concentration  for  0=2.63  of  30  ppm  in  our  atmospheric-pressure  oxyacetylene 
flame. 


C.  CH  Spectral  Features 

Vibrational  temperatures  and  spatial  profiles  are  derived  from  the  integrated  signal 
intensities  of  the  well  resolved  Rjg(8),  Rjj(8),  R2e(8)  and  R2f(8)  transitions  of  the  (0,0) 
and  (1,1)  bands  of  the  A^A-X^H  system.^®  Figure  3  shows  DFWM  and  LIF  spectra  at 
0.4  mm  above  the  burner  nozzle.  The  fine-structure  components  of  the  R(8)  transition  in 
Fig.  3  are  distinguished  as  follows:  1  =  2A5/2-2n3/2,  2  = 

doublets  are  labeled  e  and  f  The  A-X  system  is  highly  diagonal,  i.e.,  the  potential  curves 
are  nearly  identical,  causing  the  (0,0)  and  (1,1)  bands  to  have  essentially  equal  transition 
probabilities  with  values  approximately  two  orders  of  magnitude  larger  than  the  off- 
diagonal  transitions  (0,1),  (1,0),  and  (1,2).31>32  goth  the  X  and  A  states  rapidly 
approach  Hund's  case  (b)  coupling  as  rotation  increases.  Therefore,  the  sets  of  Einstein 
emission  and  absorption  coefficients  are  taken  to  be  equal  for  all  of  the  R(8)  transitions. 

HL  QUANTIFICATION  BY  DFWM 

Although  a  generalized  description  of  DFWM  including  collisional,  polarization, 
and  saturation  effects  does  not  exist  at  the  present  time,  Abrams  and  Lind  (AL)3>33,34 
have  presented  a  model  of  DFWM  that  has  given  much  insight  to  the  process,  as  has  been 
shown  by  good  agreement  with  experiments.  1 1 8,20,3 5  jjie  model  considers  a 
nondegenerate  two-level  atomic  system  in  the  presence  of  arbitrary  pump  amplitudes  and 
weak  probe  and  signal  amplitudes.  The  probe  and  signal  amplitudes  are  weak  in  the  sense 
that  they  do  not  significantly  affect  the  level  populations  whereas  the  pump  amplitudes  can 
be  well  above  saturation.  In  this  model  all  fields  have  the  same  polarization  and  their  k- 
vectors  are  configured  in  the  collinear  phase-conjugate  geometry  discussed  above.  For  the 
case  of  equal  intensity  pumps  and  low  absorption,  the  expressions  presented  by  AL  are 
easily  reduced  to  an  expression  for  the  line-center  signal  intensity  given  by 
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FIG.  2.3.  DFWM  and  LIF  spectra  for  the  R(8)  lines  of  the  CH  A^A-X^n  (0,0)  and  (1,1) 
bands  at  a  position  of  0.4  mm  above  the  burner  nozzle.  These  spectra  are  recorded 
simultaneously  and  represent  an  average  of  30  laser  shots.  The  fine-structure  components 
are  distinguished  as  follows:  1  =  ^  the  A  doublets  are 

labeled  e  and  f.  The  small  difference  in  the  Boltzmann  factors  of  the  Ri(8)  and  R2(8) 
transitions  is  enhanced  in  the  DFWM  spectrum  compared  to  the  LIF  spectrum,  because 
the  DFWM  intensities  depend  quadratically  on  the  population  difference  whereas  the  LIF 
intensities  depend  linearly. 
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is  the  line-center  absorption  coefficient,  and 


(3) 


is  the  line-center  saturation  intensity.  Here  AA^^,  =Ng-N,  is  the  population  difference 
(molecules/cm^)  in  the  absence  of  applied  fields,  is  the  square  of  the  transition  dipole 
moment  (proportional  to  absorption  cross-section)  connecting  the  lower  level  (g)  to  the 
upper  level  (e),  r„  and  are  the  population  and  coherence  decay  rates  respectively,  n  is 
the  intensity-dependent  index  of  refraction  (reference  1,  page  222,  Eq.  19),  and  L  is  the 
effective  interaction  length  of  the  radiation  fields.  All  of  the  equations  are  in  SI  units,  and 
the  constants  c,  s^,  and  h  have  their  usual  meanings.  The  model  predicts  that  the  DFWM 
line-center  intensity  will  be  quadratic  in  the  population  difference  of  the  transition. 
Equations  (l)-(3)  are  used  in  this  analysis,  where  we  have  replaced  with  the  ground 
state  population  since  Ng{y",  J’') »  Ng(v\  at  typical  flame  conditions.  The  widths  of  the 
DFWM  lines  for  the  (0,0)  and  (1,1)  bands  are  found  to  be  essentially  the  same.  Thus  we 
take  the  line-center  peak  heights  as  proportional  to  the  integrated  areas. 


IV.  QUANTIFICATION  BY  LIF 

At  high  pressures,  obtaining  rovibrational  distributions  by  LIF  becomes  difficult 
because  of  the  dominating  effect  of  collisions  on  the  fluorescence  quantum  yield. 

However  CH  vibrational  temperatures  can  be  obtained  relatively  free  from  collisional 
effects  because  of  the  highly  diagonal  character  of  the  CH  A^A-X^II  transition  and  the 
similar  energy  transfer  rates  for  the  A2A(v -0,  N')  and  A2A(v -1,  N')  states.36  if  the 
assumption  is  made  that  the  collisional  redistribution  of  population  reaches  steady  state 
over  the  time  period  of  the  laser  excitation,  then  the  following  relation  for  the  fluorescence 
signal  intensity  is  obtained  for  a  two  level  system 
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where  C  is  combined  proportionality  constant  that  includes  the  LIF  collection  efficiency, 
N(  is  the  total  species  population,  and  fg  is  the  Boltzmann  factor  for  the  rovibrational  state 
being  excited.  The  competing  population  transfer  rates  are  defined  as  follows:  is  the 

rate  of  spontaneous  emission;  Q  is  the  collision  quenching  rate;  Fis  the  vibrational  energy 
transfer  rate;  R  is  the  rotational  energy  transfer  rate;  E  is  the  electronic  energy  transfer 
rate;  is  the  rate  of  absorption;  and  is  the  rate  of  stimulated  emission.  The 

assumption  of  steady-state  populations  over  the  duration  (~  20  ns)  of  the  laser  pulse  is  a 
good  one  for  CH  in  an  atmospheric-pressure  flame  because  the  population  redistribution  is 
completed  within  a  couple  of  nanoseconds.^^ 

In  general  the  energy  transfer  rates  in  Eq.  (4)  are  dependent  on  temperature, 
collision  partner  and  quantum  state,  and  therefore,  are  not  known  for  a  wide  range  of 
conditions.  This  dependence  has  precluded  the  application  of  quantitative  LIF  at  high 
pressures  (several  hundred  Torr  and  above)  because  the  detailed  energy  transfer  rates  as  a 
function  of  theses  variables  are  not  known.  However,  Eq.  (4)  can  be  simplified  for  CH 
detection.  The  first  simplification  is  that  only  the  (0,0)  and  (1,1)  bands  need  to  be 
considered  since  they  are  approximately  two  orders  of  magnitude  larger  than  the  off- 
diagonal  (0,1),  (1,0)  and  (1,2)  bands.  In  addition,  the  spontaneous  emission,  stimulated 
emission,  and  absorption  rates  are  the  same  for  the  (0,0)  and  (1,1)  bands.^  ^  Furthermore, 
vibrational  energy  transfer  between  the  A^A(v -0)  and  A2A(v  — 1)  states  is  insignificant  for 
this  system,  however,  rotational  and  electronic  energy  transfer  rates  are  appreciable.^^ 
Since  the  monochromator  bandpass  is  sufficiently  large  to  collect  all  relevant  Q-branch 
transitions  for  both  the  A^A-X^H  (0,0)  and  (1,1)  transitions  while  discriminating  against 
R-  and  P-branch  transitions,  the  detected  LIF  signals  for  R(8)  excitation  will  not  be 
sensitive  to  rotational  redistribution  in  the  A2A(v -0)  and  A2A(v -1)  states.  Therefore,  R 
can  be  eliminated  in  Eq.  (4).  Likewise,  even  if  vibrational  energy  transfer  between  the 
A^A(y'=0)  and  A2A(v  — 1)  states  occurs,  the  monochromator  bandpass  collects  emission 
from  both  levels  with  equal  efficiencies.  Hence  V  can  also  be  eliminated  in  Eq.  (4). 
Electronic  energy  transfer,  on  the  other  hand,  transfers  population  efficiently  fi’om  the  A^A 
to  the  state  with  significantly  different  rates  for  the  A2A(v -0)  and  A2A(v -1)  levels. 
Because  the  B2Z~-X2n  fluorescence  is  outside  of  the  monochromator  bandpass,  the  rate 
of  electronic  energy  transfer  cannot  be  neglected. 

Garland  and  Crosley^^  have  measured  the  ratio  of  electronic  energy  transfer  rate 
to  the  quenching  rate  for  the  A^A(v-0)  and  A2A(v-1)  states  and  report  values  (or  E/Q  of 
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0.025  ±  0.004  and  0. 10  ±  0.02,  respectively.  The  variation  of  E/Q  as  a  function  of  N'  was 
not  statistically  significant  for  these  states.  Finally,  with  =  (IJ^  + 1/  27^  + Eq. 
(4)  reduces  to 


IuF=Cf,iy"  =  (i,J»)N, 


_ 4®) _ 

1.025  Qiy'  =  0)  + 1^\ + (27^  + 1  /  27,  + 1)]  +  ^ 


(5) 


for  LIF  via  the  A^A-X^Il  (0,0)  band,  and 

JuF  =  CUv"  =  1,7")A^, - - T - 

1.10  Q(v' =  1)+ /5[l  +  (27^  + 1  /  27,  + 1)]  + 

for  LIF  via  the  A^A-X^Il  (1,1)  band  where  A  is  the  Einstein  emission  coefficient,  A  = 
^(1,1)  =  ^(0,0),  and  B  is  the  Einstein  absorption  coefficient,  B=B(\,l)=B(0,0). 
Equations  (5)  and  (6)  are  used  to  relate  the  integrated  LIF  signal  intensities  to 
concentrations. 


(6) 


V.  RESULTS  AND  DISCUSSION 


A.  CH  Concentration  Profiles 

Emission  from  the  CH  radical  is  responsible  for  the  blue  coloration  of  hydrocarbon 
flames  and  essentially  defines  the  reaction  zone  of  these  combustion  systems.  This 
emission  originates  primarily  from  chemiluminescence  of  CH*(A2A)  produced  from  the 

reactions26 


(major) 
(major) 
(minor) 

as  opposed  to  thermally  excited  CH.  The  CH  concentration,  however,  is  much  lower 
than  the  CH  concentration,  A[CH*]/A^CH]  ~  lO’^,  and  thus  justifies  the  approximation  in 
section  III  that  AN~Ng(v",  J”).  Although  this  emission  is  abundant  and  easily  detected, 
the  CH  molecules  in  the  A^A  state  are  not  thermalized,  and  rovibrational  distributions  are 
more  indicative  of  the  exothermicity  of  the  above  reactions  rather  than  the  translational 
temperature  of  the  flame.  In  the  flame  studied  here  the  CH  emission  outlines  a  small 
conical  blue  primary  reaction  zone,  of  1.2  mm  diameter  and  1.6  mm  in  height,  and  a 
diffuse  blue-green  secondary  reaction  zone  of  approximately  6  mm  in  height. 


C2  +  0H  =  CH*  +  C0 
C2H  +  O  =  CH*  +  CO 
C2H  +  02  =  CH*  +  C02 
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Relative  CH  ground-state  concentration  information  is  obtained  by  measuring  the 
DFWM  and  LIF  integrated  intensities  of  the  four  R(8)  transitions  as  a  function  of  height  in 
the  flame.  The  Boltzmann  factors  for  the  N"=S  levels  of  the  CH  X^fl  state  are  relatively 
constant  for  temperatures  ranging  fi-om  2300-3300  K,  i.e.,  the  integrated  line  intensities  of 
transitions  originating  fi-om  these  levels  are  weak  functions  of  temperature.  Therefore,  the 
relative  [CH(v"=0,  N'-S)]  concentration  profiles  are  representative  of  the  CH  relative 
concentration  profile  and  can  be  obtained  directly  from  the  integrated  line  intensities  of 
individual  fine-structure  components  of  transitions  originating  from  v"=0,  N"=  8.  From 
Eq.  (1),  the  relative  DFWM  concentration  profile  can  be  obtained  simply  and  is  given  by 

=  ll^(y  =  y^)r' 

Likewise,  the  LIF  concentration  profile  is  obtained  from  Eq.  (5)  and  is  given  by 

N[CH(y)] 

N[CHiy  =  y^^)] 

In  both  Eqs.  (7)  and  (8),  is  the  position  above  the  burner  nozzle  of  maximum  CH 
concentration.  In  addition,  the  DC  level  of  the  LIF  signal  is  recorded  at  each  position 
which  provides  a  measure  of  the  CH*  emission  from  the  flame.  The  emission  is 
normalized  at  the  same  position  as  the  LIF  profile,  y  =  y^  (LIF),  and  is  shown  with  the 
DFWM  and  LIF  concentration  profiles  in  Fig.  4.  Each  data  point  in  Fig.  4  represents  the 
average  of  the  normalized  profile  for  all  four  fine-structure  components.  The  DFWM 
profile  has  been  shifted  by  160  pm  to  longer  distance  (approximately  20%  of  the  beam 
waist)  with  respect  to  the  LIF  profile  to  give  the  best  agreement.  Because  of  the  quadratic 
fall  off  of  DFWM  signal  strength  with  concentration,  the  DFWM  data  points  for  positions 
greater  than  2.5  mm  are  omitted.  The  first  point  in  the  emission  profile  is  higher  than  the 
DFWM  and  LIF  profiles  because  of  laser  light  scattering  from  the  burner  nozzle.  This  is 
the  only  data  point  affected  in  this  manner  because  alt  other  points  are  at  least  one  beam 
waist  from  the  nozzle. 

The  spatial  resolution  of  the  LIF  and  emission  measurements  is  determined  by  the 
magnification  of  the  imaging  system  and  the  monochromator  slits  and  is  50  pm  x  100  pm 
(width  X  height).  The  DFWM  profile  has  a  spatial  resolution  limited  by  the  overlap  region 
of  the  three  beams  each  having  a  fiill-width-at-half-maximum  of  470  pm.  Under  these 
conditions,  both  the  LIF  and  DFWM  signals  spatially  average  over  species  concentration 
and  temperature  gradients,  and  therefore  represent  the  average  CH  concentration  as  a 
function  of  height.  The  distance  between  points  in  the  spatial  profile  is  250  pm  (height)  as 
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FIG.  2.4.  Relative  CH  concentration  profile  obtained  by  DFWM,  LIF  and  CH*  (A^A) 
emission  as  a  function  of  height  above  the  burner  nozzle.  The  CH  concentration  reaches  a 
maximum  in  the  primary  reaction  zone,  region  1,  and  extends  into  the  secondary  reaction 
zone,  region  2,  as  a  result  of  excess  fuel  and  fuel  fragments  reacting  with  diffused  oxygen 
from  the  laboratory  air. 
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measured  by  a  micrometer.  The  DFWM  profile  is  not  significantly  broadened  with  respect 
to  the  LIF  profile  obtained  at  a  spatial  resolution  more  than  2  times  smaller  than  the  step 
size.  This  suggests  that  the  effective  interaction  diameter  of  the  DFWM  beams  is  on  the 
order  of 250  pm. 

There  is  good  agreement  between  the  DFWM,  the  LIF,  and  the  emission  profiles 
which  show  the  CH  concentration  increasing  to  a  maximum  in  the  primary  reaction  zone 
and  then  rapidly  decreasing  at  the  boundary  between  the  primary  and  secondary  reaction 
zones.  This  profile  defines  the  primary  reaction  zone,  region  1,  and  the  small  hump  at 
1.6  mm  in  the  CH  emission  marks  its  boundary.  Since  the  flame  is  fuel  rich,  0  =  2.64  for 
R  =  0.947,  the  CH  concentration  extends  past  the  primary  reaction  zone  to  the  secondary 
reaction  zone,  region  2,  of  the  flame  owing  to  excess  fuel  and  fuel  fi’agments  reacting  with 
oxygen  that  has  diffused  into  the  flame  from  the  surrounding  laboratory  air. 

B.  CH  Vibrational  Temperatures 

In  calculating  the  CH  vibrational  temperatures,  we  take  advantage  of  the  fact  that 
the  (0,0)  and  (1,1)  bands  have  the  same  transition  probabilities,^  ^  and  that  the  ratio  of  the 
rotational  energy  transfer  rate,  R,  to  the  to  the  quenching  rate,  Q,  for  a  given  V'  is  nearly 
the  same  for  the  v -0  and  v -1  levels  of  the  A^A  state.  For  example.  Garland  and 
Crosley^^  report  a  R/Q  value  of  2.4  ±0.1  for  the  N-6,  v -0  level  and  2.8  ±  0.2  for  the 
V-6,  V -1  level  of  the  A^A  state.  These  R/Q  values  suggest  that  the  saturation  properties 
of  the  (0,0)  and  (1,1)  transitions  are  nearly  identical.  Furthermore,  when  specific  fine- 
structure  components  of  these  bands  are  compared,  i.e.,  the  Rie(v"=0)  is  compared  to  the 
Rle(v"=l),  any  effects  arising  fi'om  polarization  will  affect  both  transitions  identically  and 
cancel.  Therefore,  from  Eq.  (1)  we  can  write 

w(R,.wi'’" = ■)  =  i)]^ 

where  Ri.2;e,/  and  F,  correspond  to  the  particular  R-branch  transition  originating  in  a 
specific  spin-orbit  state,  either  Fi(7"=iV"+l/2)  or  F2(/''=Ar"-l/2),  and  A-doublet  level, 
either  e  or  f  Notice  that  since  the  transition  probabilities,  degeneracies,  and  collisional 
widths  are  the  same,  the  DFWM  signal  ratio  depends  exclusively  on  the  square  of  the 
concentration  ratio. 

The  LIF  expressions  can  be  simplified  in  a  like  manner.  The  first  reduction  results 
from  the  fact  that  the  spectra  are  recorded  at  approximately  8  times  the  LIF  saturation 
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intensity,  which  is  discussed  in  Sec.  VC  In  this  case  the  small  electronic  energy  transfer 
rate  corrections  to  the  LIF  intensity  in  the  denominators  of  Eqs.  (5)  and  (6)  can  be 
ignored.  Furthermore,  since  the  rotational  energy  spacing  of  the  v -0  and  v -1  levels  are 
very  similar,  and  the  R/Q  ratios  are  almost  equivalent  for  a  given  N'  in  either  v -0  or  v -1, 
we  make  the  assumption  that  ^(A^Aiv -0)  =  Q(A^Ay=\).  In  this  case,  the  ratio  of  Eq. 
(6)  to  Eq.  (5)  reduces  to 


where  and  have  the  same  definitions  as  in  Eq.  (9).  Note  that  the  LBF  signal 
ratio  depends  only  on  the  concentration  ratio. 

For  an  assumed  Boltzmann  distribution,  a  vibrational  temperature  can  be  obtained 
directly  from  the  ratio  of  the  DFWM  intensities  of  the  same  fine-structure  component  of 
the  (0,0)  and  (1,1)  transitions  by 

Jw(Ri.u/.'’  -i)]  ^ 

~  ®/J 


and  from  the  LIF  intensities  by 


=  exp(-AE„BlkJy), 


(12) 


where  AEy,B  is  the  energy  difference  between  v'-l  and  v"=0  for  a  given  A  component  and 
J",  kg  is  the  BotIzmann  constant,  and  Ty  is  the  vibrational  temperature.  Spectra  like  the 
one  shown  in  Fig.  3  were  recorded  at  the  region  of  maximum  CH  concentration,  and  Eqs. 
(1 1)  and  (12)  were  used  to  obtain  CH  vibrational  temperatures.  An  average  of  six  scans 
gave  a  ratio  of  A^[CH(v"  =  l)]/V[CH(v"  =  0)]  equal  to  0.260  ±0.015  for  DFWM  and 
0.230  ±0.012  for  LIF,  which  correspond  to  vibrational  temperatures  of  2882  ±  123  K  and 
2642  ±  99  K,  respectively.  The  reported  error  bars  represent  one  standard  deviation  and 
are  more  indicative  of  the  precision  of  the  measurements  rather  than  their  accuracy 
considering  the  fact  both  the  DFWM  and  LIF  spectra  spatially  average  over  concentration 
and  temperature  in  the  flame.  This  is  perhaps  the  chief  reason  for  any  difference  between 
the  results,  since  the  DFWM  and  LIF  spectra  represent  spatial  averages  arising  from 
different  sampling  volumes.  The  relative  CH  concentration  measurements  are  not  as 
sensitive  to  this  effect  because  the  comparison  is  between  relative  values  of  the  same 
measurement  (DFWM  or  LIF)  over  the  same  sampling  volume. 
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These  vibrational  temperatures  can  be  used  in  conjunction  with  the  estimated  CH 
concentration  of  30  ppm  to  determine  the  sensitivity.  For  an  average  temperature  of 
2760  K,  30  ppm  corresponds  to  a  CH  number  density  of  8  x  lO^^  molecules/cm^  (8  x 
loll  molecules/cm^  per  quantum  state).  The  data  shown  in  Fig.  3  represents  a  signal-to- 
noise  ratio  approaching  500: 1,  resulting  in  a  sensitivity  limit  (signal-to-noise  ratio  of  1)  for 
CH  of  4x1012  molecules/cm^  (4xl0l®  molecules/cm^  per  quantum  state)  under  the  present 
conditions.  In  reporting  this  number  one  must  consider  the  fact  that  the  actual  flame 
dimensions  limit  the  interaction  length,  L,  of  Eq.  (1)  to  approximately  1  mm.  The 
coherence  length  of  the  laser,  on  the  other  hand,  is  on  the  order  of  a  centimeter,  and  at  2° 
incidence  angle,  the  effective  overlap  length  of  the  beams  is  approximately  1 5  mm.  This 
implies  that  in  an  environment  that  does  not  limit  the  interaction  length,  the  DFWM  signal 
intensity  could  increase  from  the  dependence  by  as  much  as  2  orders  magnitude.  Since 
the  DFWM  signal  increase  quadratically  with  concentration,  an  increase  in  signal  intensity 
by  two  orders  of  magnitude  results  in  only  an  order  of  magnitude  increase  in  sensitivity. 
Therefore,  a  more  reasonable  estimate  of  the  CH  detection  limit  is  4xl0ll  molecules/cm^ 
(4x1 0^  molecules/cm^  per  quantum  state). 

C.  DFWM  and  LIF  Power  Dependence 

The  effect  of  laser  power  on  DFWM  and  LIF  integrated  signal  intensities  for  the 
R(8)  transition  of  the  (0,0)  band  was  determined  for  pump  beam  energies  of  1.3  to  81.3  pJ 
(spectral  intensities  of  0. 17  to  1 1 .2  x  10^  W/cm^-cm’l).  This  range  corresponds  to  total 
energies  of  all  three  beams  from  2.8  to  183  pJ  (spectral  intensities  of  0.39  to 
25.2  X  10^  W/cm^-cm'l).  Figures  5a  and  6a  show  the  DFWM  and  LIF  integrated  signal 
intensities  as  a  function  of  total  laser  energy  (bottom  axis)  and  spectral  intensity  (top  axis). 
Figures  5b  and  6b  are  log-log  plots  for  DFWM  and  LIF,  respectively.  The  integrated 
intensities  in  these  graphs  represent  an  average  of  the  integrated  signal  intensity  of  the 
four-fine  structure  components  weighted  linearly  (LIF)  and  quadratically  (DFWM)  by  the 
relative  Boltzmann  factors  using  an  average  vibrational  temperature  of 2670  K.  Since  the 
term  energies  only  vary  by  4  cm'l  in  1014  cm'l  for  the  Fjg  f(y"=7^"+l/2)  and 
1/2)  components,^^  the  normalization  is  determined  primarily  by  the  relative  degeneracy 
factors  (2J"+1). 

The  DFWM  integrated  signal  intensities  of  Fig.  5a  are  linear  over  a  wide  range  of 
laser  intensities.  The  slope  of  the  line  in  Fig.  5b  resulting  from  a  linear  least  squares  fit  of 
the  data  between  28.1  and  154.7  pJ  (3.9  and  21.3  x  10^  W/cm^-cm"^)  is  0.87  ±  0.07 
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RLG.  2.5.  (a)  DFWM  integrated  signal  intensity  as  a  function  of  total  laser  energy 
(bottom)  and  spectral  intensity  (top),  and  (b)  Log[DFWM  integrated  intensity]  vs. 
Log[total  energy].  The  DFWM  integrated  signal  intensity  is  linear  over  a  large  range  of 
laser  energies.  The  slope  of  the  line  in  (b)  resulting  from  a  linear  least  squares  fit  of  the 
data  between  28.1  and  154.7  jxJ  (3.9  and  21.3  x  10^  W/cm^-cm"^)  is  0.87  ±  0.07 
corresponding  to  an  average  power  dependence  in  this  region  of  Idfwu 
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FIG.  2.6.  (a)  LIF  integrated  signal  intensity  as  a  function  of  total  laser  energy  (bottom) 
and  spectral  intensity  (top),  and  (b)  Log[LIF  integrated  intensity]  vs.  Log[total  energy]. 
The  line  resulting  from  a  linear  fit  to  the  data  points  below  14.1  |iJ  (1.9  x  10^  W/cm^- 
cm'l)  is  shown  in  (a)  as  well  as  a  fit  of  the  data  to  the  power  dependence  predicted  by  Eq. 
(5).  The  LIF  integrated  signal  intensities  in  (a)  show  a  deviation  from  linearity  beginning 
atl4.1  pJ(1.9xl06w/cm2  -cm'l)  resulting  from  saturation.  The  slope  of  the  line 
shown  in  (b)  for  the  data  points  below  14.1  juJ  (1.9  x  10^  W/cm^-cm'l)  is  0.98  +  0.05 
indicating  at  most  only  a  slight  degree  of  saturation  in  this  energy  regime. 


corresponding  to  an  average  power  dependence  in  this  region  of  IcFim  is 

significantly  less  than  the  P  power  dependence  for  unsaturated  DFWM.  The  LIF 
integrated  signal  intensities  in  Fig.  6a  show  a  deviation  from  linearity  beginning  at  14. 1  |iJ 
(1.9  X  10^  W/cm^-cm*!).  The  slope  of  Fig.  6b  for  the  points  below  this  energy  is 
0.98  ±  0.05  indicating  that  the  LIF  signals  are  at  most  only  slightly  saturated  in  this 
intensity  regime.  Bonczyk  and  Shirley^^  have  measured  the  LIF  saturation  intensity  for 
CH  in  an  atmospheric-pressure  oxyacetylene  flame  to  be  0.34  x  10^  W/Cm^-cm'^,  and 
Takubo,  Yano,  Matsuoka,  and  Shimazu^^  measured  0.33  x  10^  W/cm^-cm*!  in  an 
atmospheric-pressure  propane-air  flame.  These  results  are  consistent  with  the  data  of 
Figs.  6a,  and  6b  which  suggest  that  saturation  for  CH  in  an  atmospheric-pressure  flame 
occurs  at  spectral  intensities  on  the  order  of  0.5  x  10^  W/cm^-cm’l. 

Another  indication  of  saturation  is  power  broadening  of  spectral  lines.  Figure  7 
shows  the  line  widths  (fwhm)  for  both  the  DFWM  and  LIF  spectral  lines  as  a  function  of 
total  energy  (bottom)  and  spectral  intensity  (top)  with  each  data  point  representing  an 
average  over  all  four  fine-structure  components.  Both  the  DFWM  and  the  LIF  line  widths 
are  power  dependent  even  at  the  lowest  total  energies  studied  2.8  pJ  (0.39  x  10^  W/cm^- 
cm’l)  In  Fig.  7,  it  is  immediately  apparent  that  the  DFWM  line  wdths  are  substantially 
narrower  than  the  Doppler-broadened  LIF  line  widths  at  low  laser  intensities,  and  that  at 
high  laser  intensities,  the  power-broadened  line  shapes  of  both  techniques  are 
approximately  equal.  The  sub-Doppler  DFWM  line  shape  is  limited  by  the  large  laser 
bandwidth  of  0.18  cm’^  as  opposed  to  the  estimated  collisional  width  (fwhm)  of 
0.02  cm'1.3^  At  total  laser  energies  lower  than  28.1  pJ  (3.9  x  10^  W/cm^-cm'l)  the 
DFWM  line  width  is  smaller  than  the  laser  bandwidth  and  has  a  value  of  0. 12  ±  0.02  cm*l 
at  the  lowest  total  laser  energy. 

Figure  8a  shows  the  LIF  and  DFWM  line  shapes  for  the  Rie(8)  transition  at  a  total 
laser  energy  of  28.1  pJ  (3.9  x  10^  W/cm^-cm'^  ).  The  LIF  line  shape  is  Gaussian  owing  to 
the  large  Gaussian  laser  bandwidth  of  0.18  cm*l  and  Doppler  width  of  0.24  cm"^  (2760  K) 
compared  to  the  collisional  width  of  0.02  cm'^  The  DFWM  line  shape  is  fit  in  Fig.  8a  to  a 
Lorentzian-cubed  line  shape.  Figure  8b  shows  the  best  fit  results  to  the  DFWM  line  shape 
for  a  Gaussian,  Lorentzian,  and  Lorentzian-cubed  profile.  The  Lorentzian-cubed  profile 
gives  the  best  fit  to  the  line  shape  in  agreement  with  the  AL  model  predictions,^>^^ 
although  experimental  conditions  appear  to  be  outside  of  the  model. 


Total  Intensity  (x10®  W/cm^-cm'^) 


Total  Energy  (|liJ) 


FIG.  2.7.  DFWM  and  LIF  line  widths  (fwhm)  as  a  function  of  total  laser  energy  (bottom) 
and  spectral  intensity  (top).  Curves  are  drawn  as  a  guide  to  the  eye.  At  low  energy,  the 
sub-Doppler  DFWM  line  widths  are  significantly  narrower  than  the  Doppler-broadened 
LIF  line  widths,  and  at  high  laser  intensities,  both  techniques  have  comparable  power- 
broadened  line  widths.  The  DFWM  line  width  is  smaller  than  the  laser  bandwidth  (0.18 
cm*^)  for  total  laser  energies  lower  than  28.1  |iJ  (3.9  x  10^  W/cm^-cm'^)  and  has  a 
minimum  value  of  0.12  ±  0.02  cm'^  at  2.8  pJ  (0.39  x  10^  W/cm^-cm'^ ). 


FIG.  2.8.  (a)  DFWM  and  LIF  line  shapes  for  the  Rie(8)  transition  of  the  CH  A^A-X^Il 
(0,0)  band  at  28.1  (3.9  x  10^  W/cm2-cm-l).  The  DFWM  line  shape  is  fit  to  a 

Lorentzian-cubed  profile  as  predicted  by  the  AL  model,  and  the  LIF  line  shape  is  fit  to  a 
Gaussian  profile  owing  to  the  large  Gaussian  widths  (0.18  cm'^  laser  bandwidth  and  0.24 
cm‘1  Doppler  width),  (b)  Comparison  of  the  best  fit  results  to  the  DFWM  line  shape  for 
a  Gaussian,  Lorentzian,  and  Lorentzian-cubed  profile. 
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The  DFWM  integrated  signal  intensities  are  essentially  linear  in  laser  energy  over 
the  range  (two  orders  of  magnitude)  studied.  However,  if  the  line-center  peak  height  as  a 
function  of  pump  laser  intensity  is  plotted,  the  effect  of  saturation  is  more  apparent.  See 
Fig.  9.  The  line-center  signal  intensity  rolls  off  very  rapidly  to  nearly  a  constant  value  as 
predicted  by  Eq  (1).  From  Eq.  (1),  the  DFWM  line-center  power  dependence  is  of  the 
form 
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where  1  is  the  intensity  of  the  forward  pump  field  and  1^  =  1^,=  4/^  in  our  experimental 
configuration.  The  line-center  peak  heights  were  fit  to  Eq.  (13),  and  the  fit  is  also  shown 
in  Fig.  9.  The  fit  gave  an  experimentally  determined  value  of  /“  =  1 1 .7  ±  1 . 1  pJ 
(1.61  ±  0.15  X  106  W/cm2.cm-l  ). 

Contrary  to  naive  first  impressions,  the  sensitivity  of  DFWM  does  not  increase 
without  limit  as  a  function  of  laser  power.  A  measure  of  the  signal-to-noise  ratio  in  these 
experiments  is  the  DFWM  reflectivity,  which  can  be  obtained  by  rearranging  Eq. 

(1)  to  give 

R  dlO  (14) 

Ij,  (l  +  4//4i) 


Equation  (14)  represents  the  ratio  of  the  DFWM  signal  intensity  to  the  input  probe  field 
intensity,  and  this  equation  is  a  good  measure  of  the  signal-to-noise  ratio  because  scattered 
laser  light  fi’om  the  probe  beam  is  the  largest  source  of  noise.  If  only  the  laser  intensity  is 
varied  in  a  series  of  experiments,  as  was  done  for  the  power  dependence  of  Fig.  9,  Eq. 

(14)  can  be  simplified  to  give 


p  _  ^DFim  _  (iiij 


(15) 


where  t]  is  a  combined  constant  that  represents  the  DFWM  signal  intensity  for  a  given  set 
of  conditions  and  can  be  measured  experimentally.  The  DFWM  reflectivity  curve  was 
generated  using  Eq.  (15)  and  the  experimentally  determined  value  of  4^^ ,  and  the  data  of 
Fig.  9  was  divided  by  the  probe  field  intensity.  The  theoretical  and  experimental 
reflectivity  has  been  scaled  to  a  value  of  t]  defined  here  as 
( . 
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FIG.  2.9.  DFWM  line-center  signal  intensity  as  a  function  of  pump  beam  energy  (bottom) 
and  spectral  intensity  (top).  The  error  bars  are  ±  10  %  resulting  primarily  from  the  shot- 
to-shot  laser  energy  fluctuations.  The  line-center  peak  heights  were  fit  to  the  power 
dependence  predicted  by  Eq.  (13).  The  fit  (solid  curve)  gave  an  experimentally  determined 
value  of /“  =  1.7  ±  1.1  pj(1.61  ±  0.15  x  10^  W/cm2-cm-l). 
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FIG.  2.10.  DFWM  reflectivity,  as  a  function  of  //  /„,  for  the  value  of  /“  = 

1.61  ±  0.15  X  10^  W/cm^-cm'l  (1 1.7  ±  1.1  (iJ)  determined  from  the  fit  shown  in  Fig.  9. 
The  reflectivity  has  been  scaled  to  an  experimentally  determined  value  of  rj  =  1 . 14  x  10"^ 
defined  by  Eq.  (16).  The  reflectivity  is  indicative  of  the  signal-to-noise  ratio  for  these 
experiments  and  is  a  maximum  for  laser  intensities  approximately  equal  to  . 
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that  was  measured  at  //  /“  =  1  and  are  plotted  as  a  function  of  //  4^,  in  Fig.  10.  Figure 
10  shows  that  the  maximum  reflectivity,  hence  signal-to-noise  ratio,  is  obtained  at  spectral 
intensities  approximately  equal  to  .  The  highest  signal-to-noise  ratios  also  were 
observed  experimentally  at  these  energies.  As  the  laser  intensity  (energy)  increases  above 
7^,,  the  Une-center  signal  intensity  becomes  constant,  the  peak  profiles  power  broaden, 
and  the  scattered-light-noise-baseline  increases.  The  combined  result  is  that  the  signal-to- 
noise  ratio  decreases.  Working  at  powers  much  greater  than  will  not  improve  the 
detection  limit  of  DFWM,  but  will  only  diminish  the  advantages  of  the  technique,  such  as 
narrow  sub-Doppler  line  shapes. 

VI.  CONCLUSIONS 

In  this  chapter  we  report  the  first  DFWM  study  of  the  CH  radical  and  have 
demonstrated  that  trace  flame  .species  can  be  detected  at  atmospheric  pressure  with  high 
sensitivity  using  DFWM.  Accurate  CH  concentration  profiles  and  vibrational 
temperatures  with  regard  to  LEF  results  have  been  obtained  without  having  to  make 
corrections  for  collisional  processes.  At  atmospheric  pressure  the  DFWM  and  LIF 
intensities  are  comparable;  however,  at  higher  pressures  DFWM  should  prove  to  be  more 
sensitive  than  LIF  because  of  the  effect  of  quenching  collisions  on  the  LIF  signal.  An 
additional  advantage  of  DFWM  is  that  it  generates  a  coherent  signal  beam  that  can  be 
remotely  detected  with  high  efficiency.  This  attribute  enables  nonintrusive  investigation  of 
many  important  chemical  environments  that  suffer  fi'om  limited  optical  access,  high  levels 
of  emitted  radiation,  and  high  pressures,  such  as  plasmas,  flashes,  flames,  and  discharges. 
DFWM  will  be  especially  advantageous  for  the  detection  of  molecular  species  that  do  not 
fluoresce  or  whose  fluorescence  is  complicated  by  collisional  effects. 
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CHAPTER  3 

Theoretical  Development  of  the  Reduction  of  DFWM  Spectra  to 
Relative  Populations  in  the  Weak-Field  Limit 
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"....it  may  be  remarked  that  the  main  object  of  physical  science  is  not 
the  provision  of  pictures,  but  is  the  formulation  of  laws  governing  phenomena 
and  the  application  of  these  laws  to  the  discovery  of  new  phenomena.  If  a 
picture  exists,  so  much  the  better;  but  whether  a  picture  exists  or  not  is  a 
matter  of  only  secondary  importance.  In  the  case  of  atomic  phenomena  no 
picture  can  be  expected  to  exist  in  the  usual  sense  of  the  word  'picture*  by 
which  is  meant  a  model  fimctioning  essentially  on  classical  lines.  One  may, 
however,  extend  the  word  'picture'  to  include  any  way  of  looking  at  the 
fundamental  laws  which  makes  their  self-consistency  obvious.  With  this 
extension,  one  may  gradually  acquire  a  picture  of  atomic  phenomena  by 
becoming  familiar  with  the  laws  of  quantum  theory." 

P.  A.  M.  Dirac 
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I.  INTRODUCTION 


LIF  is  a  powerful  spectroscopic  technique  for  analyzing  molecular  species  in 
experiments  that  range  from  chemical  reaction  dynamics  to  plasma  physics.  A  common 
goal  of  these  experiments  is  to  extract  accurate  values  for  relative  populations  from  the 
relative  intensity  distribution  of  the  emitted  light.  Several  authors  have  shown  how  the 
emitted  light  intensity  depends  on  the  excitation/detection  geometry  (polarization),  the 
temperature,  pressure,  and  chemical  composition  of  the  experiment  (collisional 
relaxation),^’^  and  the  laser  power  employed  (saturation).^"^  In  general,  the  population 
as  well  as  the  elements  of  the  first-  and  second-rank  multipole  moments  of  the  total 
angular  momentum,  called  the  orientation  and  the  alignment,  respectively,  can  be 
determined  from  low-pressure  experiments  performed  under  single-collision  conditions. 

In  collisionally  dominated  environments,  however,  the  determination  of  these  quantities 
becomes  extremely  difficult. 

An  appealing  alternative  approach  to  the  study  of  collisionally  dominated 
environments  is  DFWM.^®  DFWM  uniquely  provides  information  regarding  the  chemical 
composition,  kinetics,  and  dynamics  of  environments  that  are  inherently  difficult  to  study 
because  of  radiant  interference  from  the  emission  of  excited  species.  ^  ^  One  of  the  most 
important  features  of  DFWM  is  the  large  four-wave  mixing  enhancement  exhibited  when 
the  laser  frequency  is  tuned  to  a  molecular  one-photon  resonance  (this  includes  electronic, 
vibrational,  and  rotational  transitions);  this  feature  makes  DFWM  a  very  sensitive 
molecular  probe.  Other  important  features  include  sub-Doppler  spectral  resolution, 
excellent  spatial  and  temporal  resolution,  imaging  capabilities,  nonintrusive  detection,  and 
remote  sensing.  Because  of  these  attributes,  the  potential  of  DFWM  spectroscopy  to 
obtain  chemical  and  dynamical  information  has  generated  much  excitement.  Like  LIF, 
however,  DFWM  requires  an  understanding  of  polarization,  collisional,  and  saturation 
effects  before  information  can  be  extracted  from  the  signal  intensities. 

At  present,  DFWM  signal  intensities  are  most  commonly  interpreted  by  applying 
the  stationary  absorber  model  proposed  by  Abrams  and  Lind,  which  considers  a 
nondegenerate  two-level  system  in  the  presence  of  arbitrarily  strong  pump  fields.  This 
model  accounts  only  for  population  contributions  to  the  DFWM  response  and  neglects 
contributions  from  orientation  and  alignment.  This  model's  predictions  have  been  applied 
with  much  success  to  relate  saturated  DFWM  signal  intensities  to  populations  in 
experiments  performed  under  collisionally  dominated  conditions.!  Note,  however, 

that  signals  from  a  single  rotational  branch  and  polarization  configuration  were  used  in 


these  experiments.  These  conditions  serve  to  minimize  the  error  introduced  from  applying 
this  model  to  molecular  systems  with  degenerate  magnetic  sublevels. 

To  date,  attempts  to  formulate  a  complete  model  that  includes  polarization, 
collisional,  and  saturation  effects  have  been  unsuccessful.  Studies  that  emphasize 
saturation  effects  ^  ^"27  are  limited  to  level  degeneracies  in  the  range  of  1  to  3.  Even 
numerical  approaches'^  are  computationally  impractical  at  present  for  level  degeneracies 
typical  of  molecular  species.  Many  authors^®, 29-39  h^ve  thus  resorted  to  perturbative 
treatments  (no  saturation)  to  gain  insight  into  the  DFWM  response  of  degenerate  systems. 

The  aim  of  this  chapter  is  to  present  a  theory  of  DFWM  that  can  be  used  to 
interpret  the  spectra  of  molecular  species.  This  goal  is  accomplished  by  deriving 
expressions  via  time-independent  diagrammatic  perturbation  theory  that  account  for  the 
DFWM  polarization,  collisional,  and  velocity  effects  in  the  weak-field  limit.  In  our 
treatment,  we  assume  that  the  DFWM  process  couples  levels  of  sharp  (definite)  angular 
momentum  J  (omitting  nuclear  spin).  We  consider  three  input  fields  of  arbitrary 
polarization  that  interact  with  an  isotropic  sample  to  produce  a  fourth  field.  Therefore, 
our  treatment  is  general  in  that  it  applies  to  molecular  species  for  which  J  is  a  good 
quantum  number.  The  general  result  is  specialized  to  apply  to  circularly  and  linearly 
polarized  fields  that  interact  in  nearly  collinear  phase  matching  geometries  in  collisional 
environments  where  the  multipole  moments  of  the  total  angular  momentum  distribution 
relax  independently  (isotropic  relaxation)  and  at  the  same  rate.  Figure  1  shows  the  three 
specific  geometries  we  treat.  In  the  following  chapter,  we  extend  the  key  results 
presented  here  to  interpret  data  taken  under  saturated  conditions. 

The  remainder  of  this  chapter  is  organized  into  four  sections  and  one  appendix.  In 
Sec.  II  we  present  expressions  for  the  DFWM  signal  intensity  as  a  function  of  input 
polarization,  collisional  dephasing,  and  experimental  geometry  under  the  assumption  that 
the  rates  of  collisional  relaxation  of  the  population,  the  orientation,  and  alignment  are  the 
same.  This  section  is  intended  as  a  guide  for  the  experimentalist.  In  Sec.  Ill  we  interpret 
the  DFWM  signal  intensity,  as  has  been  done  for  LIF,  in  terms  of  the  population,  the 
orientation,  and  the  alignment.  In  Sec.  IV  we  compare  our  results  with  those  of  other 
treatments,  and  in  Sec.  V  we  present  conclusions  of  our  findings.  In  the  Appendix 
expressions  are  presented  for  the  polarization  tensor  products  necessary  to  extend  this 
treatment  to  other  experimental  configurations. 


n.  DEPENDENCE  OF  DFWM  SIGNAL  INTENSITIES  ON  POLARIZATION, 
COLLISIONS,  AND  PHASE  MATCHING  GEOMETRY 

A.  General  Expressions  and  Discussion 

In  four-wave  mixing  three  incoming  waves  with  electric  fields  E,  (r,r),  E^{r,t), 
and  Ej(r,t),  propagation  vectors  k^,  k^,  and  k^,  and  frequencies  ©i,  ©2>  and  ©3  interact 
through  the  third-order  nonlinear  susceptibility  to  generate  a  fourth  field,  E^,  with 
propagation  vector  k^  and  frequency  ©4.  The  electric  fields  are  defined  as 

Ejirj)  =  +C.C.  ;  Ej  (1) 

where  Ej  is  the  vector  amplitude,  Cj  is  the  scalar  amplitude,  and  is  the  normalized 
(Sj'Bj  =1)  polarization  unit  vector  of  the  electric  field  labeled  j.  The  conventions  for 
expressing  the  unit  vectors  are  given  in  the  Appendix. 

For  fully  resonant  DFWM  ©,=©2=©3=©4=©,  and  we  assume  that  the 
excitation  bandwidth  is  sufficiently  narrow  compared  with  the  density  of  states  (including 
Doppler  broadening)  of  the  absorbing  molecules  so  that  the  interaction  is  exclusively 
between  the  degenerate  magnetic  sublevels  of  the  two  levels  involved  in  the  one-photon 
resonant  transition.  These  levels  are  assumed  to  be  characterized  by  total  angular 
momentum  J.  Figure  2  shows  a  schematic  energy-level  diagram  of  the  process  described 
above.  In  this  figure  we  restrict  the  letters  g  and  g'  to  refer  to  two  degenerate  magnetic 
sublevels  (either  the  same  or  different)  of  the  lower  level  (usually  the  ground  electronic 
state),  and  the  letters  e  and  e'  refer  to  two  degenerate  magnetic  sublevels  (either  the  same 
or  different)  of  the  upper  level  (usually  another  electronic  state).  In  the  following 
discussion  we  adopt  the  standard  nomenclature  for  the  four  fields  involved  in  four-wave 
mixing:  1  and  3  refer  to  the  two  pump  fields,  2  refers  to  the  probe  field,  and  4  refers  to  the 
signal  field. 

The  general  expression  for  X^^^  that  applies  here  consists  of  96  terms  (48  for  the 
upper  level  and  48  for  the  lower  level)  that  completely  describe  the  interaction  of  the 
electric  fields  with  the  molecular  system. ^^>^1  These  terms  differ  in  the  time  ordering  of 
the  interaction  fields  and  the  permutations  of  the  quantum  states  involved.  Evaluating  all 
of  these  terms  is  a  formidable  task  in  general;  however,  in  a  fiilly  resonant  four-wave 
mixing  experiment,  the  intensity  at  a  particular  resonance  in  the  spectrum  is  dominated  by 
only  sixteen  terms..  This  result  is  well  established  in  the  literature  on  resonant  coherent 
CARS.^^»^^  These  terms  are  represented  using  double-sided  Feynman  diagrams  in  Fig.  3. 
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FIG.  3.2.  Energy-level  diagram  for  DFWM  of  a  degenerate  two-level  system.  In  the 
figure  i  =  l,j  =  3  or  i  =  3J  -  1. 
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FIG.  3.3.  Double-sided  Feynman  diagrams  used  to  construct  the  DFWM  third-order 
susceptibility.  The  first  set  of  eight  resonant  contributions  is  formed  by  substituting  i  =  1, 
7  =  3  and  the  second  by  i  =  3, 7  =  1.  The  energy-resonant  denominators  (Eq.  (23)) 
associated  with  each  set  of  diagrams  is  indicated. 
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For  fully  resonant  DFWM  the  summation  over  all  state  permutations  reduces  to  a 
summation  over  the  magnetic  sublevels  Mg',  Mg,  and  Mg>  of  the  two  levels  involved. 

Furthermore,  if  the  initial  magnetic  sublevel  distribution  is  isotropic,  i.e.,  all  degenerate  M 
levels  are  equally  populated  and  no  phase  relation  exists  between  them,  the  summation  is 
equally  weighted.  Taking  these  aspects  into  consideration,  we  write  the  DFWM  signal 
intensity  as 


where 


£(o) ) = [zf,  (® ) + z|,  (ti> ) + 1;,  (0) ) + z;,  (ffl )] 


(3) 


is  the  total  line  shape  function  and 


is  the  total  geometric  factor.  In  Eq.  (2)  Ij  is  the  intensity  of  the  electric  field  labeled  j,  Ng 
and  Ng  are  the  total  populations  of  the  levels  g  and  e,  respectively,  in  the  absence  of 
applied  fields,  and  Bg^{jg,Jg)  is  the  Einstein  absorption  coefficient  that  connects  the  level 
with  total  angular  momentum  Jg  to  the  level  with  Jg.  In  Eq.  (3)  the  )  are  complex 
line  shape  functions.  In  Eq.  (4)  the  Gp(e4,e,,  8^,82 are  geometric  factors  and 
depend  solely  on  the  polarization  unit  vectors  of  the  electric  fields  Sj  and  the  total  angular 
momentum  quantum  numbers  Jg  and  Jg.  The  dimensionless  weighting  factors  )  and 
)  of  Eq.  (4)  are  defined  as 


£((0) 


(5a) 


and 


(5b) 


and  they  satisfy  the  condition 


(6) 


Equation  (2)  is  the  key  result  of  this  section  and  expresses  the  DFWM  signal  intensity  as  a 
product  of  a  concentration  part,  AN^,  a  one-photon  molecular  part,  a  line  shape  part, 
|£(0  )|  ,  and  a  laboratory-frame  geometric  part,  |Gp  .  These  quantities  relate  to  different 
aspects  of  the  experiment,  and  the  DFWM  signal  intensity  is  directly  proportional  to  them. 
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The  Einstein  absorption  coefficients  describe  the  strength  of  the  interaction  of  the 
molecule  with  the  excitation  fields  in  the  molecular-frame.  These  coefficients  can  be 
related  to  other  molecular  parameters  such  as  the  absorption  cross-section,  the  oscillator 
strength,  the  line  strength,  the  spontaneous  emission  lifetime,  and  the  transition-dipole 
moment.^2  7^5  literature  provides  an  extensive  theoretical  and  experimental  data  base  of 
these  parameters.  The  physical  interpretation  of  the  geometric  factors 
Gf  (54, 81,8^  ,82;^^, y*)  is  given  in  Sec.  III.  The  crucial  point  here  is  that  these  factors  are 
just  real  numbers  that  account  for  the  geometry  of  the  interaction  of  the  molecule  with  the 
electric  fields  in  the  laboratory  frame.  These  terms  are  given  in  Table  3. 1  for  the  allowed 
polarization  configurations  of  circularly  and  linearly  polarized  light  for  the  phase  matching 
geometries  of  Fig.  1.  The  angle  0  in  Fig.  1  is  small,  usually  less  than  2°,  so  we  assume 
that  all  of  the  electric  fields  involved  in  the  DFWM  process  (excitation  and  signal) 
propagate  along  the  space-fixed  Z-axis.  The  collinear  beam  approximation  in  turn  restricts 
the  electric  fields  to  lie  in  the  space-fixed  XY-plane,  and  we  use  the  following 
conventions:  (R)  for  right  circularly  polarized,  (L)  for  left  circularly  polarized,  X  for  (j)  =  0 
linearly  polarized,  and  Y  for  (j)  =  7t/2  linearly  polarized  (Fig.  4).  Finally,  the  line  shape 
functions  Z,"^((o  )  are  defined  as 


f  /(v)^'v 

J  [(*y-*j)  v-/r,] 


»  (7) 


where  /(v)  is  the  normalized  (Maxwell-Boltzmann)  velocity-distribution  function,  is 
the  energy  difference  between  the  e  and  g  levels,  and  is  the  homogeneous  full-width 
(no  Doppler  broadening)  of  a  dipolar  transition  between  the  e  and  g  levels.  The  sub-  and 
superscripts  on  L  refer  to  the  first  terms  of  Eq.  (7).  Specifically  for  the  line  shape  function 
^"2(0  ),  w  is  the  quantum-state  label  of  r„  where  l/r„  is  the  total  lifetime  of  the  rt*^  level, 
and  j  is  the  field  label,  either  1  or  3,  that  corresponds  to  the  k  vector  difference 
^j2  -  (.^j  The  velocity  integration  of  Eq.  (7)  has  been  performed  by  many 
authors  ^^>29-3 1,40,43  gnj  jg  jjot  discussed  here. 

Only  collisional  relaxation  is  considered  (no  spontaneous  emission)  in  Eq.  (7),  and 
we  assume  that  the  collisional  relaxation  of  the  population,  the  orientation,  and  the 
alignment  can  be  represented  by  a  single  rate  r„.  This  "single  relaxation"  assumption  is 
generally  valid;  however  in  experiments  where  efficient  energy  transfer  collision  partners 
such  as  water  are  absent,  levels  characterized  by  low  values  of  J  may  relax  with  multiple 
rates.  This  topic  is  discussed  in  Sec.  IIIC. 
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Table  3.1.  Geometric  factors  for  circular  and  linear  polarization  combinations  with  the  following  conventions:  (R)  for  right 
circular  polarization;  (L)  for  left  circular  polarization;  X  for  (|)  =  0  linear  polarization;  and,  Y  for  ([)  =  7t/2  linear  polarization. 
Permuting  R  L  or  X  Y  does  not  change  the  values  of  the  geometric  factors. 
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In  Eq.  (4),  the  geometric  factors  are  weighted  by  the  line  shape  fiinctions  defined 
in  Eq.  (7).  Therefore,  the  total  DFWM  polarization  dependence  for  a  given  experiment  is 
sensitive  to  more  than  the  polarizations  of  the  excitation  fields;  it  is  also  sensitive  to  the 
phase  matching  geometry,  the  velocity  distribution  of  the  absorbing  molecules,  and  the 
collisional  dynamics  caused  by  the  environment.  This  interdependence  is  an  important 
aspect  of  DFWM. 

For  many  cases  the  two  geometric  factors  of  Eq.  (4)  are  equal  or  equally  weighted. 
In  these  situations  the  total  geometric  factor  becomes 

and  is  independent  of  the  velocity  distribution  of  the  absorbing  molecules  and  the 
collisional  dynamics  caused  by  the  environment.  All  of  this  information  is  then  contained 
exclusively  in  the  total  line  shape  function  £(co  )!  This  is  a  powerful  result  because  the 
DFWM  polarization  dependence  has  been  disentangled  from  the  collisional  dynamics  of 
the  experimental  environment  and  depends  only  on  the  polarization  states  of  the  incident 
fields  and  the  total  angular  momentum  of  the  levels  involved.  Equation  (8)  has  been 
evaluated  for  all  allowed  polarization  configurations  for  circularly  and  linearly  polarized 
light,  and  the  results  are  given  in  Table  3.2  as  a  function  of  transition  type  using  the  same 
conventions  as  in  Table  3.1. 

The  total  geometric  factors  of  Table  3.2  can  be  related  to  what  is  more  commonly 
referred  to  as  four-wave  mixing  line  strength  factors 
by 

where  Sjj^  is  the  molecular  rotational  line  strength  (Honl-London  factor)  for  a  one- 
photon  transition^  and  Gl{z^,z^,z^,Z2,J g,J e)  is  defined  in  Eq.  (8).  In  the  general  case, 
population  distributions  cannot  be  extracted  using  the  line  strength  factors  of  Eq.  (9), 
because  the  total  geometric  factor  depends  on  the  line  shape  functions,  i.e.,  the  overall 
DFWM  polarization  response  depends  on  more  than  the  polarization  of  the  excitation 
fields. 

In  our  discussion  of  polarization  effects,  it  is  useful  to  derive  expressions  that 
relate  how  observed  DFWM  signals  vary  for  a  given  transition  as  the  polarization  of  the 
excitation  fields  is  varied.  This  variation  in  DFWM  signals  is  best  expressed  in  terms  of  a 
polarization  ratio  defined  as 
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Table  3.2.  Total  geometric  factors  for  FB-DFWM  and  for  BB-.  PC-DFWM  when  Chapter  3  Eq.  (8)  is  valid.  The  conventions 
are  the  same  as  Table  3.1,  and  permuting  R  o  L  or  X  Y  does  not  change  the  values  of  the  total  geometric  factors. 
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which  is  merely  the  ratio  of  the  DFWM  signal  intensities  for  two  different  polarization 
configurations  in  which  all  other  aspects  of  the  experiment  are  the  same.  The  three 
fundamental  polarization  ratios  for  circularly  and  linearly  polarized  light  that  can  be 
derived  using  the  values  of  Table  3.2  in  Eq.  (10)  are  given  in  Table  3.3.  Note  that  more 
ratios  can  be  formed  by  permuting  o  Z  and  X  <->  Y,  but  they  are  identical  to  those 
given  in  Table  3.3.  The  polarization  ratios  rapidly  approach  a  high-Z limit,  and  therefore, 
these  values  generally  describe  the  DFWM  signal  variation,  as  the  polarization 
configuration  is  changed  for  those  systems  for  which  Eq.  (8)  is  valid. 

We  turn  to  a  discussion  of  specific  experimental  conditions  in  which  we  point  out 
the  cases  for  which  Eq.  (4)  reduces  to  Eq.  (8).  Here,  we  invoke  the  grating  picture  of 
DFAVM  to  gain  some  insight  as  to  the  relative  magnitude  of  the  line  shape  functions 
L"j2  (co  ) .  The  grating  picture  of  DFWM  involves  the  interference  of  the  excitation  fields. 
When  the  excitation  fields  are  tuned  to  a  one-photon  resonance,  the  interference  between 
two  of  the  fields,  E2  and  Ej  (j=l  or  3),  leads  to  a  spatial  modulation  in  the  (complex) 
refi-active  index  and  hence  forms  a  grating  that  can  scatter  the  third  field  into  a  fourth  field 
(DFWM  signal).  grating  is  characterized  by  the  grating  it- vector, 

=  (kj  -  ^2  )>  which  describes  the  grating's  orientation  and  spacing  in  the  laboratory 
frame.  The  grating  spacing  Dji  is  given  by 
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2sin(0/2) 


(11) 


where  |AA:y2|*s  the  magnitude  of  the  grating  ^-vector  and  0  is  the  angle  between  the  two 
propagation  vectors,  kj  and  k^,  of  the  fields  forming  the  grating  (0  <  0  <  7t).  When  0  is 
small  the  grating  spacing  is  large,  and  when  0  is  large  the  grating  spacing  is  small.  If  the 
absorbers  involved  in  forming  the  grating  are  moving,  such  as  in  a  gas,  the  grating  with  the 
smallest  spacing  is  "washed  out"  to  a  higher  degree  when  Doppler  broadening  is 
significant.45  This  argument  is  useful  in  predicting  the  general  DFWM  polarization 
dependence  and  is  substantiated  by  evaluating  the  integrals  of  Eq.  (7).43 
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B.  Phase  Matching  Geometry  1:  Forward  Box  DFWM  (FB-DFWM) 


We  are  prepared  to  consider  the  polarization  effects  for  specific  phase  matching 
geometries  and  the  relative  sensitivity  of  these  effects  to  the  collisional  environment.  For 
the  FB-DFWM  geometry  of  Fig.  1,  all  fields  are  nearly  copropagating,  0  s  0,  and 
1^121  =  1^321  =  0.  Therefore  all  the  gratings  formed  in  the  FB-DFWM  phase  matching 
geometry  have  large  grating  spacings,  =  A2 

A(cD)  =  A(a))  ;  A((o)  =  A(cd).  (12) 

Substituting  Eq.  (12)  into  Eqs.  (5a)  and  (5b)  shows  that  geometric  factors  that 
correspond  to  the  interchange  of  the  two  pump  fields  (1  and  3)  are  equally  weighted.  The 
total  geometric  factor  has  the  simple  form  of  Eq.  (8)  for  all  transition  types  and 
polarization  configurations.  Therefore  Tables  3.2  and  3.3  contain  all  of  the  relevant 
polarization  information  for  FB-DFWM.  As  shown  in  Table  3.3,  the  polarization  ratios 
are  equal  for  any  entry  with  orthogonally  polarized  pump  fields  (8j  •  63  =  0).  This  result 
distinguishes  the  FB-DFWM  configuration  fi'om  the  other  configurations  as  we  shall  see 
below.  This  distinction  is  in  addition  to  the  fact  that  FB-DFWM  is  not  a  sub-Doppler 
spectroscopic  technique. 

C.  Phase  Matching  Geometries  2  and  3:  Backward  Box  and  Phase  Conjugate 
DFWM  (BB-  and  PC-DFWM) 

For  the  BB-DFWM  and  PC-DFWM  geometries  of  Fig.  1,  fields  1  and  2  are  nearly 
copropagating,  |M,2|  =  0,  and  fields  3  and  2  are  nearly  counterpropagating, 

I  Afc32 1  =  2^  =  471  /  A, .  Therefore  two  types  of  gratings  are  formed  in  these  phase  matching 
geometries:  one  has  a  large  grating  spacing,  D,2  ^  00;  and  the  other  has  a  small  spacing, 
D32  ^  A,  /  2 .  Allowing  for  more  wash  out  of  the  small-spaced  grating  compared  to  the 
large-spaced  grating  when  the  molecules  are  moving,  we  may  write  the  following 
inequalities; 

A(a))>Z|2(0)  ;  A(a))^A(co).  (13) 

Unlike  for  FB-DFWM,  the  total  geometric  factors  for  BB-  and  PC-DFWM  that 
correspond  to  the  interchange  of  the  two  pump  fields  (1  and  3)  are  generally  not  weighted 
equally.  Hence  the  overall  DFWM  polarization  response,  i.e.,  the  total  geometric  factor 

of  Eq.  (4),  will  be  sensitive  to  the  collisional  environment 
because  the  weighting  factors  11^3(0) )  and  IF3,  (co  )  are  dependent  on  the  relative  relaxation 
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rates  of  the  g  and  e  levels.  This  dependence  complicates  the  interpretation  of  DFWM 
spectra  for  these  phase  matching  geometries,  and  the  weighted  average  expressed  in  Eq. 
(4)  generally  must  be  performed. 

In  some  cases,  however,  the  polarization  dependence  for  BB-  and  PC-DFWM  can 
be  disentangled  from  collisional  effects.  These  cases  arise  when 

or  when 

^g=^e  =>  )  =  A*2  (®  )  5  )  =  ^2(®  )  •  (14b) 

For  these  cases  Eq.  (4)  reduces  to  Eq.  (8),  and  again,  the  DFWM  polarization  response 
only  depends  on  the  polarization  states  of  the  incident  fields  and  the  total  angular 
momentum  of  the  levels  involved.  Equation  (14a)  is  satisfied  for  all  Q-branch  (A/  =  0) 
transitions  and  for  P-branch  (A/  =  -1)  and  R-branch  (A/  =  +1)  transitions  when  the  two 
pump  fields  are  of  the  same  polarization  (e,  •  Ej  =  1).  Therefore,  the  simple  results  of 
Tables  3.2  and  3.3  can  be  used  to  evaluate  BB-  and  PC-DFWM  signal  intensities  for  Q- 
branch  transitions  and  for  P-  and  R-branch  transitions  when  e,  •  83  =  1  in  addition  to  FB- 
DFWM  signal  intensities.  For  P-  and  R-branch  transitions  with  orthogonally  polarized 
pump  fields  (e,  •  63  =  0),  however,  Eq.  (14b)  must  be  satisfied  for  the  results  of  Tables  3.2 
and  3.3  to  apply,  i.e.,  the  upper  and  lower  levels  must  relax  at  the  same  rate.  If  neither 
Eq.  (14a)  nor  (14b)  is  satisfied,  the  DFWM  polarization  response  depends  on  the 
dynamics  of  the  collisional  environment  under  study.  This  drawback  is  offset  because 
collisional  relaxation  information  can  be  obtained  from  polarization-ratio  measurements 
without  having  to  resolve  spectral  line  shapes.  ^  ^>46 

Grating  wash  out  leading  to  the  inequality  in  Eq.  (13)  often  occurs  for  Doppler- 
broadened  systems  in  which  the  homogeneous  fiill-width  of  the  spectroscopic 
transition  between  levels  e  and  g  is  smaller  than  the  Doppler  width  Aco  ^  (2r^g  <  A®  ^). 
Experiments  conducted  under  low-pressure,  high-temperature  conditions,  such  as 
atmospheric-pressure  (and  lower)  flames  and  plasmas,  are  examples  of  Doppler-broadened 
systems.  Again  invoking  the  grating  picture  of  DFWM  for  this  case  yields  Zfj  >  Zfj  and 
1^2  >1^2  ,  ^  other  words,  the  small  spaced-grating  is  washed  out  more  than  the  large¬ 
spaced  grating.  Substituting  these  inequalities  into  Eq.  (4)  shows  that  P-  and  R-branch 
transitions  with  orthogonally  polarized  pump  fields  are  sensitive  to  differences  in  the 
relaxation  rates  of  the  g  and  e  levels.  To  emphasize  this  point  consider  the  case  when 
Tg  <  r^.  Given  two  orthogonal  polarization  states  e'  and  e",  we  have 


(15a) 


for  R-branch  (AJ  =  +1)  transitions.  If  >  T,,  the  inequalities  of  Eqs.  (15a)  and  (15b)  are 
reversed.  These  inequalities  hold  for  frequency-integrated  DFWM  signal  intensities  as 
well.  Analytic  expressions  for  the  line  shape  functions  in  the  Doppler-broadened  limit  are 
found  in  references  10  and  43  and  can  be  used  to  quantify  the  inequalities  of  Eqs.  (15a) 
and  (15b). 

Recall  that  for  P-  and  R-branch  transitions  with  the  total  geometric  factor 

is  invariant  to  the  exchange  of  the  polarizations  of  the  pump  fields.  This  result  is 
important  because  even  though  most  flame  and  plasma  experiments  are  Doppler- 
broadened,  diagnostic  applications  typically  involve  the  investigation  of  molecular  species 
in  which  rotational  energy  transfer  collisions  are  the  dominant  type  of  collision. 
Furthermore  the  energy-level  spacings  of  the  g  and  e  levels  of  these  species  are  usually 
similar,  and  as  a  result,  they  have  similar  rotational  energy  transfer  rates.  For  these 
experiments,  e  F^,  and  to  a  first  approximation,  the  simple  results  of  Tables  3.2  and  3.3 
describe  the  DFWM  polarization  response.  This  response  is  discussed  in  detail  in  Chapter 
4  for  the  case  of  the  CH  radical  in  an  atmospheric-pressure  flame. 

We  conclude  this  discussion  by  considering  systems  in  the  homogeneously 
broadened  limit,  i.e.,  2T^  »A(0p.  Experiments  performed  in  supercritical  water^^  that 
are  characterized  by  high  pressure  (tens  of  atmospheres)  and  low  temperature  (~  400  K) 
closely  approximate  this  limit.  In  this  limit  and  Z;^  =  Z^j  because  the  absorbers 

are  effectively  relaxed  before  they  have  a  chance  to  move,  and  thus  neither  the  large-  nor 
small-spaced  gratings  are  washed  out.  Therefore  for  homogeneously  broadened  systems 
Eq.  (4)  reduces  to  Eq.  (8),  and  the  results  presented  in  Tables  3.2  and  3.3  apply.  For  a 
more  precise  evaluation  of  Eq.  (4),  Wandzura'^^  obtained  analytic  expressions  for  the  line 
shape  functions  in  this  limit  that  can  be  used  to  calculate  the  DFWM  polarization 
dependence. 


D.  Other  Signal  Contributions 


The  discussion  above  assumes  that  no  other  process  contributes  to  the  DFWM 
signal.  For  some  experimental  conditions,  however,  other  types  of  laser-induced  grating 
phenomena^^  can  coherently  scatter  light  along  the  DFWM  signal  direction.  The  most 
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likely  candidate  is  a  thermal  grating  that  results  from  the  localized  absorption  of  the  laser 
energy  and  subsequent  heating  of  the  medium.  This  absorption  leads  to  a  spatial-density 
modulation  in  gases  that  can  produce  signals  on  the  order  of  or  greater  than  the  DFWM 
response  under  some  conditions.^^>^^>^^  Another  intensity  grating  effect  observed  at  high 
laser  intensities  is  an  electrostrictive  grating^l  that  results  from  an  electric-field-induced 
modulation  in  the  gas  density.  These  types  of  effects,  however,  can  result  only  from  the 
interference  of  two  fields  that  have  the  same  polarization  (intensity  gratings);  fields  with 
orthogonal  polarizations  do  not  produce  such  effects  (polarization  gratings)  because  there 
is  no  spatial  modulation  of  the  field  intensity  (only  its  polarization  is  spatially 
modulated).^^’^^  If  the  molecule  is  optically  active  (chiral),  however,  then  circular- 
dichroism  induced  thermal  gratings  can  be  formed  with  orthogonally  polarized  fields. 

If  other  types  of  gratings  are  present,  polarization  dependences  different  from 
those  discussed  above  are  observed.  In  particular,  polarization  configurations  that  involve 
large-spaced  intensity  gratings  have  anomalous  signal  intensities  if  additional  intensity 
grating  contributions  are  present.  A  good  way  to  test  for  this  effect  is  to  measure  the 
polarization  ratio  of  a  Q-branch  transition  with  orthogonally  polarized  pump  fields  in  the 
BB-  or.  PC-DFWM  phase  matching  geometry.  If  the  polarization  ratio  is  not  unity, 
another  process  is  likely  to  be  contributing.  Rahn  and  Brown^^  used  similar 
polarization  techniques  to  determine  that  thermal  gratings  make  significant  contributions 
to  OH  spectra  taken  in  atmospheric-pressure  flames.  If  other  gratings  are  determined  to 
be  contributing,  the  YXXY  polarization  configuration  best  discriminates  against  their 
contribution  for  all  experimental  conditions,  but  the  YYXX  and  RRLL  configurations  are 
also  good  choices  for  Doppler-broadened  systems  with  grating  wash  out  because  only 
small-spaced  intensity  gratings  are  present. 

£.  Influence  of  Hyperfine  Structure 

A  full  treatment  considering  the  influence  of  hyperfine  structure  is  possible,^ 
but  the  additional  complexity  required  to  include  such  effects  is  not  warranted^.^S  for  the 
conditions  described  in  this  chapter,  i.e.,  the  spectroscopy  of  molecular  systems  in 
collisionally  dominated  environments.  This  conclusion  can  be  justified  by  invoking  the 
vector  model  of  angular  momenta.  In  the  absence  of  collisions,  the  rotational  angular 
momentum  J  couples  to  the  nuclear  angular  momentum  I  to  form  the  total  angular 
momentum  F  so  that  the  prepared  direction  of  J  in  the  laboratory-frame  is  lost  to  some 
degree  (depolarized).  This  directional  blurring  obviously  influences  the  geometric  factors 


related  to  the  distribution  of  J  in  the  laboratory  frame.  If  the  magnitude  of  J  is 
substantially  larger  than  the  magnitude  of  I  (typically  I  takes  values  up  to  about  5/2),  J 
and  F  will  be  close  to  parallel,  and  the  hyperfine  depolarization  will  be  small.  This 
condition  is  usually  satisfied  for  DFWM  spectra  of  molecular  species  in  which  the 
magnitude  of  J  for  the  levels  of  interest  is  much  larger  than  the  magnitude  of  I.  However 
for  low  values  of  J  (magnitudes  of  J  close  to  or  smaller  than  those  of  I),  the  hyperfine 
depolarization  can  be  significant. 

If  J  is  effectively  decoupled  from  I  by  collisions,  the  hjrperfine  depolarization  is 
small  even  for  low  J  values.  We  used  time-independent  diagrammatic  perturbation  theory 
to  derive  the  results  of  this  chapter,  which  implies  that  the  results  are  valid  for  low- 
intensity,  steady-state  conditions  in  which  population  transfer  is  assumed  to  be  small. 
Inherent  in  this  approximation  is  the  fact  that  the  population  relaxation  and  collisional 
dephasing  times  must  be  much  shorter  than  the  temporal  duration  of  the  laser  pulse. 

These  collisional  times  are  typically  much  shorter  than  the  precessional  period  of  J  about 
F  (usually  tens  of  nanoseconds  or  longer),  which  implies  that  hyperfine  depolarization  will 
not  be  significant  for  collisionally  dominated  systems.  A  quick  experimental  check  is  to 
determine  if  the  hyperfine  structure  is  spectrally  resolved  in  either  the  BB-  or  PC-DFWM 
(sub-Doppler)  phase  matching  geometries.  If  it  is  not  resolved,  hyperfine  interactions  will 
not  significantly  affect  the  results  presented  here.  This  conclusion  follows  because 
unresolved  hyperfine  structure  suggests  that  the  hyperfine  splitting  is  small  compared  with 
the  collisional  broadening  of  the  hyperfine  components,  which  in  turn  implies  that  the 
precessional  period  of  J  about  F  is  slow  compared  with  the  time  it  takes  to  collisionally 
decouple  J  from  I. 

m.  INTERPRETATION  OF  THE  DFWM  SIGNAL  INTENSITY  IN  TERMS  OF 
MULTIPOLE  MOMENTS 

A.  Derivation  and  General  Expressions 

In  DFWM  two  molecular  levels  vdth  total  angular  momentum  quantum  numbers 
Jg  and  Jg,  respectively,  are  coupled  by  the  input  fields.  Describing  the  resulting  coherence 
between  these  levels  in  terms  of  the  total  angular  momentum  distribution  is  convenient. 
The  most  complete  description  of  the  angular  momentum  distribution  is  in  terms  of  state 
multipoles  that  represent  the  populations  of  the  levels  and  the  coherences  existing  between 
them.^^>^^  The  state  multipoles  are  spherical  tensors  of  rank  K  and  component  Q 
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{-K  <Q<K).  The  Q  =  0  components  describe  the  projection  of  J  onto  the  space-fixed 
Z-axis,  and  the  2  0  components  describe  the  projection  of  J  onto  the  space-fixed  XY- 
plane.  Here  the  monopole  term  (K=0)  is  proportional  to  the  population.  All  odd  rank 
multipoles  (dipole,  octopole,  etc.)  describe  the  orientation  of  the  angular  momentum,  and 
all  even  rank  multipoles  (quadrupole,  hexadecapole,  etc.)  are  related  to  the  alignment  of 
the  angular  momentum.  Below  we  show  that  the  highest  allowed  value  of  K  for 
unsaturated  DFWM  is  2.  Therefore,  we  take  the  K=0, 1,  and  2  terms  to  describe  the 
(scalar)  population,  (dipolar)  orientation,  and  (quadupolar)  alignment,  respectively. 

We  outline  a  method  for  calculating  the  third-order  nonlinear  polarization  for  fiilly 
resonant  DFWM  of  molecular  gases  in  collisionally  dominated  environments.  The 
expressions  are  derived  using  diagrammatic  perturbation  theory^^j'^lj^^j^^  and  are 
evaluated  using  a  spherical  tensor  formalism^>^^»^l  The  techniques  used  are  standard; 
consequently,  only  the  key  elements  of  the  derivation  are  presented. 

The  molecular  levels  are  coupled  by  three  input  electric  fields.  The  intensity  of 
each  electric  field  Ej(r,t)  is  given  by 


(16) 


where  the  angular  brackets  denote  a  cycle  average  and  j  is  the  field  label,  is  the 
permittivity  of  free  space,  c  is  the  speed  of  light,  and  Sj  is  the  scalar  amplitude  of  the 
electric  field.  The  result  of  the  interaction  of  the  three  electric  fields  with  the  molecular 
sample  is  the  generation  of  the  third-order  electric  polarization  We  note  that 

many  conventions  are  used  for  expressing  P^^\r,t).  Therefore  for  clarity,  we  follow  the 
convention  of  Butcher  and  Cotter^^  throughout  and  write  SI  units  as 


where  is  the  vector  amplitude,  and  Pofwm  fh®  scalar  amplitude,  and  84  is  the 

normalized  (84  *  64  =  1)  polarization  vector.  P^^fu  is  defined  in  terms  of  the  DFWM 
third-order  nonlinear  susceptibility  as 

4 


(18) 


where  the  symbol  :  refers  to  tensor  contraction.  The  DFWM  signal  intensity  is 
proportional  to  the  cycle  average  of  the  absolute  square  of  the  third-order  nonlinear 
polarization.  In  analogy  with  Eq.  (16),  the  DFWM  signal  intensity  is  proportional  to  the 
absolute  value  squared  of  the  scalar  amplitude  i.e.. 


# 


^  (19) 

Therefore  our  concern  is  to  calculate  pofww  where 

^^DFWM  ~  ^FWM  '^4  ~  ^0  ^  Xz)fWM(~®  4>®  I 3>  ~®  2  )  (^9) 

4 

and  Xd™m(~®4>®i>®  3»~®2)  is  the  scalar  form  of  the  third-order  susceptibility.  See 
reference  62,  page  27. 

The  significant  contributions  to  the  DFWM  third-order  nonlinear  susceptibility  are 
represented  using  double-sided  Feynman  diagrams  in  Fig.  3.  Note  that  sixteen  diagrams 
are  necessary  to  account  for  all  of  the  time  orderings  that  correspond  to  the  energy-level 
dia;^ram  of  Fig.  2.  These  diagrams  represent  the  following  expression^^  for 
Xdfiwi/(~®4>®i>®3>~®2)  in  the  perturbative  (weak-field)  limit: 

XdTOw(“‘®4>®1>®3j~®2)“  ^3  [^(~®4>®1>®3>~®2)'*'^(“®4>®3>®1>~®2)]5  (^l) 

where 

C(-©4,(O,,(l)^,-02) 

=  ^{  X  ){e  |e, -nig  )(^  -5,) 

allM 

X^  X^  |ei-ii|e)(52-^3)  |.(22) 

allM 

and 

^  =(pS-P«)|  [®<,-®  +  *,  v-iT,J[(A^-*2)-v-Jr^J[G)<,-®  +  A:4-v-/r^J  I  ,(23a) 

P«  ){  [w ,  -  ©  +  ^2  •  V  +  iT,^,  fik,  -  ATj  )  •  V  -  zT^  J[w ,  -  ©  +  Ac4  •  V  -  zT^^  ]  |  ,  (23b) 

^2  =  (p^^i-~piv)|  [®o-®+^-v-/r^^][(A:y-^2)-v-/r^,][©o-©+^4-v-/r^^]  |  ,(23c) 
and 


60 


t 


m 


# 


»,  =  (P?i’-p2’){  [®.-®+*i  V  +  ;Tj(*,-*,).v-,r„][o,-<i)+*,.v-,T,,]  I  .  (23d) 

In  Eq.  (22),  A^is  the  total  number  of  absorbers,  the  ket  |m)  represents  the  total  molecular 
wave  function  for  the  quantum  level  characterized  by  total  angular  momentum  J^,  n  is  the 
electric  dipole  moment  operator,  and  A,  and  represent  energy-resonant 

denominators.  In  Eqs.  (23a)-(23d),  is  the  initial  density  matrix  element  for  the 
magnetic  sublevel  of  the  level  n  and  refers  to  the  initial  probability  of  the  system  being 

in  that  sublevel,  h(o^  is  the  energy  difference  between  the  levels  e  and  g,  kj  is  the 
propagation  vector  of  the  field  label  j,  v  is  the  velocity  vector  of  the  absorbing  molecule, 
and  r„^  represents  the  total  dephasing  rate  of  the  coherence  between  the  magnetic 
sublevels  and  A/^. 

We  seek  to  evaluate  Eqs.  (21)-(23)  using  a  spherical  tensor  formalism  and  to 
express  the  DFWM  third-order  nonlinear  polarization  in  terms  of  its  multipole 
components.  This  method  provides  a  well-developed  and  efficient  way  of  using  the 
inherent  symmetry  of  the  system  and  enables  dynamic  and  geometric  factors  to  be 
separated  fi'om  each  other.  The  evaluation  is  similar  to  that  of  coherent  LIF  by  Grreene 
and  Zare  (GZ)^  with  the  exception  that  four  photons  are  involved.  This  greater 
complexity,  however,  does  not  significantly  increase  the  difficulty  of  solving  the  problem, 
which  illustrates  the  power  of  the  spherical  tensor  method. 

The  first  step  in  the  evaluation  is  to  recognize  that  the  summation  over  projections 
onto  the  several  degenerate  sublevels  of  the  level  n  constitutes  a  single  projection  operator 
indicated  by 

P=  2i|«)(w|  =1,  (24) 

which  is  scalar  (tensor  of  rank  0).  Substituting  Eq.  (24)  into  the  first  term  of  Eq.  (22) 
gives 

M, 

where 

O  =  (e:  •  uJ(e^.  •  u ^(e;  •  ^J(E,  •  u,.),  (26) 

and  the  energy-resonant  denominators  have  been  omitted  for  simplicity.  The  transition 
dipole  moment  operators  have  been  subscripted  for  bookkeeping.  The  tensor  product  in 
Eq.  (26)  has  four  dipole  (rank  1)  terms  that  can  be  organized  in  various  ways  to  produce  a 
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tensor  of  rank  4.  In  DFWM,  however,  the  process  initiates  and  terminates  in  the  same 
level,  denoted  here  as  where  is  the  projection  of  the  total  angular 

momentum  J„  (quantum  number  on  the  space-fixed  Z-axis,  and  a„  represents  all  other 
quantum  numbers,  i.e.,  electronic  and  vibrational.  The  only  operator  O  that  can  connect 
the  ket  to  the  bra  |  is  a  scalar,  c,  i.e., 

(a.4M,|0|a,J,M.)  =  (a,y.M,|c|a,J,M,)  =  c  (27) 

because  any  higher  rank  tensor  will  cause  the  matrix  element  to  vanish.  A  tensor  of  rank 
zero  can  be  formed  by  contracting  tensors  of  higher  and  equal  rank.  For  example  two 
tensors  of  rank  1  (vectors)  make  a  tensor  of  rank  0  by  the  vector  dot  product.  Similarly, 
two  vectors  of  rank  two  (quadrupoles)  can  also  be  contracted.  Recall  that  Eq.  (26)  has 
four  dipole  (rank  1)  terms  that  can  be  organized  in  various  ways,  for  instance  as  a  tensor 
product  of  three  dipoles  to  form  tensors  up  to  rank  3  and  a  single  dipole  of  rank  1 . 
However  this  tensor  product  cannot  be  contracted  to  form  a  scalar.  Therefore  we 
organize  Eq.  (26)  in  terms  of  tensor  products  of  two  dipoles  to  form  tensors  up  to  rank  2 
which  can  be  contracted  to  give  a  scalar.  The  resulting  expressions  are  readily  interpreted 
using  the  grating  picture  of  DFWM.  We  note  here,  however,  that  higher-order  tensors 
will  contribute  to  the  DFWM  signal  as  the  laser  intensity  approaches  saturation.  In  a 
perturbative  treatment,  an  increase  in  the  laser  intensity  is  accounted  for  by  including 
higher  order  terms  in  the  perturbative  expansion,  i.e.,  X^>  ••••  X”.  These 

additional  terms  contain  matrix  element  products  with  6,  8,  10, ...  «  + 1  factors  and  can 
form  tensors  up  to  rank  («  + 1)  /  2 . 

The  spherical  tensor  convention  adopted  here  is  that  of  Zare.^  Below  ®  refers  to 
the  tensor  product  and  x  refers  to  scalar  multiplication.  We  begin  by  expressing  the 
vector  dot  product  involved  in  the  electric  dipole  interaction  in  spherical  tensor  notation: 

(e.^)  =  -V3[s<'>®4<'>f  (28) 

Q 

where  the  general  expression  for  the  scalar  product  is 

0 + 1)“2  .  (29) 

Q 

In  Eqs.  (28)  and  (29)  each  vector  is  written  in  terms  of  its  spherical  tensor  components. 
See  the  Appendix.  The  product  of  two  dipole  terms  is 
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(e*  •  n,)(e,.  •  H J  =  3  ®  ®  [e/^  0  j 

=  3((ll)i^(ll)^:0|(ll)0(ll)00)2][[8;^'^®8/'^f^  [^i/'^0n/'^p|“\(3O) 

K  ° 

where  the  operators  have  been  regrouped  using  the  recoupling  transformation  of  four 
angular  momenta.  The  summation  over  K  in  Eq.  (30)  is  from  K=0  to  K=2  because  the 
tensor  product  of  two  dipoles  can  form  a  scalar  (the  dot  product),  a  vector  (the  cross 
product),  or  a  tensor  of  rank  2.  The  factor  in  angular  brackets  is  a  recoupling  coefficient 
and  involves  a  9-j  symbol,  three  of  whose  arguments  are  zero.  This  coefficient  has  been 
evaluated  (GZ  Eq.  (37))  and  has  the  general  form 

((/x')^''(^')^"o|(^)o(r.s:')oo)=(2^+i)'2(2.s:'+i)'2(2^"+i)2 .  (31) 

The  operator  0  defined  in  Eq.  (26)  can  be  represented  using  Eqs.  (30)  and  (31)  as 

0=  ^(2.^ +  1)2  (2^' +1)2 
K,K' 


X 


and  with  one  more  recoupling  transformation  all  of  the  polarization  factors  can  be  grouped 
together.  Performing  this  transformation  gives 

K,K\K” 


X 


(1) 


^A4 


('>  /o>  <') 

P2 


(33) 


which  can  be  simplified  using  Eq.  (31)  and  expanded  in  terms  of  its  tensor  components  by 
Eq.  (29)  to  yield 


# 
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Equation  (34)  is  the  same  as  Eq.  (39)  of  GZ  for  LIF  with  the  important  distinction  that 
four  different  electric  fields  are  involved.  Recalling  that  only  a  scalar  operator,  i.e.,  a 
tensor  of  rank/T'-O,  can  connect  the  ket  to  the  bra  (a„J„M„  |,  we  have  K=K' 

because  only  two  tensors  of  equal  rank  can  be  contracted  to  form  a  scalar.  Therefore  the 
summation  over  K,  K',  and  ^^''of  Eq.  (34)  is  replaced  with  a  single  summation  over  K. 
Substituting  into  Eq.  (25)  gives  a  key  result,  namely. 


Z(“ 


.J.K 


>•4  [I- 


"  J  Jo 


P4  0P, 


where  K  can  have  values  of  0, 1  and  2  (the  range  of  K  and  K'),  and  the  shorthand  notation 
for  the  molecular  wave  functions  |/i)  has  been  replaced  with  Equation  (35) 

displays  the  disentanglement  of  the  geometric  aspects  (polarization  tensors)  of  the 
problem  from  the  dynamics  (transition  moment  tensors)  of  the  electric  dipole  interaction. 

The  first  factor  of  Eq.  (35)  represents  the  scalar  contraction  of  two  polarization 
tensors  to  form  a  scalar  and  is  essentially  in  its  most  fundamental  form.  The  explicit 
evaluation  of  this  factor  is  straightforward  and  is  shown  in  the  Appendix.  The  second 
factor  of  Eq.  (35),  on  the  other  hand,  can  be  reduced  to  a  more  fundamental  form. 
Applying  the  Wigner-Eckart  theorem  (Zare  Eq.  (5.14))  yields  after  some  simplification 


=  (2J^-nY‘{a, 


J  J  K-^w- 


The  last  factor  of  Eq.  (36)  is  the  reduced  matrix  element  of  a  compound  tensor  operator 
and  can  be  evaluated  using  the  following  relation  (Zare  Eq.  (5.74)) 


where  the  term  in  curly  brackets  is  a  6-j  symbol.  Applying  Eq.  (37)  to  (36)  gives  the 
following  result  after  some  simplification 


aJ  a'J'  aV"  ^ 

X  (  a^y^||p<'>||  a'J'  )(  a'J'  ||p;>||  aJ  )(  aJ  ||p('^||  a"J"  }(  a"J"  ||pf^||  a^J^).(38) 

For  fully  resonant  DFWM  we  have  assumed  that  the  excitation  bandwidth  is  sufficiently 
narrow  compared  to  the  density  of  states  (including  Doppler  broadening)  of  the  absorbing 
species  so  that  the  interaction  is  exclusively  between  the  degenerate  magnetic  sublevels  of 
the  two  levels  involved  in  the  one-photon  resonant  transition.  Therefore  the  electric 
dipole  operators  in  Eq.  (38)  can  only  connect  <->  J,,  i.e.,  J"  J^,  J  ^  J^,  and 
J'  ->  J^,  and  Eq.  (38)  simplifies  to 


J,  K 


J" 


(39) 


where  the  labels  on  the  electric  dipole  moment  operator  have  been  dropped.  The  square 
of  the  reduced  matrix  element  (a^yj  molecular  line  strength 

S{aigJg,a^J^)  of  the  ^  o  e  transition,  i.e., 


5(a^7,;a,7,)  =  |(a,7,||  =|(a,7j|  P^'^||a,7,>|^ 

Equation  (40)  is  often  expressed  (Zare  pp.  283-290  and  312-315)  as 


(40) 

(41) 


where  Sjj  is  the  strength  of  the  <->  vibronic  band  (in  many  cases  simply  the 
product  of  the  Franck-Condon  factor(s)  and  the  square  of  the  electronic  transition 
moment),  and  is  the  rotational  line  strength  (H6nl-London  factor). 
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Equations  (29),  (35),  and  (39)  are  used  to  express  the  key  result 


O 


M 
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K  ^  Q  =  -K 


.2  K 


The  first  factor  of  Eq.  (42)  represents  the  dependence  of  the  DFWM  signal  intensity  on 
the  dynamics  of  the  electric  dipole  interaction,  i.e.,  the  greater  the  line  strength  of  the 
molecular  transition,  the  larger  the  DFWM  signal  intensity.  The  second  factor  of  Eq.  (42) 
represents  the  dependence  of  the  DFWM  signal  intensity  on  the  level  degeneracies  and  the 
polarization  vectors  of  the  electric  fields. 

The  final  step  in  the  derivation  of  the  macroscopic  polarization  is  to  average  over 
the  initial  molecular  distribution.  For  an  isotropic  gas  all  of  the  magnetic  sublevels  of  a 
given  level  are  equally  populated  (p,^  =  and  no  phase  relation  exists  between  levels. 
Therefore  the  average  consists  of  integrating  over  the  velocity  distribution  of  the 
absorbing  molecules,  which  is  defined  as  follows 

N  =  J N(v)c/\  =  nJ  fiy)d\  .  (43) 

In  most  experiments  /(v)  is  the  normalized  Maxwell-Boltzmann  velocity  distribution. 
Using  Eqs.  (20)-(23)  and  the  identity 

I’T’O)  0 ^  [^(1)  0  j(i)  (44) 


the  scalar  amplitude  of  the  macroscopic  polarization  becomes 
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K  '■  *''  0=-K 


,2  K 


e=-/: 


i2  /: 


.*(.)^_  (.)F>  [ -o)^  (.)r'^> 

J-G 


a:  ^  Q=-K 

where  the  line  shape  factors  have  the  form 
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I  [(*^-*j)-v-»r,(^)]|  [a),-a)+*^-v-»T^]  [co, -<d +*,•  v+iT,]  |[a),-Q) v-iT^] 


(46) 


In  Eq.  (45),  the  K=0,  1  and  2  terms  describe  the  (scalar)  population,  (dipolar)  orientation, 
and  (quadupolar)  alignment,  respectively,  and  in  Eq.  (46)  the  phenomenological  relaxation 
rates  r,,(Ar)  of  the  spherical  tensor  representation  have  replaced  the  r„,„  rates.  The  r„(Ar) 
are  defined  as  follows:  r„(0)  is  the  relaxation  rate  of  the  global  population  of  the 
level,  and  r„(Ar  =  1)  and  r„(Ar  =  2)  are  the  relaxation  rates  of  the  molecular  orientation 
and  alignment  of  the  level,  respectively.  Furthermore  we  considered  only  collisional 

relaxation  of  the  molecular  system  (no  spontaneous  emission),  and  we  made  the 
reasonable  assumption  that  for  a  dipolar  transition  =  F^^  =  i  ®  >  that  only 

one  optical  relaxation  rate  need  be  considered.^^  As  written,  the  specific  collisional 
relaxation  rates  are  dependent  on  temperature  but  independent  of  velocity. 

Equation  (45)  together  with  Eq.  (46)  represents  the  general  solution  to  the  DFWM 
third-order  polarization  for  molecular  systems  (no  hyperfine  structure)  excited  by  electric 
dipole  radiation  in  the  weak-field  limit  and  can  be  evaluated  explicitly  for  the  phase 
matching  geometry  and  polarization  configuration  of  interest.  Physically,  the  third-order 
DFWM  polarization  of  Eq.  (45)  is  interpreted  as  resulting  from  the  contributions  of  four 
terms:  the  first  term  represents  the  contribution  of  the  diffraction  of  wave  3  from  a 
ground-state  grating  formed  by  fields  1  and  2,  the  second  term  represents  the  contribution 
of  the  diffraction  of  wave  3  from  an  excited-state  grating  formed  by  fields  1  and  2,  the 
third  term  represents  the  contribution  of  the  diflfraction  of  wave  1  from  a  ground-state 
grating  formed  by  fields  3  and  2,  and  the  fourth  term  represents  the  contribution  of  the 
diffraction  of  wave  1  from  an  excited-state  grating  formed  by  fields  3  and  2.  In  addition 
the  multipole  components  can  be  interpreted  as  ground-  and  excited-state  population, 
orientation,  alignment  gratings  for  Ar=0, 1,  and  2,  respectively.  This  treatment 
demonstrates  that  in  the  most  general  case  the  DFWM  signal  may  be  regarded  as  arising 
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from  the  contributions  of  twelve  different  gratings!  The  gratings  can  be  distinguished  by 
spacing  or  Ak^^),  by  the  level  in  which  the  grating  is  formed  (ground  or  excited), 
and  by  the  multipole  nature  of  the  grating  (population,  orientation,  or  alignment). 

To  complete  our  derivation,  Eq.  (45)  is  substituted  into  Eq.  (19)  to  yield  the 
expression^'^  for  the  DFWM  signal  intensity: 

K 


I DFim  ^ 


K 

+  I^^{&,K)G{Jg,J,,K)F(B^,B^,B^,B2,K) 

K 

K 


(47) 


where  the  following  substitutions  have  been  made: 

(2J„+1) 


(48) 
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{^eJe\ 

|p('> 
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(49) 


G(J,r-K)  =  (2J+\) 


\J  J  kV 

ii  1  A ' 


(50) 


and 
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F(e4,8„e^.,82;/i:)= 

i2=-F 


In  Eq.  (50)  the  G(J,J';K)  factors  are  simply  6-j  symbols  and  depend  only  on  K,  Jg,  and 
Jg.  These  values  are  listed  in  Table  3.4  as  a  function  of  transition  type.  Conversely,  the 
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of  polarization  configuration  and  K  value. 
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■^(^456, -,6^,82;^)  factors  ofEq.  (51)  depend  on  the  field  polarizations  but  not  on  and 
Jg.  The  polarization  tensor  products  ofEq.  (51)  are  evaluated  in  the  Appendix,  and  the 
results  for  specific  polarization  cases  treated  in  Sec.  II  are  given  in  Table  3.5. 

B.  Equal  Relaxation  of  the  Multipole  Moments 

In  Sec.  n  we  assumed  that  the  relaxation  is  isotropic  and  that  all  of  the  state 
multipoles  relax  at  the  same  rate.  For  many  experiments,  population  relaxation  rates  such 
as  quenching  are  much  larger  than  collisional  rates  that  only  perturb  the  magnetic  sublevel 
distribution.  Therefore,  a  single  relaxation  rate  might  be  expected  to  characterize  the 
level.  Some  experimental  evidence  supports  this  "single  relaxation"  assumption, 

Under  the  "single  relaxation"  assumption  the  line  shape  factors  ofEq.  (46)  become 
K  independent  and  take  the  form  ofEq.  (7)  of  Sec.  HA,  i.e.,  ,  K)  ->  (co ) . 

Removing  the  line  shape  factors  from  the  summations  ofEq.  (47)  gives  the  familiar 
expression  ofEq.  (2)  with  the  total  geometric  factor  defined  in 

terms  of  an  expansion  in  K  as 

M®) 

+ 1^2  (®  )X  .  63 ,  S2  ;i:) 

K 

^,J;,K)F{z^,z^,z,,z^-,K) 

K 

+  Z|2(co)^G(J^,  J^;Ar)F(84,83,e,,82;Ar)  > 

K  J 


A^2(®)  ^G(7g,y,;Z:)F(84,8„83,82;A:) 


K 

where  the  dependence  of  the  line  shape  factors  on  collisions  and  phase  matching  is  the 
same  as  in  Sec.  II.  Equation  (52)  reduces  to  Eq.  (4)  of  Sec.  HA  because 
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This  identity  is  proven  in  Appendix  A  of  reference  30.  Equation  (53)  results  from  the  fact 
the  first  four-photon  matrix  element  product  in  Eq.  (22)  differs  from  the  second  product 
only  in  the  ordering  of  the  matrix  elements.  Because  each  matrix  element  is  simply  a 
complex  number,  it  is  immediately  obvious  that 

X  }(e  je, -njg  ) 

allM 

=  )(^  }(^  (54) 

allM 

Note  however  that  the  ordering  of  the  angular  momentum  coupling  is  different  for  the  two 
four-photon  matrix  element  products  of  Eq.  (54).  From  Eqs.  (4),  (52),  and  (53),  the 
geometric  factors  are  defined  as 

K 

^'^G{J^,J^-,K)F(z^,Zj,Zi,z^\K).  (55) 

K 

Equation  (55)  was  used  to  generate  the  results  found  in  Tables  I-III. 

The  geometric  factor  Gp^'^^{z^,z^,z.^,Z2’,Jg,Jf,<a)  ofEq.  (52)  is  interpreted  as 
being  proportional  to  the  total  contribution  of  the  2*^  multipole  moment  of  the  total  angular 
momentum  distribution  to  the  DFWM  signal  amplitude  along  84 .  We  have  plotted  the 
G^^^^(e4, 81,83, 82;  J^,J,;co)  for  experimental  conditions  where  I,"2((d)  =  Z^2(®)  and  the 
upper  and  lower  levels  relax  at  the  same  rate,  i.e.,  for  the  cases  where  the  total  geometric 
factor  is  decoupled  from  the  collisional  effects  and  phase  matching  geometry  of  the 
experiment.  These  results  are  shown  in  Figs.  5-8.  Note  that  inherent  in  Eq.  (52)  is  the 
assumption  that  the  relaxation  is  isotropic  and  all  of  the  state  multipoles  relax  at  the  same 
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HG.  3.6.  Evaluated  Gj^'^>(Y,X,Y,X;/^,/,;(0)  =  G;^^>(Y,Y,X,X;7^,7,;(0) 

=— GJ^^^(R,R,R,R;7^,7,;co)  factors  as  a  function  of  Jg  for  (a)  P-branch,  (b)  Q-branch, 

and  (c)  R-branch  transitions.  Permuting  X  Y  or  R  L  does  not  change  the  plots.  In 
the  figure  -A  -O  - ,  and  correspond  to  X=  0, 1,  and  2,  respectively. 


rate.  The  conditions  that  Figs.  5-8  describe  are  very  specific;  however,  the  discussion 
presented  below  can  be  extended  to  most  experiments. 

Figures  5  and  6  show  the  )  values  as  a  function  ofjg 

when  the  fields  have  identical  polarizations,  either  linear  (X  or  Y)  or  circular  (R  or  L). 
According  to  the  grating  picture,  interference  of  the  grating-forming  fields  produces  a 
spatial  modulation  of  laser  intensity  with  spacing  given  by  Eq.  (1 1).  These  intensity 
gratings  are  usually  referred  to  as  population  gratings,  but  as  can  be  seen  in  Figs.  5  and  6, 
the  DFWM  signal  amplitude  has  significant  contributions  from  higher-order  moments.  For 
instance,  the  population  and  alignment  contributions  to  the  DFWM  signal  amplitude  for  a 
Q-branch  transition  excited  with  linear  polarized  light  (see  Fig.  5b)  are  nearly  equal!  This 
and  other  effects  can  be  understood  by  employing  known  concepts  from  linear 
spectroscopy  to  develop  a  qualitative  picture  of  DFWM. 

The  first  step  is  to  recognize  that,  although  the  intensity  is  spatially  modulated,  the 
field  polarization  in  these  regions  is  preserved.  Thus  we  can  think  of  the  initial  isotropic 
distribution  of  absorbers  interacting  with  polarized  light  in  spatially  distinct  regions.  For 
linearly  polarized  light,  we  would  expect  to  create  a  spatial  modulation  of  the  population 
and  alignment  of  the  absorbers,  i.e.,  population  and  alignment  gratings.  As  shown  in  Fig. 
5,  only  population  and  alignment  gratings  contribute.  Likewise  for  circularly  polarized 
light,  we  would  expect  to  create  a  spatial  modulation  of  the  population,  orientation,  and 
alignment  of  the  absorbers,  i.e.,  population,  orientation,  and  alignment  gratings.  Figure  6 
shows  that  this  expectation  is  indeed  met. 

Another  interesting  case  involves  linear  polarized  light  in  which  both  pump  fields 
have  the  same  polarization  but  are  orthogonally  polarized  with  respect  to  the  probe  field 
(see  Fig.  7).  The  interference  produced  in  this  case  is  spatially  uniform  in  intensity  but 
spatially  modulated  in  polarization.  This  polarization  grating  is  typically  referred  to  as  a 
coherence  grating  because  spatial  population  modulation  does  not  exist.  Under  these 
circumstances  only  orientation  and  alignment  contributions  are  expected,  as  confirmed  in 
Fig.  7. 

In  Fig.  7  orientation  gratings  are  found  to  dominate  for  A/  =  ±1  transitions  and 
alignment  gratings  dominate  for  A/  =  0  transitions.  Furthermore,  inspection  of  Figs.  6-8 
shows  that  this  conclusion  is  generally  true.  This  behavior  can  be  understood,  again,  if  we 
consider  the  interaction  of  an  isotropic  distribution  of  molecules  interacting  with  polarized 
light.  After  the  absorption  of  linearly  polarized  light  the  system  will  be  aligned,  and  after 
the  absorption  of  circularly  polarized  light  the  system  will  be  oriented  and  aligned.  Q- 
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FIG.  3.7.  Evaluated  Gp^’^\Y,X,X,Y;Jg,J/,(0)  factors  as  a  function  of  Jg  for  (a)  P- 

branch,  (b)  Q-branch,  and  (c)  R-branch  transitions.  Permuting  X  «->  Y  does  not  change 
the  plots.  In  the  figure  -O  -  and  correspond  to  X  =  1  and  2,  respectively. 


branch  transitions  are  more  easily  aligned  than  oriented,  and  the  opposite  is  true  for  P-  and 
R-branch  transitions. 

Finally,  the  most  significant  aspect  of  Figs.  5-8  is  that  the  relative  multipole 
moment  contributions  are  very  dependent  on  the  value  of  y  for  low  lvalues,  but  rapidly 
approach  a  high-J  limit,  i.e.,  for  J  >  4  the  total  geometric  factor  and  its  relative 
contributions  are  essentially  independent  of  J.  This  behavior  means  that  the  J  dependence 
of  the  DFWM  signal  intensity  for  a  given  branch  (P,  Q,  or  R)  taken  in  a  single  polarization 
configuration  will  primarily  reflect  the  J  dependence  of  the  one-photon  absorption 
coefficient  Bgg  as  well  as  any  J-dependent  relaxation  effects.  Molecules  are  often  in  the 
high-J  limit,  especially  at  flame  and  plasma  temperatures  at  which  lvalues  ranging  from  5 
to  50  are  typically  populated. 

C.  Unequal  Relaxation  of  the  Multipole  Moments 

We  have  assumed  thus  far  that  the  multipole  moments  relax  at  equal  rates,  which  is 
a  good  approximation  in  flame  experiments  where  efficient  energy  transfer  collision 
partners,  such  as  water,  are  present.  We  have  recently  observed^^  experimental 
conditions  for  which  these  rates  are  not  the  same,  particularly  for  low  values  of  J,  in 
collisional  environments  in  which  light  species  such  as  helium  are  the  dominant  collision 
partners.  Such  species  have  slow  inelastic  collisional  rates  (quenching)  compared  with 
elastic  collisional  rates  (depolarizing).  The  presence  of  unequal  relaxation  can  be 
determined  by  polarization-ratio  measurements. 

In  the  event  that  the  multipole  moments  of  the  angular  momentum  distribution  do 
not  relax  equally,  the  expressions  for  the  DFWM  signal  intensity  presented  in  Sec.  II  do 
not  apply.  If  the  relaxation  is  isotropic,  however,  the  multipole  moments  relax 
independently,  and  the  expressions  presented  in  Sec.  IIIA  are  valid.  The  effects  of 
unequal  relaxation  has  been  discussed  by  Ducloy  and  Bloch^®  and  is  the  topic  of  a  future 
publication.^^  Finally  if  the  relaxation  is  anisotropic,  for  example,  the  gas  is  in  a  static 
electric  or  magnetic  field,  then  the  situation  is  considerably  more  complicated.^^  Note, 
however,  that  the  interpretation  of  DFWM  signal  intensities  at  this  level  of  detail  is  not 
necessary  for  most  applications. 


IV.  COMPARISON  WITH  PREVIOUS  WORK 


In  what  follows  we  attempt  to  place  our  work  in  the  context  of  previous 
perturbative  theoretical  treatments  of  DFWM.  Lam  and  Abrams  (LA)29  were  among  the 
first  to  describe  the  polarization  effects  for  a  degenerate  two-level  system.  They  used  a 
density  matrix  approach  to  express  the  DFWM  response  as  a  summation  of  matrix  element 
products  that  represented  different  coherences  corresponding  to  contributions  fi'om 
normal  population,  cross-population,  and  Zeeman  coherence  mechanisms.  Although  this 
model  provides  an  intuitive  picture  for  the  origin  of  the  DFWM  polarization  effects,  its 
application  depends  on  the  coordinate  system  chosen.  This  coordinate-frame  dependence 
also  makes  the  inclusion  of  collisional  relaxation  difficult,  and  LA  acknowledged  that  an 
accurate  description  of  the  DFWM  collisional  effects  in  systems  with  degenerate  sublevels 
necessitates  the  introduction  of  the  irreducible  (spherical  tensor)  representation  of  the 
density  matrix.^^>52 

Such  an  approach  was  presented  by  Ducloy  and  Bloch  (DB),30  who  discussed  the 
DFWM  response  in  terms  of  its  irreducible  tensor  components.  In  their  tensorial  density 
matrix  treatment,  DB  considered  the  degeneracy  of  the  resonant  levels,  unequal  relaxation 
of  the  magnetic  sublevels,  spontaneous  emission  (optical  pumping),  thermal  motion,  and 
the  effect  of  pump-probe  angular  separation  (noncollinear  geometries).  In  later  papers 
Berman,  Steel,  Khitrova,  and  Liu  (BSKL)^^  and  Alekseev^  ^  derived  similar  expressions 
that  also  included  the  effects  of  hyperfine  structure.  Together  the  work  of  DB,  BSKL, 
and  Alekseev  illustrates  the  fundamental  aspects  of  DFWM  in  the  perturbative  regime  and 
represents  general  solutions  of  the  problem.  Unfortunately  these  studies  have  not  been 
available  to  most  experimentalists  because  the  final  results,  expressed  in  density  matrix 
notation,  are  quite  complicated.  Recently  Kupiszewska  and  Whitaker  (KW)^^  addressed 
this  problem  by  specializing  the  treatment  of  BSKL  to  stationary  diatomic  molecules 
described  by  Hund's  case  (a)  coupling  with  no  hyperfine  structure.  The  simplified 
expressions  presented  by  KW  provide  some  insight  into  the  polarization  effects  in  DFWM 
for  molecular  systems;  however,  the  specialization  of  the  treatment  necessarily  limits  the 
applicability  of  the  results. 

An  alternative  approach  to  interpret  DFWM  signal  intensities  is  to  use 
diagrammatic  perturbation  theory  (double-sided  Fe5mman  diagrams)  to  evaluate  the 
DFWM  polarization  and  collisional  effects.^O, 41,58  jhe  diagrammatic  method  is 
equivalent  to  the  density  matrix  formalism  but  has  the  advantage  that  the  third-order 
nonlinear  polarization  can  be  evaluated  directly  without  the  need  to  calculate  lower-order 
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processes.  Using  diagrammatic  methods,  Attal-Tretout,  Monot,  and  Muller-Dethlefs 
(TMD)33  and  Aben,  Ubachs,  van  der  Zwan,  and  Hogervorst  (AUZH)^^  have  derived 
expressions  for  resonant  CARS  spectroscopy.  This  work  focused  on  obtaining  four-wave 
mixing  line  strength  factors  for  diatomic  molecules  (omitting  nuclear  spin)  to  be  used 
directly  in  the  analysis  of  experimental  signal  intensities.  Both  TMD  and  AUZH  note  that 
the  CARS  expressions  could  be  extended  to  describe  DFWM,  i.e.,  fully  resonant  CARS. 
We  noted  in  Sec.  II  that  the  line  strength  factors  themselves  are  dependent  on  the 
relaxation  caused  by  the  environment  and  not  just  angular  momentum  considerations. 
Consequently,  this  suggestion  is  correct  only  with  qualification. 

The  general  CARS  expressions  involve  four  distinct  frequencies  o,,  Oj,  (O3,  and 
0)4,  and  four  independent  states  |a),  |6},  |c),  |t/)  (labeled  a,  n,  h,  n'  in  TMD  and  a,  b,  c,  d 
in  AUZH).  In  the  fully  resonant  case,  ffl,=c)2=o)3=(i)4=(D,  and  the  four  independent 
states  converge  to  four  distinct  magnetic  sublevels,  i.e.,  |a)  |g),  |i)  -» \e),  |c)  |^'), 

and  =  |e').  This  mapping  leads  to  the  inclusion  of  eight  diagrams  (eight  more  if 
N^^O)  and  four  distinct  resonant  four-photon  matrix  element  products,  namely, 

|el-li|^’)(^'|e3-n|g)(^'|8;-n|e  ){e  |e, ),  (56a) 

allM 

2]  (^'|e3'p|^')(g’|e*-p|e  )(e  je,  -pjg  )[g  (56b) 

allM 

X  hi  h3-pk},  (56c) 

allM 

and 

2]  h’hi -^k’X^  1^3 X^  (^^d) 

allM 

Equations  (56a)  and  (56b)  have  the  same  value  when  summed  over  all  M,  but  the  ordering 
of  the  angular  momentum  coupling  is  different,  and  thus  they  represent  different 
contributions  from  the  zeroth,  first,  and  second  rank  multipole  moments.  The  same  is  true 
for  Eqs.  (56c)  and  (56d).  The  proper  accounting  of  the  multipole  moment  contributions, 
to  our  knowledge,  has  not  been  done  using  diagrammatic  perturbation  theory  but  has  been 
done  using  density  matrix  approaches. In  addition  Eqs.  (56a)  and  (56c)  appear 
quite  similar,  but  have  distinctly  different  values  even  when  summed  over  all  M,  and  the 


same  is  true  for  Eqs.  (56b)  and  (56d).  Therefore  care  must  be  taken  in  extending  CARS 
expressions  to  DFWM. 

Bervas,  Le  Boiteux,  Labrunie,  and  Attal-Tretout  (BBLT)^^  and  Freidman-Hill, 
Rahn,  and  Farrow  (HRF)^^  recently  reduced  the  CARS  expressions  to  apply  to  DFWM. 
The  expressions  of  BBLT  contradict  our  results  for  BB-  and  PC-DFWM  but  agree 
fortuitously  with  our  results  for  FB-DFWM.  We  offer  the  following  explanation  for  these 
results.  In  Fig.  4  of  BBLT,  diagrams  (a)-(d)  describe  the  primary  contributions  to  BB- 
and,  PC-DFWM  in  the  Doppler-broadened  limit,  and  diagrams  (a)-(h)  describe  the 
primarily  contributions  to  FB-DFWM.  Because  BBLT  do  not  distinguish  between  g  and 
g'  or  e  and  e’,  all  of  the  diagrams  appear  to  be  represented  by  the  same  matrix  element 
product,  namely  Eq.  (8)  of  BBLT.  Diagrams  (a)  and  (c)  of  BBLT  Fig.  4  correspond  to 
our  Eq.  (56a),  diagrams  (b)  and  (d)  of  BBLT  Fig.  4  correspond  to  our  Eq.  (56c),  diagrams 
(e)  and  (g)  of  BBLT  Fig.  4  correspond  to  our  Eq.  (56b),  and  diagrams  (Q  and  (h)  of 
BBLT  Fig.  4  correspond  to  our  Eq.  (56d).  In  terms  of  the  notation  developed  in  this 
chapter,  the  result  of  BBLT  predicts  the  following  proportionality  for  the  BB-  and  PC- 
DFWM  signal  intensity: 


I  QC 


2 


(57a) 


which  is  independent  of  the  difference  in  the  relaxation  rates  of  the  g  and  e  levels  for  all 
polarization  configurations,  whereas  our  treatment  predicts  that 


which  is  dependent  on  the  relative  relaxation  rates  of  the  g  and  e  levels  whenever  the  two 
geometric  factors  differ  (see  Sec.  II).  Furthermore,  our  result,  Eq.  (57b),  agrees  with  Eq. 
(45)  of  DB  derived  using  a  density  matrix  approach  and  experimental  DFWM  polarization 
ratios  obtained  in  our  laboratory  for  the  CH  radical.  The  expressions  for  FB-DFWM 
presented  by  BBLT  agree  with  our  results  because  Zfj  )  =  ^2(® )  A*2  (^)-  ^2(®  )• 

In  our  notation,  BBLT  predict  the  following  proportionality  for  the  FB-DFWM  signal 
intensity: 


2 

G^(e4„E„63,e2;7g,J.){Zfj((D)  +  Zi'2((D))  +  Gp(e4,E3,e„ej;Jg,J,))(Z|2(©)  +  Z;j(co))|  ,  (58a) 


which  is  equivalent  to  our  result 


G^(E„E„E3,E2;J^,J.)(If2((D)  +  ^;2(o)))  +  G^(E„e3,E„E2;J^,J.))(Z^3(0))  +  I|2((o))|^  (58b) 


81 


f 


# 


« 


when  Zf2((D)  =  Z|2((o)  and  HRF  considered  three-level  or  "crossover" 

resonances  that  occur  when  two  transitions  that  share  a  common  level  become  resonant 
because  of  a  simultaneous  Doppler  shift.  HRP  show  that  "crossover"  resonances  can 
influence  DFWM  signal  intensities  even  if  they  are  not  spectroscopically  resolved  from 
their  parent  transitions.  Such  considerations  could  be  important  in  spectrally  congested 
regions  with  overlapping  branches. 

In  summary,  the  work  to  date  can  be  classified  in  two  categories:  first,  density 
matrix  approaches  that  emphasize  the  complex  relation  between  polarization,  collisional 
relaxation,  and  phase  matching  at  the  expense  of  simplicity;  and  second,  diagrammatic 
approaches  that  focus  on  expressing  the  polarization  effects  (without  specifying  multipole 
contributions)  analyticly  while  treating  relaxation  effects  superficially.  The  emphasis  of 
this  chapter  has  been  on  the  key  aspects  of  interpreting  molecular  DFWM  spectra  in 
collisionally  dominated  environments  because  these  types  of  experiments  are  those  for 
which  DFWM  has  found  the  largest  application.  In  these  experiments,  hyperfine  structure, 
spontaneous  emission,  and  in  most  circumstances  unequal  relaxation  of  the  multipole 
moments  of  the  angular  momentum  distribution  can  be  neglected  so  that  relatively  simple 
expressions  are  obtained  without  unnecessarily  trivializing  the  treatment.  Thus  we  have 
attempted  to  combine  the  best  aspects  of  the  two  approaches  by  presenting  analytic 
expressions  with  the  relevant  collision,  velocity,  and  polarization  considerations  for  the 
most-utilized  experimental  configurations.  These  expressions  can  be  used  directly  in  the 
evaluation  of  DFWM  spectra  for  these  specific  cases.  Finally  the  polarization  tensor 
products  evaluated  in  the  Appendix  can  be  used  to  extend  the  general  expressions  of  Sec. 
niA  to  other  experimental  geometries  of  interest. 

The  above  results  were  derived  in  the  weak-field  (no  saturation)  limit.  An 
appealing  aspect  of  DFWM  is  that  in  the  saturation  regime,  the  DFWM  signal  intensity 
becomes  relatively  insensitive  to  the  specific  value  of  the  relaxation  rates.  ^  In  this 
regime,  the  problem  of  extracting  relative  population  distributions  essentially  reduces  to 
knowing  the  absorption  coefficients.  This  aspect  is  discussed  in  detail  in  Chapter  4. 

V.  CONCLUSIONS 

In  the  present  chapter  we  derived  expressions  via  time-independent  diagrammatic 
perturbation  theory  that  account  for  the  DFWM  polarization,  collisional,  and  velocity 
effects  in  the  weak-field  limit  (no  saturation).  In  our  treatment,  we  assumed  that  the 
DFWM  process  couples  levels  of  sharp  (definite)  angular  momentum  J.  Three  input  fields 
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of  arbitrary  polarization  interact  with  an  isotropic  sample  to  produce  a  fourth  field.  The 
general  result  (Sec.  Ill)  was  specialized  to  apply  to  circularly  and  linearly  polarized  fields 
that  interact  in  nearly  collinear  phase  matching  geometries  in  collisional  environments 
where  the  multipole  moments  of  the  total  angular  momentum  distribution  relax 
independently  (isotropic  relaxation)  and  at  the  same  rate.  These  specialized  expressions 
(Sec.  II)  generally  apply  to  DFWM  experiments  performed  in  coUisionally  dominated 
environments. 

In  Sec.  II  we  showed  that  the  DFWM  signal  intensity  for  coUisionally  dominated 
systems  is  proportional  to  the  square  of  the  concentration  difference  of  the  levels  involved 
in  the  one-photon  resonant  transition,  the  fourth  power  of  the  one-photon  transition 
strength,  the  square  of  a  total  line  shape  function,  and  the  square  of  the  total  geometric 
factor.  The  total  geometric  factor  was  shown  to  depend  not  only  on  the  polarization  of 
the  input  fields  but  also  on  the  specifics  of  the  experiment,  i.e.,  collisional  relaxation, 
velocity,  and  phase  matching.  This  interdependence  complicated  the  interpretation  of  the 
DFWM  polarization  response  but  also  enabled  qualitative  information  about  the  collisional 
relaxation  caused  by  the  environment  to  be  obtained  from  polarization  measurements. 

Finally  in  Sec.  HI,  we  showed  that  in  the  most  general  case  the  DFWM  signal  may 
be  regarded  as  arising  from  the  contributions  of  twelve  different  gratings.  These  gratings 
are  distinguished  by  spacing  (AA,2  or  the  level  in  which  the  grating  is  formed 

(ground  or  excited),  and  by  the  multipole  nature  of  the  grating  (population,  orientation,  or 
alignment).  Therefore  simple  population  vs.  coherence  grating  pictures  of  DFWM  are  not 
adequate  descriptions.  For  instance,  DFWM  experiments  in  which  all  fields  are  the  same 
polarization  (either  linear  or  circular)  have  contributions  (in  addition  to  population)  from 
higher-order  moments  of  the  total  angular  momentum  distribution,  namely,  alignment  for 
linear,  and  orientation  and  alignment  for  circular.  We  also  showed  that  the  relative 
multipole  moment  contributions  depend  strongly  on  J  for  low  7  values  but  rapidly 
approach  a  high-7  limit.  In  this  limit  the  7  dependence  of  the  DFWM  signal  intensity  for  a 
given  branch  (P,  Q,  or  R)  taken  with  a  single  polarization  configuration  primarily  reflects 
the  7  dependence  of  the  one-photon  absorption  coefficient  Bgg  as  well  as  any  7-dependent 
relaxation  effects. 

We  have  discussed  the  key  aspects  of  DFWM  in  reference  to  molecular  species  in 
coUisionally  dominated  environments.  We  presented  explicit  expressions  for  the  DFWM 
signal  intensity  that  can  be  used  with  molecular  absorption  and  relaxation  data  to  obtain 


the  relative  population  distributions.  Absorption  data  in  the  form  of  cross  sections,  line 
strengths,  emission  coeflScients,  etc.,  are  avmlable  for  most  molecules  of  interest. 
However  obtaining  an  accurate  and  complete  set  of  relaxation  rates  for  a  given  molecule 
in  a  specific  environment  is  a  formidable  task. 
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CHAPTER  4 


Effects  of  Field  Polarization  in  the  Reduction  of  DFWM  Spectra  to 
Relative  Populations  in  the  Presence  of  Weak  and  Strong  Fields 


There  are  no  absolutes  in  combustion  just  a  lot  of  pretty  sures. 


Personal  Opinion 


I.  INTRODUCTION 


Abrams  and  Lind  (AL)1  have  presented  a  model  of  degenerate  four- wave  mixing 
(DFWM)  which  has  shown  good  agreement  with  experiments.  2-4  The  AL  model 
considers  a  nondegenerate  two-level  homogeneously  broadened  system  in  the  presence  of 
two  counterpropagating  pump  fields  of  arbitrary  intensity  and  weak  probe  and  signal 
fields.  The  probe  and  signal  intensities  are  weak  in  the  sense  that  they  do  not  significantly 
affect  the  level  populations.  In  this  model  all  fields  are  monochromatic,  are  of  the  same 
polarization  state,  and  are  configured  in  the  phase-conjugate  geometry  (PC-DFWM) 
discussed  in  Chapter  3.5 

One  of  the  most  significant  predictions  of  the  AL  model  is  that  the  DFWM  signal 
becomes  relatively  insensitive  to  collisions  when  the  intensity  of  the  pump  fields  is 
increased  so  that  the  population  difference  of  the  two-level  system  oscillates  with  a  rate 
(Rabi  frequency)  greater  than  the  relaxation  and  dephasing  rates.  When  the  pump  field 
intensity  approaches  this  value  (7^,)  the  DFWM  signal  is  said  to  be  saturated.  The 
collisional  independence  of  saturated  DFWM  has  been  verified  experimentally  for  many 
species  including  NO,^>^  NH,!^  andCH.^2  Note  however  that  the  AL  model 

should  be  viewed  only  as  a  qualitative  description  of  these  experiments  because  in  most 
cases  the  lasers  employed  are  not  monochromatic,  the  molecules  are  moving,  and  the 
probe  field  intensity  is  a  significant  fi'action  of  the  pump  field  intensity.  In  addition, 
molecules  possess  polarization-dependent  spectroscopic  properties  such  as  those 
discussed  in  Chapter  3  that  are  not  described  by  the  AL  model. 

Meacher,  Charlton,  Ewart,  Cooper,  and  Alber  (MCECA)15-15  h^ve  developed  a 
theory  for  DFWM  with  broad-bandwidth  lasers  that  has  shown  qualitative  agreement  with 
experiments.  mCECA  model  is  similar  to  the  AL  model  in  that  intense  pump  fields 

and  a  weak  probe  field  interact  with  stationary  nondegenerate  two-level  absorbers  in  a  PC- 
DFWM  phase  matching  geometry;  however,  the  pump  bandwidths  are  assumed  to  be 
much  larger  than  the  homogeneous  line  width  (including  power  broadening)  of  the  optical 
transition.  In  this  model  the  probe  bandwidth  can  be  any  size.  The  results  show  that 
increasing  the  pump  bandwidth  leads  to  an  increased  and  a  slightly  different  DFWM 
power  dependence  compared  to  the  AL  model. 

Lucht,  Farrow,  and  Rakestraw  (LFR)1^  have  addressed  the  effects  of  molecular 
motion  and  strong  (saturating)  probe  field  intensity.  LFR  numerically  evaluated  the  time- 
and  space-dependent  density  matrix  equations  directly  for  a  nondegenerate  two-level 
system  with  molecular  motion  in  which  any  or  all  of  the  incident  fields  can  be  saturating. 
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LFR  calculated  DFWM  intensities  and  line  shapes  for  conditions  of  interest  for  flame  and 
plasma  diagnostics,  and  compared  the  results  of  their  calculations  to  experimental  results 
ofNO  in  100  Torr  of  helium  excited  by  a  narrow-bandwidth  (0.004  cm‘1)  laser.  Some  of 
the  findings  of  LFR  relevant  to  our  work  are  as  follows:  first  the  presence  of  molecular 
motion  slightly  increases  7^°,;  second,  when  the  probe  field  intensity  is  one  quarter  of  the 
pump  field  intensity  the  DFWM  saturation  behavior  is  very  similar  to  the  case  where  the 
probe  field  is  held  constant  at  a  value  much  less  than  7^,;  and  third,  the  DFWM  signal  is 
nearly  independent  of  collisions  when  the  pump  field  intensity  is  approximately  21°^,. 

LFR  discussed  their  work  in  the  context  of  previous  saturation  treatments  ^  ^"24  so 
further  comparisons  will  not  be  made  here.  The  work  of  AL,  MCECA,  and  LFR  greatly 
enhance  the  understanding  of  DFWM  saturation;  however,  because  these  models  assume  a 
nondegenerate  two-level  system,  the  effect  of  polarization  in  saturated  DFWM  remains 
unresolved. 

Our  aim  is  to  determine  the  effect  of  polarization  in  appl3dng  nondegenerate  two- 
level  models  to  experiments  involving  molecules;  in  particular  we  wish  to  determine  how 
polarization  affects  extracting  relative  internal-state  distributions  from  DFWM  signals. 

Our  approach  is  to  measure  DFWM  signals  as  a  function  field  polarization  and  intensity 
and  to  interpret  these  data  in  the  framework  of  the  AL  and  MCECA  models  using  the 
ideas  developed  in  Chapter  3.  In  Sec.  II  we  discuss  the  pertinent  experimental  details.  A 
brief  discussion  of  the  AL  and  MCECA  models  is  given  in  Sec.  Ill  with  an  emphasis  on 
the  key  elements  required  to  interpret  DFWM  signals.  In  Sec.  IV  we  present  DFWM  data 
of  the  CH  radical  in  an  atmospheric-pressure  flame  as  a  function  of  field  polarization  and 
intensity,  interpret  the  results,  and  present  methods  for  extracting  relative  population 
distributions  from  saturated  DFWM  signals.  In  Sec.  V  we  present  conclusions  of  our 
findings. 

n.  EXPERIMENTAL  DETAILS 
A.  Experimental  Apparatus 

Figure  1  shows  the  experimental  configuration  in  which  three  fields  of  the  same 
fi’equency  overlap  in  the  reaction  zone  of  an  atmospheric-pressure  oxyacetylene  flame  to 
generate  a  fourth  field  (the  DFWM  signal).  The  two  coaxial  and  counterpropagating 
pump  fields,  denoted  E,  (forward  pump)  and  (backward  pump),  are  crossed  at  a  small 
angle,  0  =  2°,  by  a  probe  field,  E^.  This  configuration  is  commonly  referred  to  as  the 
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FIG.  4.1.  Experimental  configuration  for  the  measurement  of  DFWM  signals.  Spectra  are 
recorded  with  unfocused  beams  (diameter  ~  0.5  mm)  having  energies  of  16  pJ  for  the 
pump  beams  (1  and  3)  and  4  pJ  for  the  probe  beam  (2).  The  flame  consists  of  a  standard 
welding  torch  fitted  with  a  0.94  mm  nozzle.  The  abbreviations  used  in  the  figure  are 
defined  as  follows;  M:  mirror,  P:  polarizer,  L:  lens,  PH;  pinhole,  PEM;  photoelastic 
modulator,  BS;  beamsplitter,  WP;  half-wave  plate,  DFR;  double  Fresnel  rhomb,  EM; 
energy  meter,  PMT;  photomultiplier  tube. 


DFWM 
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phase-conjugate  geometry  (PC-DFWM).  The  PC-DFWM  signal  field,  E^,  is  coherently 
generated  and  propagates  6  m  to  a  photomultiplier  (Hamamatsu  R2393P)  where  it  is 
detected. 

Excitation  is  provided  by  a  KrF  excimer-pumped  dye  laser  system  (Lambda  Physik 
EMG  102  MSC  and  FL  2002)  operated  over  the  wavelength  range  of 425-450  ran  with  a 
bandwidth  of  Au^=0.18  ±  0.02  cm"l  (fwhm)  and  a  pulse  width  of  At^=16.3  ±  1  ns 
(fwhm).  The  output  beam  of  the  dye  laser  is  spatially  filtered  to  produce  a  Gaussian  beam 
profile.  The  fwhm  diameter  of  the  beams  is  measured  using  the  procedure  outlined  in 
reference  25  and  is  typically  -0.5  mm.  The  beam  diameter  can  be  determined  to 
±  0.02  mm  .  The  laser  energy  is  computer  (PC/XT)  controlled  using  a  photoelastic 
modulator  (Hinds  International  PEM-CF4)  between  two  crossed  polarizers  that  allows  the 
laser  energy  to  be  adjusted  between  1  and  800  pJ  with  no  beam  deviation.  The  energy  of 
the  probe  beam  is  one  quarter  of  the  forward  pump  beam,  and  the  energy  of  the  backward 
puinp  beam  is  adjusted  to  the  value  of  the  forward  pump  beam  using  a  double  Fresnel 
rhomb  and  a  polarization  beamsplitter  which  reflects  S-polarized  and  transmits  P-polarized 
light.  We  measure  the  beam  energy  by  placing  an  energy  meter  (Molectron  J4-09)  after 
the  polarization  beamsplitter  which  allows  the  energy  to  be  measured  directly  and 
accurately  to  a  fraction  of  a  microjoule. 

The  DFWM  signal  is  amplified  (Lecroy  WIOOBTB)  before  being  processed  by 
boxcar  a  averager  (SRS  SR250).  The  laser  energy  is  processed  directly  by  a  boxcar 
averager  (SRS  SR250).  Up  to  four  channels  are  read  at  each  wavelength,  digitized,  and 
averaged  over  30  laser  shots.  The  data  is  then  stored  on  a  computer  (PC/XT)  for  further 
analysis.  If  not  otherwise  stated,  the  reported  errors  include  statistical  and  experimental 
uncertainties  and  represent  la  deviation. 

B.  Flame 

The  atmospheric-pressure  oxyacetylene  flame  consists  of  a  standard  welding  torch 
fitted  with  a  0.94  mm  diameter  nozzle  operated  in  the  open  air.  The  flows  of  acetylene 
and  oxygen  are  maintained  with  mass  flow  controllers  (MKS)  at  380  cc/min  and  360 
cc/min,  respectively.  This  results  in  an  oxygen  to  acetylene  volumetric  flow  ratio, 
O2/C2H2,  of  0.947  and  an  equivalence  ratio,  C>,  of  2.64  (fuel  rich).  The  temperature  is 
2774  ±  140  K  at  0.4  mm  above  the  burner  nozzle  and  the  CH  concentration  is  estimated 
to  be  30  ppm.  12  At  these  concentrations  (-8  x  lO^  molecules/cm^  per  quantum  state)  no 
CH  absorption  (<  1%)  is  observed. 
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C.  CH  Spectral  Features  and  Collisional  Rates 


Figure  2  shows  a  DFWM  spectrum  of  the  CH  radical  taken  at  0.4  mm  above  the 
burner  nozzle.  The  integrated  signal  intensities  of  the  well-resolved  (0,0)  band  transitions 
of  the  A^A-X^n  system^b  are  used  in  our  analysis.  The  A-X  system  is  highly 
diagonal, and  both  the  X  and  A  states  rapidly  approach  Hund's  case  (b)  coupling  as 
rotation  increases.  The  radiative  lifetimes  of  the  v -0  level  of  the  A^A  state  are  nearly 
independent  of  rotational  level.^^*^^  In  recent  calculations,  however,  Luque  and 
Crosley29  show  that  there  is  about  a  4%  variation  in  the  radiative  lifetimes  for  N'  quantum 
numbers  ranging  from  N-2  to  N-20.  Therefore  we  use  the  calculations  of  Luque  and 
Crosley  and  the  rotational  line  strengths  of  Kovacs^2,33  calculate  the  Einstein 

absorption  coefficients  used  in  our  analysis.^^  Experimental  term  energies^^  are  used  for 
J"  quantum  number  values  up  to  12.5.  Term  energies  for  higher  J"  values  are  calculated 
using  the  spectroscopic  constants  of  reference  35. 

High-temperature  collisional  data  for  the  A^A(v-0)  state  of  the  CH  radical  are 
shown  in  Table  4. 1 .36-41  ju  Table  4. 1,  refers  to  rotational  energy  transfer  and  Q  refers 
to  quenching.  Electronic  and  vibrational  energy  transfer  are  not  included  because  these 
rates  are  negligible  compared  to  R  and  When  only  cross  sections  are  reported,  those 
of  CO  and  H2  are  used  to  calculate  the  rate  coefficient.  These  two  molecules  are 
considered  because  Matsui,  Yuuki,  and  Sahara^^  have  determined  by  mass-spectrometric 
analysis  that  the  reaction  zone  of  oxyacetylene  torch  burners  consists  primarily  of  CO  and 
H2.  The  total  collisional  transfer  rate  coefficient  7Jj  (s*^)  from  energy  level  /  to  energy 
level  j  of  species  p  due  to  collisions  with  several  species  labeled  q  is  given  by 

T.=yNva^‘’  (1) 


where  is  the  number  density  (molecules/cm^)  of  species  q,  is  the  cross-section 
(cm2),  and  is  the  average  relative  speed  (cm/s)  of  species  p  and  q.  The  average 

relative  speed  is  given  by 


%  = 


^2 


Ttp 


M  J 


(2) 


where  kg  is  the  Boltzmann  constant,  7  is  the  temperature,  and  is  the  reduced  mass  for 
species  p  and  q.  To  our  knowledge  the  highest  temperature  data  available  are  at  2000  K. 
The  temperatures  of  interest  in  our  investigations  are  approximately  2774  K.  Therefore 
the  rates  are  scaled  to  this  temperature  by  making  a  IjT  density  correction  and  a 


FIG.  4.2.  DFWM  spectrum  of  selected  Q-  and  R-branch  transitions  of  the  CH  A^A- 
X^n  (0,0)  and  (1,1)  bands  recorded  using  the  YYXX  polarization  configuration.  In 
the  figure  only  the  (0,0)  band  transitions  are  labeled.  The  fine-structure  components 
are  indicated  on  the  branch  designations  by  subscript  1  =  ^^5/2-^113/2  and  subscript 
2  =  ^^3/2“^ni/2.  The  upper  state  and  lower  state  A  doublets  are  labeled  by  e  and  f; 
when  they  are  the  same  for  both  levels,  ee  is  abbreviated  by  e  and  fif  is  abreviated  by 
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429.0  429.5  430.0  430.5  431.0  431.5 

Wavelength  (nm) 


Table  4.1.  CH  energy  transfer  rates. 


Ref 

CH  Source 

Temperature 

R/Q 

0^ 

37 

Flame 

2040  K 

6,  14 

2.4 +  0.1,  4.1 +0.4 

36 

Flame 

2000  K 

11 

4b 

5.1  +0.4 

38 

Pyrolysis 

1300  K 

average 

— 

5.6  ±0.4 

39 

Flame 

1600  K 

6 

3.6  ±  0.5 

3.6  ±2.4 

39 

Flame 

1600  K 

3,  12 

2.7  ±0.5,  3.8 +  0.5 

40 

Flame 

1800  K 

average 

— 

4.7  ±0.3 

41 

Flame 

1700  K 

average 

— 

1.8  ±  0.4c 

^Estimated  for  CH  at  2774  K  in  units  of  10*  s*l.  ^Error  not  reported.  ^Estimated  error. 


velocity  correction.  Such  a  method  is  accurate  if  the  cross-sections  do  not  change 
significantly  with  velocity.  Strictly  speaking,  the  T^J  (s'^)  have  units  of  a  first-order 
chemical  kinetic  rate  coefficient;  however,  they  are  more  commonly  referred  to  as  rates  in 
nonlinear  optics.  Therefore  in  what  follows  we  will  refer  to  the  collisional  transfer  rate 
coefficients  T]j  as  rates. 

D.  Saturation  Measurements 

In  Chapter  1  saturation  measurements  were  performed  by  obtaining  spectral  scans 
at  specific  laser  intensities.  The  line-center  signal  intensities  (peak  heights)  were  easily 
distinguished  fi^om  the  baseline  noise.  Because  of  the  time  required  to  take  a  spectral  scan 
at  each  laser  intensity  the  data  accumulation  time  was  quite  long  (~2  hours).  We  observe 
for  our  present  experimental  conditions  that  the  major  source  of  noise  is  scattered  light. 
We  are  not  able  to  observe  a  coherent  nonresonant  background  signal  in  these 
experiments.  Furthermore,  the  scattered  light  results  primarily  fi-om  beam  2  (probe) 
because  it  scatters  light  directly  along  the  detection  axis.  A  lesser  contribution  from  beam 
3  (backward  pump)  is  also  observed  due  to  distortion  and  refi’action  of  the  beam  after 
passing  through  the  flame.  Scattered  light  from  beam  1  (forward  pump)  is  found  to  be  on 
the  order  of  the  electronic  noise  of  the  detection  system,  which  is  much  less  than  the  noise 
from  the  scattered  light  of  beams  2  and  3. 

In  the  present  chapter  we  vary  the  laser  energy  with  the  laser  frequency  held  at  line 
center.  Specifically,  the  KrF  excimer  laser  trigger  is  phase  locked  with  the  50  kHz 
oscillations  of  the  photoelastic  modulation  crystal.  Varying  the  laser  trigger  delay  varies 
the  retardation  of  the  light  exiting  the  crystal,  and  hence,  the  amount  of  light  that  passes 
through  the  second  polarizer.  The  laser  energy  is  varied  under  computer  coritrol,  and  the 
total  light  intensity  reaching  the  detector  and  laser  energy  are  recorded  each  laser  shot. 
Then  in  a  separate  scan  the  forward  pump  beam  is  blocked,  and  a  scattered  light  signal  is 
obtained  as  a  function  of  laser  energy.  Both  traces  are  sorted  as  a  function  of  laser  energy 
(corrects  for  shot-to-shot  fluctuations)  and  binned  at  intervals  of  0.5  |xF.  The  upper  and 
lower  solid  symbol  traces  in  Fig.  3  are  examples  of  the  total  signal  and  scattered  light 
signal,  respectively.  The  scattered  light  signal  (linear  in  laser  energy)  is  subtracted  from 
the  total  signal  to  obtain  the  DFWM  signal  as  a  function  of  laser  energy  (open  symbols  in 
Fig.  3).  This  method  allows  an  entire  saturation  curve  to  be  obtained  in  approximately 
five  minutes  with  each  data  point  consisting  of  an  average  of  60-90  laser  shots. 


FIG.  4.3.  Line-center  saturation  curve  for  the  YYXX  linear  polarization  configuration  of 
the  Qief(^)  transition  as  a  function  of  saturation  intensity  (bottom),  and  pump  field 
spectral  intensity  (top).  The  upper  and  lower  solid  symbol  traces  in  the  figure  are  the  total 
signal  (-■-)  and  scattered  light  (-A-)  signal,  respectively.  The  scattered  light  signal  (linear 
in  laser  intensity)  is  subtracted  from  the  total  signal  to  obtain  the  DFWM  signal  as  a 
function  of  laser  intensity  (-0-)- 
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In  calculating  the  laser  spectral  intensity  from  the  measured  beam  energies,  we  use 
Au^=0.18  ±  0.02  cm’^,  ATj^=16.3  ±  1  ns,  and  A  =  where  dj.^  is  the  "top  hat" 

diameter  of  a  Gaussian  beam^^  (dj^j  equals  1.2  d^^„  and  typically  dj^f^=  0.5  mm).  For 
example,  typical  beam  energies  employed  in  these  experiments  are  16  |xJ  for  the  forward 
and  backward  pump  beams  and  4  pJ  for  the  probe  beam.  These  values  correspond  to 
1.93  MW/cm^-cm'l  and  0.48  MW/cm^-cm'l,  respectively.  In  what  follows  is  the 
theoretical  line-center  saturation  intensity,  /"  is  the  experimental  line-center  saturation 
intensity,  and  <  >  is  the  average  experimental  line-center  saturation  intensity  for 

different  polarization  configurations  of  single  transition. 

E.  Polarization  Measurements 

The  experiment  is  configured  so  that  any  linear  polarization  configuration  is 
obtained  by  rotating  the  half-wave  plates  shown  in  Fig.  1 .  After  polarization  rotation  each 
pump  beam  is  directed  by  one  additional  mirror.  Differences  in  S  and  P  reflectance  for  all 
the  mirrors  is  approximately  3%.  The  beamsplitter  in  the  probe  axis,  which  is  used  to 
collect  the  DFWM  signal,  reflects  P-polarized  light  at  66.4  +  6%  times  the  efficiency  for 
S-polarized  light.  The  probe  mirror,  however,  reflects  P-polarized  light  at  a  higher 
efficiency  than  S-polarized  light  which  results  in  less  than  5%  difference  in  the  probe  field 
intensity  at  the  flame  as  its  polarization  is  rotated.  The  polarization  purity  of  all  of  the 
electric  fields  is  greater  than  98%. 

Shot-to-shot  polarization  measurements  are  performed  by  placing  the  photoelastic 
modulator  in  place  of  one  of  the  half-wave  plates  and  removing  the  polarizer  from  the 
detection  axis.  The  laser  delay  relative  to  the  photoelastic  modulator  is  adjusted  on  a 
shot-to-shot  basis  so  that  every  other  shot  corresponds  to  either  0°  or  180°  retardation 
(half-wave  plate).  The  DFWM  signal  and  laser  energy  are  averaged  over  30  laser  shots 
for  each  polarization  configuration  and  for  each  spectral  position.  The  shot-to-shot 
measurements  enable  polarization  ratios  to  be  obtained  in  real-time  at  a  single  laser 
intensity.  These  measurements  are  used  to  check  the  accuracy  of  polarization 
measurements  taken  a  few  minutes  apart. 

In  what  follows  Y  represents  the  linear  polarization  state  that  corresponds  to  S- 
polarized  light,  i.e.,  light  polarized  perpendicular  to  the  plane  of  the  laser  table,  and  X 
represents  P-polarized  light.  In  addition  we  use  the  notation  established  in  Chapter  3 
where  a  polarization  configuration  is  given  as  84e,8382  where  Zj  is  the  polarization  vector 


of  the  electric  field  labeled  j.  For  example  the  YYXX  polarization  configuration 
corresponds  to  and  £,  being  Y  polarized  and  and  E^  being  X  polarized. 

F.  Rotational  Temperature  Measurements 

In  our  rotational  temperature  analysis  resolved  Qj  and  Q2  transitions  of  the 
a2A(v -0)-X2n(v"=0)  system  are  used.  For  CH  the  Q  branch  covers  a  small  spectral 
region  and  is  sufficiently  resolved  to  allow  the  analysis  of  individual  DFWM  transitions. 
Therefore  rotational  temperatures  are  obtained  with  minimal  scan  times.  The  XXXX  and 
YXYX  polarization  configurations  have  been  used  in  previous  experiments,  3,6-12  but  we 
chose  the  YYXX  configuration  because  it  offers  many  advantages.  First  the  YYXX 
polarization  configuration  greatly  discriminates  against  thermal  grating  contributions 
because  only  small-spaced  intensity  gratings  are  present.^’^^  We  have  not  observed  any 
direct  evidence  to  suggest  that  thermal  gratings  are  produced  in  this  flame;  however,  this 
polarization  configuration  assures  that  their  interference  is  negligible.  Second,  as  stated  in 
Sec.  lie,  the  primary  source  of  noise  in  these  experiments  is  scattered  light  from  the  probe 
beam  (-90%)  and  the  backward  pump  beam  (-10%)  with  a  negligible  contribution  (on  the 
order  of  the  electronic  noise)  fi'om  the  forward  pump  beam.  In  the  YYXX  polarization 
configuration,  the  signal  beam  is  crossed  polarized  to  both  of  the  "noisy"  beams,  and 
therefore  the  YYXX  configuration  offers  true  zero  background  detection  when  a  linear 
polarizer  is  placed  in  the  detection  axis.  Contrary  to  first  impressions,  the  PC-DFWM 
signal  intensity  using  the  YYXX  polarization  configuration  is  equal  to  the  YXYX 
configuration.  Both  these  configurations  give  approximately  1/9  the  signal  intensity  of  the 
XXXX  configuration  (Table  3.3),  but  the  thermal-grating  and  scattered-light 
discrimination  make  up  for  the  loss  of  signal. 

m.  SATURATION  MODELS 

The  AL  model  considers  a  nondegenerate  two-level  homogeneously  broadened 
system  in  the  presence  of  arbitrary  pump  field  intensities  and  weak  probe  and  signal 
intensities.  In  this  model  all  fields  are  assumed  to  be  monochromatic.  In  the  limit  of  weak 
probe  and  signal  intensities  a  full  solution  of  the  density  matrix  equation  for  the  interaction 
of  the  pump  fields  with  the  two-level  system  is  possible.  The  density  matrix  solution  for 
two-level  nondegenerate  saturable  absorbers  in  the  presence  of  arbitrary  field  intensities  is 


well  known  in  semi-classical  laser  theory.  An  excellent  discussion  of  the  topic  is  given  in 
reference  45. 

For  electric  fields  defined  as 


where  Ej  is  the  vector  amplitude,  is  the  scalar  amplitude,  Sj  is  the  normalized 
(e^  •  E*  =  1)  polarization  unit  vector,  and  Ij  is  the  intensity  of  the  electric  field  labeled  j, 
the  density  matrix  solution  in  the  rotating-wave  approximation  for  the  steady-state 
population  (number  density)  difference  =  Ng-N^\s 
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The  Rabi  frequency  |Q^  |  associated  with  the  field  labeled  j  is  defined  with  respect  to  the 
field  intensity/^  as 

|2 

-Ij,  (5) 

and  the  saturation  Rabi  frequency  =  (For^)'''^  is  defined  with  respect  to  the  line- 
center  saturation  intensity  1°^,  as 
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In  Eqs.  (3)-(6)  g  refers  to  the  lower  level,  e  refers  to  the  upper  level,  AN^  -Ng-  N°  is 

||2 
ItgJ  is 

the  square  of  the  transition  dipole  moment  connecting  levels  g  and  e,  |Q|  is  the  Rabi 
frequency  associated  with  the  pump  fields  (|0|  =  |n,|  =  la|)>  ^  **  the  intensity  of  the  pump 
fields  (/  =  /,=  /j),  k  is  the  magnitude  of  the  propagation  vector  of  the  fields,  A  =  o)  -  co  „ 
is  the  spectral  shift  from  the  resonance  frequency  (o^,  and  F^  and  F^^  are  the  population 
relaxation  and  coherence  dephasing  rates,  respectively.  All  of  the  terms  are  in  SI  units, 
and  the  constants  e^,  fi  and  c  have  their  usual  meanings. 

In  Eq.  (4)  we  see  that  the  result  of  the  interference  of  the  two  counterpropagating 
pump  fields  is  a  sinusoidally  varying  electric  field  intensity  that  spatially  "bums  holes"  into 
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the  population  difference  spaced  one-half  wavelength  apart.^^  The  two-level  system  is 
said  to  be  saturated  when  the  population  difference,  =^Ng-N,,  is  zero  over  the  entire 
interaction  volume.  Inspection  of  Eqs.  (4)  and  (6)  shows  that  saturation  occurs  when  the 
optical  pumping  rate  denoted  by  the  Rabi  frequency  is  much  faster  than  the  population 
relaxation  and  coherence  dephasing  rates. 

Having  obtained  the  system  response  in  the  presence  of  strong  pump  fields, 
perturbation  theory  is  used  to  include  the  weak  probe  and  signal  fields.  The  final  result  is 
an  analytic  model  of  the  DFWM  signal.  AL  solved  the  slowly  varying  envelope 
approximation  (SVEA)  equations  for  the  probe  and  signal  fields.  In  the  limit  of  no 
absorption,  the  SVEA  equations  are 


dz 


=  0, 
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•4  _ 


dz 


=  /3g* 
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where  P  is  the  nonlinear  coupling  coefficient  defined  in  reference  2,  page  222,  Eq.  (16). 
For  the  case  of  equal  intensity  pumps,  the  expression  for  the  DFWM  signal2>47,48  jg 
given  by 

^DPmf  ~  |3|  L 1^=  R^i  I 2 »  (8) 
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r,=r„(i+4//c)''’.  (11) 

In  Eqs  (7)-(l  1)  is  the  AL  model  reflectivity,  a^is  the  line-center  attenuation 
coefficient,  L  is  the  effective  interaction  length  of  the  excitation  fields,  7  =  7,  =  /j  is  the 
intensity  of  the  pump  fields,  and  7^  is  the  intensity  of  the  probe  field.  Equation  (8)  defines 
the  cube  of  a  power-broadened  Lorentzian  with  an  actual  half-width  given  by 

r.(2=^-i)'". 
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The  formulation  of  the  MCECA  model  is  similar  in  concept  to  the  AL  model  with 
the  exception  that  the  excitation  fields  are  not  monochromatic  (coherent).  The  MCECA 
model  considers  a  medium  of  two-level  nondegenerate  absorbers  transversed  by  two 
broad-bandwidth  (incoherent)  copolarized  pump  fields  of  equal  and  arbitrary  intensity  and 
a  weak  probe  field  of  any  bandwidth.  Furthermore  MCECA  assume  that  the  fields  are 
statistically  independent  (chaotic)  and  uncorrelated  (delayed  with  respect  to  each  other  by 
a  distance  much  greater  than  the  coherence  length).^^"^  1  The  bandwidths  of  the  pump 
fields  have  a  Lorentzian  spectral  shape  with  a  fwhm  of  A© ,  =  Ao  3  =  2Z> .  Most 
importantly,  this  model  assumes  that  the  bandwidths  exceed  all  other  rates  that  determine 
the  time  evolution  of  the  system,  i.e.,  b  » |n|,r„,r,^ .  In  this  case  the  bandwidth- 
dependent  saturation  Rabi  frequency  is  given  by  and  the  average  density 

matrix  population  terms  vary  on  a  time  scale  of  l/T^  whereas  the  much  shorter  correlation 
time  of  the  fields  varies  as  Mb.  Under  these  conditions,  all  products  of  field  and 
population  terms  may  be  treated  as  uncorrelated.  This  aspect  is  the  key  element  of  the 
treatment. 

In  contrast  to  the  AL  model  expression  of  Eq.  (4)  which  involves  coherent  fields, 
the  interference  of  the  counterpropagating  pump  fields  considered  by  MCECA  does  not 
produce  a  steady-state  standing  wave  pattern  because  the  fields  are  incoherent  and 
uncorrelated.  Replacing  the  cos^  (kz)  term  of  Eq.  (4b)  with  its  average  value  of  1/2  and 
assuming  that  b  »  ,  we  obtain  the  expression  for  the  for  the  steady-state  population 

difference  AA^  =  Ng-N^  of  the  MCECA  model,  namely. 


M  = 


AA/„ 


l+ml/rj 


(12a) 


with 


2 

^  = - 7-T, 


(12b) 


where  is  defined  by  Eq.  (7)  with  We  see  in  Eq.  (12)  that  stronger  fields 

are  required  to  saturate  the  two-level  system  when  broad-bandwidth  lasers  are  used 
because  the  Rabi  fi’equency  must  now  exceed  (bF^y^^  which  is  much  larger  than  (F^Fij)'^^. 

Analytic  expressions  for  the  DFWM  signal  of  the  MCECA  model  are  obtained 
when  the  bandwidth  of  the  probe  field,  A©  2  =  2/),  is  also  large  in  the  sense  that 
|Q|^  «  pb  and  p  »  F^,F^g .  Specifically  for  p=b  and  for  small  detunings  compared  to  the 
laser  bandwidth  (A  «  6),  we  have  the  simplified  MCECA  expression  for  the  DFWM 
signal,  namely. 
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r^=s(i+v//cf,  (15) 

r„=i(i+2//cf,  (16) 

and 

r„  =  *(i+(2+v)//Cf.  (17) 

In  Eqs.  (13)-(17)  is  the  ratio  of  the  population  relaxation  rate  to  the  coherence 
dephasing  rate  =  r„/r,2 ),  and  the  remaining  terms  have  been  defined  previously.  For 
7  ^  3/j",,  Eq.  (13)  approximates  a  power-broadened  Lorentzian-cubed  profile  like  the  AL 
model  but  with  a  spectral  wdth  determined  by  the  laser  bandwidth.  For  higher  laser 
intensities,  the  line  shape  power  broadens  more  rapidly  than  a  Lorentzian-cubed  profile 
and  develops  a  dip  at  line  center.  The  intensity  range  over  which  Eq.  (13)  applies  depends 
on  the  laser  bandwidth  and  the  population  relaxation  rate  caused  by  the  environment. 
Because  b  »  |n|and|n„l  =ir„  we  find  that  the  condition 


(18) 


needs  to  be  satisfied  for  Eq.  (13)  to  be  valid. 

The  collisional  relaxation  and  dephasing  rates  for  our  experimental  conditions  must 
be  determined  in  order  to  apply  the  models  discussed  in  this  section.  We  base  the 
discussion  on  the  available  relaxation  data  for  the  CH  radical  (Sec.  IIC),  but  the  discussion 
is  also  applicable  to  other  flame  and  plasma  species.  For  molecules  in  collisionally 
dominated  environments,  the  rotational  levels  involved  in  the  one-photon  transition  are 
effectively  coupled  by  rotational  energy  transfer  (R)  collisions  to  other  rotational  levels. 
Furthermore,  the  excited  rotational  levels  are  coupled  to  the  ground  rotational  levels  by 
quenching  (Q)  and  spontaneous  emission  (Agg).  For  this  case  the  relaxation  of  the 
molecular  levels  approximates  a  four-level  system.  ^2  poj.  such  a  system  the  population 
relaxation  and  coherence  dephasing  rates  are  defined^^  as 
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r«=5(r,+r,)+C,  (20) 

where  Fg  and  are  the  population  relaxation  rates  of  levels  g  and  e,  and  is  the  pure 
dephasing  rate  of  the  dipolar  transition. 

In  our  analysis  we  take  R  to  be  the  same  for  the  A^A  (v -0)  and  X^n  (v"=  0)  levels 
because  the  rotational  spacings  are  very  similar  and  define  the  ground  and  excited 
population  relaxation  rates  as  Fg  =  and  F^  =  g  respectively.  With  these 
definitions  we  estimate  an  average  population  relaxation  rate  of  F„=  2  x  10^  using  the 
values  of  Table  4.1.  Assuming  that  the  pure  dephasing  rate  is  small  compared  to  R,  we 
have  Fo  =  F^g  =  0.98  for  R/Q  =  3.5).  Furthermore  because  the  characteristic 
collisional  relaxation  times  are  ~  40  times  shorter  than  the  laser  pulse,  i.e.,  Ax^.  » 

1/Fg  =  1/F^ ,  steady-state  conditions  are  assumed,  and  we  use  the  AL  and  MCECA 
expressions  discussed  in  this  section.  53  in  addition  we  expect  the  MCECA  model  to 
apply  for  laser  intensities  up  to  approximately  twice  the  saturation  intensity  because 
b/T^  =  9  in  Eq.  (18).  In  Eqs.  (8)  and  (13)  AN^  is  replaced  with  the  ground  state 
population  because  N° »  N°  for  our  flame  conditions. 

IV.  RESULTS  AND  DISCUSSION 

A.  Polarization  Ratios  as  a  Probe  of  Collisional  Dynamics 

In  Chapter  3  we  discussed  how  DFWM  polarization  measurements  provide 
information  about  the  collisional  relaxation  caused  by  the  environment.  Line-center 
saturation  curves  like  those  of  Fig.  3  for  the  Rif(8)  and  Qief(8)  transitions  were  obtained 
for  different  linear  polarization  configurations.  The  line-center  saturation  curves  of  the 
Rlf(8)  transition  are  shown  in  Fig.  4.  The  absolute  signal  intensities  dramatically  differ  for 
the  four  polarization  configurations.  To  emphasize  the  differences  in  signal  intensities, 
polarization  ratios  as  a  function  of  laser  intensity  are  made  from  the  saturation  curves. 
Figure  5  shows  polarization  ratios  of  the  Rif(8)  and  Qief(8)  transitions  as  a  function  of 
pump  field  intensity.  The  solid  symbols  in  Fig.  5  are  experimental  data  and  dashed  lines  in 
Fig.  5  are  the  predicted  polarization  ratios  using  Table  3.3. 

Figure  5  a  shows  the  XYYX/XXXX  polarization  ratio  of  the  Rif(8)  transition,  and 
Fig.  5b  shows  the  YYXXATXYX  polarization  ratio  of  the  Qief(8)  transition.  The 
predicted  polarization  ratios  in  Figs  5a  and  5b  are  in  good  agreement  with  the  data  at 
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FIG.  4.4.  Line-center  saturation  curves  for  different  linear  polarization  configurations  of 
the  Rif(8)  transition  as  a  function  of  average  saturation  intensity  <  /“  >=  1.77  MW/cm^- 
cm’i  (bottom),  and  pump  field  spectral  intensity  (top).  In  the  figure  -A-,  and 
correspond  to  the  XXXX,  YXYX,  YYXX  and  XYYX  polarization  configurations, 
respectively.  The  absolute  signal  intensities  dramatically  differ  for  the  four  polarization 
configurations,  but  the  saturation  intensities  (Table  II)  obtained  from  the  AL  model 
(dashed  curves)  and  MCECA  model  (solid  curves)  are  nearly  equal  for  all  polarization 
configurations. 
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FIG.  4.5.  The  solid  points  are  the  experimental  polarization  ratios  and  the  dashed  lines  are 
the  predicted  YYXX/YXYX  polarization  ratios  using  Table  3.3.  (a)  XYYX/XXXX 
polarization  ratio  for  the  Rif(8)  transition  as  a  function  of  average  saturation  intensity 
<  -C/  >=  1'77  MW/cm2-cm-l.  The  average  value  of  the  data  between  0. 1  and  0.4  <  7“  > 

^  sat 

IS  0.20  ±  0.06  which  is  in  good  agreement  with  the  calculated  value  of  0.24.  (b) 
YYXXAi'XYX  polarization  ratio  for  the  Qief(8)  transition  as  a  function  of  average 
saturation  intensity  <  7“  >=  1.10  MW/cm^-cm'l.  The  average  value  of  the  data  between 
0. 1  and  0.4  <  7“  >  is  0.92  ±  0. 14  which  is  close  to  the  expected  value  of  1 .0.  (c) 
YYXX/YXYX  polarization  ratio  for  the  Rj£(8)  transition  as  a  function  of  average 
saturation  intensity  <  7“  >=  1.77  MW/cm^-cm*^  The  average  value  of  the  data  between 
0.1  and  0.4  <  7“  >  is  1.04  +  0.09  which  is  close  to  the  expected  value  of  1.0. 
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weak  pump  field  intensity.  Specifically,  the  average  value  of  the  data  between  0. 1  and  0.4 
<  /“  >  is  0.20  +  0.06  in  Fig  5a  and  0.92  +  0. 14  in  Fig.  5b.  These  values  are  close  to  the 
predicted  polarization  ratios  of  0.24  and  1.0,  respectively.  Recall  that  the  values  of  Table 
3.3  were  derived  assuming  that  all  of  the  multipole  moments  of  the  total  angular 
momentum  distribution,  i.e.,  the  population,  the  orientation,  and  the  alignment,  relax  at  the 
same  rate  and  that  no  other  process  contributes  the  observed  signal. 

We  have  observed  Q-branch  polarization  ratios  different  from  the  values  of  Table 
3.3  for  the  OH  radical  in  the  burnt  gas  region  of  atmospheric-pressure  hydrogen-oxygen 
flames.  We  have  attributed  this  effect  to  the  presence  of  thermal  gratings.  With 
helium  dilution,  however,  thermal  gratings  are  no  longer  present,  but  anomalous 
polarization  ratios  are  still  observed  for  levels  with  low  J  quantum  numbers.  Nearly 
degenerate  four-wave  mixing  (NDFWM)  scans  confirmed  that  in  the  helium-diluted  flames 
the  rates  of  relaxation  of  the  multipole  moments  are  not  equal  .  ^6  Therefore  the  good 
agreement  between  the  data  at  low  intensity  in  Figs.  5a  and  5b  with  the  calculated  values 
implies  that  the  mutipole  moments  are  relaxing  equally  and  that  other  laser-induced 
phenomena  such  as  a  thermal  grating  are  not  present. 

In  PC-DFWM  of  Doppler  broadened  systems,  polarization  ratios  of  P-  and  R- 
branch  transitions  are  dependent  on  the  relative  values  of  the  relaxation  rates  of  ground 
and  excited  levels  when  the  pump  fields  are  cross  polarized,  i.e.,  when  the  YYXX  and 
YXYX  polarization  configurations  are  used.  We  estimate  a  collisional  full  width,  2r^^,  of 
4x10^  s"l  (0.02  cm*l)  and  a  Doppler  width,  Aco  p,  of  46  x  10^  s’^  (0.24  cm*l)  for  CH  at 
2774  K.  Our  experiment  corresponds  to  a  Doppler-broadened  system  because  A©  ^  > 

2r,g .  For  such  a  system  the  YYXX/YXYX  polarization  ratio  of  the  Rif(8)  transition  has 
a  value  of  1.4  when  «  T^,  0.70  when  »  F^,  and  1.0  when  F^  =  F^. 

Figure  5c  illustrates  the  YYXX/YXYX  polarization  ratio  for  the  Rif(8)  transition, 
and  the  dashed  line  of  Fig.  5c  is  the  prediction  of  table  3.3  which  corresponds  to  F^  =  F^. 
The  average  value  of  the  data  between  0.1  and  0.4  <  /“  >  is  1.04  +  0.09  which  suggests 
that  the  relaxation  rates  of  the  ground  and  excited  levels  are  approximately  equal.  This 
conclusion  is  consistent  with  the  energy  transfer  data  available  for  the  CH  radical.  The 
rotational  spacings  of  the  v=0  vibrational  levels  of  the  X  and  A  states  of  CH  are  nearly  the 
same,  and  the  rotational  energy  transfer  rate  of  the  v'=0  level  of  the  A  state  is 
approximately  3.5  times  faster  than  quenching.  Furthermore  if  the  quenching  collisions  do 
not  trap  population  in  higher  vibrational  states  (rapid  energy  transfer  to  v"=0),  we  expect 
that  F^  s  F^  according  to  the  four-level  model  referred  to  in  Sec.  III. 
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An  interesting  aspect  of  Fig.  5  is  that  the  experimental  polarization  ratios  are 
within  30%  of  the  perturbative  calculations  up  to  twice  the  saturation  intensity.  In  the 
weak-field  limit  the  DFWM  signal  is  dependent  on  the  spatial  anisotropy  of  the  total 
angular  momentum  distribution  induced  by  the  three  input  fields.  Therefore  this  close 
agreement  implies  that  the  relative  anisotropy  of  the  total  angular  momentum  represented 
by  the  total  geometric  factors  persists  to  a  large  degree  as  the  laser  intensity  increases. 

This  effect  is  similar  to  that  expected  for  saturated  laser  induced  fluorescence  measured 
with  different  polarization  configurations  (Fig.  5  of  reference  57). 

B.  Saturation  Behavior 

The  line-center  saturation  intensities  resulting  from  nonlinear  least  squares  fits  of 
saturation  data  to  the  power  dependencies  of  the  AL  model,  Eq.  (8)  and  the  MCECA 
model,  Eq.  (13)  for  A  =  0  and  =  0.98  are  given  in  Table  4.2.  In  Table  4.2,  the  total 
error  increases  significantly  as  the  ratio  of  the  DFWM  signal  to  the  scattered  light  intensity 
decreases.  The  AL  model  and  the  MCECA  model  exhibit  very  similar  saturation  behavior 
for  moderately  high  laser  intensities  (/-/",).  This  statement  is  validated  by  the  fact  that 
the  fitted  saturation  curves  (Fig.  4)  and  saturation  intensities  (Table  4.2)  of  the  two 
models  are  very  similar.  Therefore  it  is  not  surprising  that  the  AL  model  has  been  able  to 
qualitatively  reproduce  the  saturation  behavior  for  experiments  employing  broad- 
bandwidth  (multi-mode)  lasers.  12  The  error  associated  with  the  MCECA  saturation 
intensities,  however,  is  approximately  30%  smaller  than  the  error  associated  with  the  AL 
saturation  intensities.  The  reduced  error  results  from  the  fact  that  the  fits  to  MCECA 
model  are  more  sensitive  to  the  value  of  the  saturation  intensity  (higher  order  power 
dependence)  and  are  better  correlated  with  the  data. 

For  a  single  transition,  the  fitted  MCECA  saturation  intensities  of  Table  4.2  are 
approximately  equal  regardless  of  the  polarization  configuration.  For  instance,  the 
saturation  intensities  for  the  four  polarization  configurations  of  the  Rif(8)  transition  are 
within  30%  of  the  average  value  (/,“)  of  1.77  MW/cm^-cm'l.  This  similarity  is  more 
easily  seen  in  Fig.  6  which  shows  the  DFWM  reflectivity  at  line  center  for  the  four 
polarization  configurations  of  the  Rif(8)  transition.  The  MCECA  line-center  reflectivity  is 
at  a  maximum  at  0.6 1°^,  and  in  Fig.  6,  the  reflectivity  for  each  polarization  configuration 
peaks  at  approximately  the  same  pump  field  intensity.  This  effect  is  independent  of  the 
fact  that  the  absolute  magnitude  of  the  reflectivities  differ  for  the  various  polarization 
configurations. 


FIG.  4.6:  Line-center  reflectivity  curves  for  different  linear  polarization  configurations  of 
the  Rif(8)  transition  as  a  function  of  average  saturation  intensity  <  /“  >=  1.77  MW/cm^- 
cm'l  (bottom),  and  pump  field  spectral  intensity  (top).  In  the  figure  -A-,  and 
correspond  to  the  XXXX,  YXYX,  YYXX  and  XYYX  polarization  configurations, 
respectively.  The  solid  curves  are  nonlinear  least  squares  fits  of  the  data  to  the  MCECA 
model,  Eq.  (14).  All  of  the  reflectivity  curves  peak  at  approximately  the  same  pump  field 
intensity. 
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Table  4.2.  Experimentally  determined  line-center  saturation  intensities. 


Transition 

ALC 

(MW/cm^-cm'l) 

MCECA  C 
(MW/cm^-cm*^) 

OleK-l) 

YYXX 

1.88  ±1.19 

1.76  ±0.74 

Rlf(4) 

YYXX 

2.47  ±  0.83 

2.21  ±0.54 

YXYX 

1.06  ±0.38 

1.12  ±0.25 

Qlel<8) 

YYXX 

1.00  ±0.51 

1.08  ±0.32 

Rli(8) 

XXXX 

1.88  ±0.33 

1.79  ±0.21 

Rlf(8) 

YXYX 

2.30  ±0.43 

2.07  ±  0.29 

Rli(8) 

YYXX 

1.68  ±0.41 

1.72  ±0.25 

Rlf(8) 

XYYX 

1.21  ±1.34 

1.51  ±0.81 

We  suggest  the  following  explanation  for  the  weak  sensitivity  of  the  DFWM 
saturation  intensity  to  polarization.  In  the  framework  of  the  MCECA  discussed  in  Sec. 
in,  the  DFWM  signal  saturates  when  the  population  difference  of  Eq.  (12)  induced  by  the 
counterpropagating  pump  fields  approaches  zero  over  the  entire  interaction  volume.  At 
this  point  the  probe  field  cannot  perturb  the  population  (make  gratings)  and  the  DFWM 
reflectivity  approaches  zero.  The  intensity  at  which  the  spatial  population  difference  of 
the  two-level  system  approaches  zero  is  determined  by  Because  for 

the  MCECA  model,  this  intensity  is  related  to  the  collisional  relaxation  caused  by  the 
environment  and  the  Einstein  absorption  coefficient  (||i^^|  =  (3eji^ /%)Bg^).  For  a 

nondegenerate  two-level  system  this  intensity  is  unambiguous. 

The  transition  probabilities  for  a  molecular  (degenerate)  transition,  however,  are 
dependent  on  the  magnetic  sublevels  involved,  the  type  of  transition,  and  polarization  state 
of  the  excitation  field.  In  other  words,  for  a  given  transition  and  field  polarization,  some 
magnetic  sublevels  interact  strongly  with  the  field  and  others  do  not.  Those  sublevels 
interacting  strongly  with  the  field  saturate  more  rapidly  than  those  that  do  not.  Therefore 
to  describe  the  saturation  properties  of  a  molecular  system  each  magnetic  sublevel  must  be 
treated  independently,  i.e.,  there  is  no  single  Rabi  frequency  for  a  degenerate  system. It 
is  common  in  rate  equation  models  of  molecular  processes,  however,  to  use  Einstein  A 
and  B  coefficients^^  which  represent  an  equally  weighted  average  of  the  transition 
probabilities  of  all  the  magnetic  sublevels.^^>^^"^®  Furthermore  if  the  initial  magnetic 
sublevel  distribution  is  isotropic  the  average  optical  pumping  rate,  i.e.,  the  rms  Rabi 
frequency,  is  independent  of  the  field  polarization. 

Therefore  to  a  first  approximation  the  average  intensity  at  which  the  population 
difference  approaches  zero  will  be  related  to  the  Einstein  absorption  coefficient. 
Furthermore  we  expect  this  intensity  to  be  similar  for  the  different  polarization  cases 
provided  that  the  magnetic  sublevels  relax  with  equal  rates,  i.e.,  if  the  population, 
orientation,  and  alignment  relax  equally.24  Sec.  IVA  we  showed  this  to  be  a 
reasonable  assumption  for  CH  in  this  flame.  The  saturation  data  presented  in  Table  4.2 
validates  this  expectation  because  the  saturation  intensities  for  different  polarization 
configurations  involving  the  same  transitions  are  similar;  the  values  are  within  a  factor  of 
two  of  one  another.  In  addition  the  ratio  of  the  Qief(4)  and  Rij(4)  saturation  intensities  is 
0.80  +  0.39,  and  the  ratio  of  the  Qief(8)  and  Rif(8)  saturation  intensities  is  0.63  +  0.21. 
These  values  are  close  to  the  inverse  ratio  of  the  B^^  coefficients  for  these  transitions  of 
1.01  and  0.71,  respectively  Note  that  the  polarization  dependence  of  the  Q-  and  R- 
branch  transitions  differ  sharply.  In  constructing  these  saturation  ratios  we  used  the 


YYXX  saturation  intensities  of  the  MCECA  model,  chose  transitions  that  shared  a 
common  ground  level,  and  assumed  that  effects  arising  from  differences  in  the  excited 
level  degeneracies  and  collisional  rates  are  small. 

C.  Numerical  Comparisons 

The  calculated  AL  and  MCECA  saturation  intensities  for  the  Rif(8)  transition 
using  r„  =  r,2  =  2  X  109  s-1  and  0.405  D  are  0.0031  MW/cm2  and  0.028  MW/cm^, 
respectively.  The  bandwidth-corrected  MCECA  saturation  intensity  is  approximately  an 
order  of  magnitude  larger  than  the  AL  saturation  intensity;  however,  the  MCECA 
saturation  intensity  is  still  an  order  of  magnitude  smaller  than  our  average  experimentally 
determined  value  of  0.32  ±  0.08  MW/cm^.  Order-of-magnitude  agreement  was  also 
observed  in  the  MCECA  experiments.  There  are  several  possible  reasons  for  the 
discrepancy  between  the  calculated  and  measured  saturation  intensities;  although  no  single 
reason  appears  to  be  able  to  account  for  such  a  large  difference. 

The  fact  the  MCECA  model  does  not  include  effects  arising  from  the  spatial  beam 
properties  could  account  for  some  of  the  difference  (possibly  up  to  a  factor  of  ten).  The 
LFR  calculations  discussed  in  Sec.  I,  however,  showed  good  agreement  with  their 
experimental  results  using  a  narrow-bandwidth  laser  even  though  the  calculations 
neglected  the  spatial  characteristics  of  the  excitation  beams.  Therefore  it  is  unlikely  that 
the  Gaussian  nature  of  the  excitation  beams  is  solely  responsible.  The  LFR  calculations 
did  show,  however,  that  the  saturation  intensity  increases  as  the  ratio  2r^^  to  A©  ^ 
decreases.  Therefore  the  fact  that  the  molecules  are  moving  may  account  for  some  of  the 
difference  (less  than  a  factor  of  2).  It  is  also  possible  that  we  have  underestimated  the 
values  of  the  relaxation  rates  by  employing  the  simple  7’’''^  temperature  extrapolation  and 
by  neglecting  the  effects  of  velocity  changing  collisions.  Using  larger  relaxation  rates 
would  increase  the  calculated  saturation  intensity  (possibly  a  factor  of  2). 

Another  indication  of  saturation  is  power-broadening  of  spectral  line  shapes. 

Figure  7  shows  the  average  DFWM  line  width  (fwhm)  of  the  Rie(8),  Rif(8),  R2e(8)»  and 
R2f(8)  transitions  as  a  function  of  if  .  Each  data  point  in  Fig.  7  represents  a  separate 
spectral  scan  from  which  an  average  value  was  obtained.  The  solid  line  is  the  theoretical 
MCECA  model  prediction  of  Eq.  (13)  using  the  laser  bandwidth  of  Au  =  0.18  cm'^  and  no 
adjustable  parameters.  The  experimental  data  is  plotted  vs.  ///“,  and  the  theoretical 
curve  is  plotted  vs.  ///,^  where  /“  and  1°^,  are  the  experimental  and  theoretical 
saturation  intensities,  respectively.  The  simplified  form  of  the  MCECA  model  expressed 
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FIG.  4.7.  Average  DFWM  line  width  (fwhm)  of  the  Rie(8),  Rif(8),  R2e(8),  and  R2f(8) 
transitions  as  a  function  of  .  Each  data  point  represents  a  separate  spectral  scan 

from  which  an  average  value  was  obtained.  The  solid  line  is  the  theoretical  MCECA 
model  prediction  of  Eq.  (13)  using  the  laser  bandwidth  of  Au  =  0.18  cm"l  and  no 
adjustable  parameters.  The  experimental  data  is  plotted  vs.  ///“  and  the  theoretical 
curve  is  plotted  vs.  ijr^,  where  7“  and  7" ,  are  the  experimental  and  theoretical 
saturation  intensities,  respectively. 
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in  Eq.  (13)  is  strictly  valid  for  small  detunings  (A  <  6);  however,  the  agreement  of  the 
model  with  the  data  is  quite  good  even  for  laser  intensities  outside  the  range  of 
applicability  defined  in  Eq.  (18).  In  addition  Eq.  (13)  approximates  a  power-broadened 
Lorentzian-cubed  line  shape  for  laser  intensities  less  than  3  4,,,  and  this  line  shape 
reproduces  the  experimental  data.^^  Therefore  an  estimate  of  the  degree  of  saturation  of 
any  transition  can  be  determined  fi-om  the  power-broadened  line  width  and  Eq.  (13). 

A  comparison  can  also  be  made  between  the  calculated  and  measured  DFWM 
reflectivities.  Steel  and  Lamb^j^l  have  empirically  modified  the  AL  reflectivity  expression 
for  homogeneously  broadened  systems  to  apply  to  Doppler-broadened  systems,  i.e., 

Aco  o  »  2r^g.  The  resulting  expression  also  showed  good  agreement  with  their 
experimental  data.  Steel  and  Lamb  express  the  Doppler-broadened  reflectivity,  Rp,  as 

where  u  is  the  most  probable  speed  of  the  Maxwell-Boltzmann  velocity  distribution  and 
(A  =  0)  is  the  line-center  DFWM  reflectivity  of  Eq.  (9)  for  homogeneously  broadened 
media.  We  note  that  Eq.  (21)  is  consistent  with  the  ratio  of  perturbative  expressions 
(Chapter  3,  Eqs.  (46)-(5 1))  evaluated  in  the  Doppler-broadened  limit  (Aco  p  »  2r^^)  and 
homogeneously  broadened  limit  (A©  P«2r^)p2 

MCECA  have  also  presented  an  approximate  expression  for  the  Doppler- 
broadened  reflectivity  (Eq.  (22)  of  reference  14).  Comparing  that  expression  to  Eq.  (14), 
we  find  that  the  Doppler-broadened  DFWM  reflectivity  using  broad-bandwidth  lasers,  i.e., 
A©  o  >  2i  »  2T^g,  scales  relative  to  the  homogeneously  broadened  reflectivity  as 

=  8  ( A  =  0),  (22) 

where  Rmcbca  the  line-center  DFWM  reflectivity  of  Eq.  (14)  for  homogeneously 

broadened  media. 

The  calculated  Doppler-broadened  reflectivity  for  the  MCECA  model  at  /  =  7"^, 
and  A  =  0  for  the  Rif(8)  transition  is  1 .0  x  10*^  using  Eqs.  (14)  and  (22).  This  value  is 
calculated  for  T=  2114  K,  =  2  x  10^  s't,  AA^  =  1.1  x  10^2  molecules/cm^ 

(population  of  the  F2f(8)  level),  0.405  D,  and  Z,  =  1  mm  (diameter  of  the  flame). 

The  calculated  reflectivity  is  close  to  the  reflectivities  measured  at  the  experimental 
saturation  intensity  (7  =  7“  )  which  ranged  from  0.2  to  1 .2  x  10*^.  See  Fig.  6.  The 
smallest  and  largest  values  in  the  range  correspond  to  data  taken  using  the  XYYX  and 
XXXX  polarization  configurations,  respectively.  Because  of  the  numerous  calibrations 


required  to  evaluate  the  reflectivity,  the  measured  values  are  only  accurate  within  a  factor 
of  three.  Based  on  the  agreement  between  the  calculated  and  measured  reflectivities,  Eq. 
(22)  should  provide  an  order-of-magnitude  estimate  for  the  feasibility  of  detecting 
molecular  species  in  Doppler-broadened  media  with  broad-bandwidth  lasers. 

D.  Relative  Branch  Intensities 

The  spectral  scan  shown  in  Fig.  2  was  taken  using  the  YYXX  polarization 
configuration  at  1.7  times  the  experimental  saturation  intensity  of  the  Qief(8)  transition. 
Using  Eq.  (13)  we  find  that  the  line  widths  of  the  transitions  shown  in  Fig.  2  represent  a 
range  of  approximately  0.9  to  2.3  /“ .  By  varying  by  ±  10%  for  various  //  7°^,  values 
and  numerically  integrating  Eq.  (13)  we  find  that  the  MCECA  model  predicts  a  to 
dependence  for  the  integrated  DFWM  signal  over  the  range  of  0.9  to  2.3  7^, .  This 
absorption  coefficient  dependence  is  consistent  with  the  experimental  results  of  Farrow, 
Rakestraw,  and  Dreier.^^  These  authors  employed  a  multi-mode  laser  similar  to  that  used 
here  and  found  that  the  integrated  DFWM  signal  exhibits  a  reduced  absorption  coefficient 
dependence  for  saturating  laser  intensities. 

We  have  plotted  in  Fig.  8  the  ratio  of  the  integrated  intensities  for  Qi  and  Rj 
transitions  that  share  a  common  initial  state.  In  plotting  the  relative  branch  intensities  as  a 
function  of  Jg  we  assume  that  the  population  relaxation  and  dephasing  rates  do  not  differ 
significantly  even  though  the  Q  branch  involves  Jg  and  Je-Jg  and  the  R  branch  involves  Jg 
and  Jg=Jg-^\.  The  dashed  curve  in  Fig.  8  is  the  saturated  two-level  (STL)  model 
prediction  for  the  Q-to-R-branch  ratio.  We  have  used  a  simple  dependence  because 
the  average  degree  of  saturation  of  the  data  shown  in  Fig.  8  is  1.26  7“  which  corresponds 
to  Bgf  according  to  the  integrated  MCECA  model.  The  dotted  curve  in  Fig.  8  is  the 
Bgfip(B^,B^,B^,B2,Jg,Jef  diagrammatic  perturbation  theory  (DPT)  prediction  where 
Gj  (84 , 8, ,  83 , 82 ;  ,  y, )  is  the  total  geometric  factor  of  Table  3.2.  Recall  that  if  =  F^ 

these  values  apply  to  R-branch  PC-DFWM  signals.  Inspection  of  Fig.  8  indicates  that 
neither  model  is  able  to  account  for  the  fact  that  the  Q-branch  transitions  are  weaker  than 
the  R-branch  transitions  for  high  values  of  Jg  even  though  the  Q-Branch  absorption 
coefficients  are  larger.  See  Fig.  9. 

The  work  of  Farrow,  Rakestraw,  and  Dreier^^  demonstrated  that  there  is  a 
reduced  absorption  coefficient  dependence  (~Bgg)  at  saturating  laser  intensities,  and  our 
polarization  ratio  data  (Fig.  5)  showed  that  the  relative  anisotropy  of  the  total  angular 
momentum  distribution  represented  by  the  total  geometric  factors  persists  in  the  presence 


FIG.  4.8.  Ratio  of  the  integrated  intensities  for  Qj  and  Rj  transitions  that  share  a 
common  initial  state  as  a  function  of  7g.  The  dashed  curve  is  the  saturated  two-level 

model  (STL)  prediction  for  the  Q-to-R-branch  ratio,  the  dotted  curve  is  the 
B^,Gj(e4,e,,E3,e2;7^,7e)^  diagrammatic  perturbation  theory  (DPT)  prediction,  and  the 
solid  curve  is  the  empirical  (WZR)  prediction. 


FIG.  4.9.  (a)  The  Einstein  absorption  coefficients  as  a  function  of  7g  for  the  Rj-,  Qp.  and 
Pj-branch  transitions  of  CH  labeled  -0-,  and  -A-,  respectively.  The  absorption 
coefficients  for  the  R2-,  Q2--  and  P2-branch  transitions  are  nearly  the  same  because  both 
the  X  and  A  states  rapidly  approach  Hund's  case  (b)  coupling  as  rotation  increases,  (b) 
Total  geometric  factors  of  Table  3.2  for  the  YYXX  polarization  configuration  as  a 
function  of  7g  with  the  same  conventions  as  (a). 
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of  saturating  fields.  Therefore,  we  have  also  plotted  in  Fig.  8 

which  is  the  solid  curve  (WZR).  In  Fig.  8  it  is  apparent  that  both  a  reduced  absorption 
coefficient  dependence  and  a  total  geometric  correction  are  necessary  to  account  for  the 
observed  Q-to-R-branch  ratios.  Note  at  Jg  =  4.5  the  DPT  and  WZR  curves  are  equal,  and 
both  curves  are  in  excellent  agreement  with  the  observed  ratio.  The  agreement  at  Jg  =  4.5 
results  fi’om  the  fact  the  absorption  coefficients  for  the  Qi(4)  and  Ri(4)  transitions  are 
equal,  so  the  ratio  is  determined  exclusively  by  the  geometric  correction.  The  agreement 
at  other  values  is  not  as  good  because  the  Qj-  and  Ri-branch  absorption  coefficients 
are  not  equal,  and  therefore  the  transitions  are  not  saturated  to  the  same  degree. 

Once  more  the  simple  introduction  of  the  geometric  factors  is  able  to  explain  the 
observed  signal  ratios  with  an  accuracy  in  the  range  of  10-30%.  We  noted  in  Chapter  3 
that  at  saturating  laser  intensities  higher  order  moments  of  the  angular  momentum 
distribution  contribute  to  the  DFWM  signal;  however,  these  data  suggest  that  the  first 
three  moments  are  sufficient  nevertheless  to  predict  relative  signal  strengths  at  laser 
intensities  up  to  approximately  twice  the  saturation  intensity. 

E.  Relative  Population  Distributions 

In  general  DFWM  signals  must  be  corrected  for  differences  in  collision  rates, 
absorption  coefficients,  and  geometric  factors  before  relative  population  distributions  can 
be  determined.  Previous  saturated  DFWM  temperature  determinations  employing  broad- 
bandwidth  excitation,  however,  have  shown  good  agreement  with  thermocouple^  and 
coherent  anti-Stokes  Raman  scattering  ^  measurements  without  applying  collisional 
and  geometric  corrections;  the  only^-dependent  corrections  made  in  these  determinations 
were  for  relative  absorption  coefficients  {Idfwm  ^  B^).  In  these  experiments  the  laser 
intensity  was  in  excess  of  the  saturation  intensity,  and  the  analysis  only  involved  DFWM 
signals  from  members  of  a  single  branch.  We  offer  the  following  explanation  for  the 
apparent  insensitivity  of  the  DFWM  signal  to  collisions  and  geometric  corrections  in  these 
experiments  and  present  the  results  of  a  rotational  temperature  analysis  of  the  CH  radical 
performed  under  similar  conditions. 

The  reduced  sensitivity  to  collisions  is  discussed  first.  The  LFR  calculations  (Sec. 
I)  for  narrow-band  DFWM  of  Doppler-broadened  systems  predict  that  the  line-center 
DFWM  signal  is  nearly  independent  of  collisions  at  approximately  2  4" , .  This  assertion  is 
substantiated  by  their  experimental  data  (Fig.  5  of  reference  3  and  Fig.  15  of  reference  16) 
taken  with  a  narrow-band  excitation  source.  We  can  understand  this  result  by  considering 
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the  how  collisions  affect  the  DFWM  signal  for  Doppler-broadened  systems.  Using  the 
total  collisional  dependence  of  the  Doppler-broadened  reflectivity  of  Eq.  (21),  we  find  that 
the  line-center  DFWM  signal  of  the  AL  model  has  a  zero-order  collisional  rate 
dependence  at  approximately  0.5/^,.  This  value  is  in  qualitative  agreement  with  the  LFR 
calculated  value  of  2 7^°,.  Without  employing  the  additional  collisional  dependence  ofEq. 
(21)  for  Doppler-broadened  systems,  the  line-center  DFWM  signal  of  the  AL  model  is 
independent  of  collisions  only  when  I »  7®^,.  This  result  is  does  not  qualitatively  agree 
with  the  LFR  calculations. 

By  comparison  we  find  using  the  total  collisional  dependence  of  the  Doppler- 
broadened  reflectivity  of  Eq.  (22)  that  the  line-center  DFWM  signal  of  the  MCECA  model 
has  a  zero-order  collisional  rate  dependence  at  approximately  0.67^° , .  The  applicability  of 
Eq.  (22)  to  our  experiments  is  validated  by  the  fact  that  the  Doppler-broadened  correction 
is  necessary  to  account  for  the  magnitude  of  our  experimentally  measured  reflectivities. 
Hence  it  is  reasonable  to  assume  that  broad-bandwidth  DFWM  signals  of  Doppler- 
broadened  systems  exhibit  a  reduced  collisional  dependence  at  laser  intensities 
approximately  equal  to  27°^,. 

The  independence  of  DFWM  signals  to  geometric  corrections  results  from  that  the 
fact  that  the  total  geometric  factors  are  essentially  constant  for  members  of  a  single 
branch.  As  is  shown  in  Fig.  9,  the  total  geometric  factors  differ  markedly  for  different 
branches,  are  J  dependent  for  lowJ^  values,  but  are  constant  for  each  branch  for  values  of 
Jg  greater  than  three  or  four.  We  expect  that  once  the  geometric  factors  approach  a  high- 
/  limit  then  the  only  remaining  7-dependent  corrections  to  saturated  DFWM  signal 

intensities  for  members  of  a  single  branch  would  arise  from  differences  in  absorption 
coefficients.  We  suggest  that  this  behavior  is  the  primary  reason  for  the  remarkable 
success  of  applying  simple  absorption  coefficient  corrections  to  determine  relative 
population  distributions  from  saturated  DFWM  signals. 

We  further  emphasize  this  important  point  by  presenting  the  results  of  a  saturated 
DFWM  rotational  temperature  analysis.  In  our  rotational  temperature  determination  we 
use  Q-branch  transitions  recorded  using  the  YYXX  PC-DFWM  experimental 
configuration  with  the  laser  intensity  maintained  at  approximately  1.7  times  the 
experimental  saturation  intensity  of  the  Qief(8)  transition.  The  average  degree  of 
saturation  of  all  the  Q-branch  transitions  used  in  the  analysis  is  1.67“  with  a  range  of 
approximately  0.9  to  2.3  7“  as  determined  from  the  power-broadened  line  widths.  This 
average  degree  of  saturation  corresponds  to  dependence  for  the  integrated  DFWM 
signal  intensities  according  to  the  MCECA  model.  A  Boltzmann  plot  of  integrated  Q- 
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branch  transitions  for  2.5  <  <  20.5  is  shown  in  Fig.  10.  The  lower  line  (STL) 

represents  a  least-squares  linear  regression  analysis  using  data  corrected  for  the 
dependence,  and  the  upper  line  (WZR)  represents  data  corrected  for  the 

dependence.  In  both  linear  regression  analyses  data 
corresponding  to  3.5  <  <  20.5  of  resolved  Qj  and  Q2  transitions  are  used  and  no 

collisional  corrections  are  applied. 

We  find  that  the  geometric  corrections  only  change  the  fitted  temperature  by  0.5%. 
This  result  is  expected  since  the  total  geometric  factors  of  the  Q-branch  are  nearly 
constant  over  the  range  3.5  <  <  20.5.  The  data  at  the  smallest  term  energies  (Jg=2.S- 

4.5),  however,  deviate  significantly  from  the  data  at  larger  term  energies.  Inspection  of 
Fig.  9  shows  that  the  Q-branch  the  absorption  coefficients  are  7  dependent  for  low  values 
ofJg.  We  showed  for  the  Q-to-R-branch  ratios  in  Sec.  IVD  that  the  applied  corrections, 
i.e.,  are  only  accurate  to  10-30%  when  the  relative 

absorption  coefficients  and  geometric  factors  differ.  Therefore  using  DFWM  signal  fi'om 
members  of  a  single  branch  and  lvalues  greater  than  three  or  four  greatly  reduces  the 
error  in  determining  the  rotational  temperature  because  both  the  absorption  coefficients 
and  total  geometric  factors  are  essentially  J  independent. 

The  degree  of  saturation  of  the  data  used  in  the  analysis  ranges  from  approximately 
0.9/^^,  at  ,^=3.5  to  2.3  7“  ato^=20.5.  This  range  of  saturation  intensities  corresponds  to 
approximately  a  60%  decrease  in  the  population  relaxation  rate,  r„,  fi'om  *7g=3.5  to 
Jg=20.5  using  cc  bT^jB^^ .  A  decrease  in  the  relaxation  rate  with  increasing  Jg  is 
expected  for  CH  in  this  flame  because  and  R  generally  decreases  as  the 

spacing  between  adjacent  rotational  levels  increases.  Hence  if  the  DFWM  signal  were 

sensitive  to  collisions,  a  60%  variation  in  the  collisional  rate  would  significantly  affect  the 
temperature  obtained  from  the  Boltzmann  analysis.  For  example  the  temperature  would 
increase  by  approximately  160  K  (>5%)  if  the  DFWM  signal  were  inversely  proportional 
to  the  collisional  rate. 

The  temperature  obtained  fi'om  the  linear  regression  analysis  of  the  data  corrected 
for  the  0^(84, e,,8j,82;./g,i/e)^  dependence  is  2774  ±  58  K.  This  value  is  in 
consistent  (within  5%)  with  our  previously  reported  ^  2  cjj  vibrational  temperatures  of 
2642  ±  99  K  and  2882  ±  123  K  obtained  by  LIF  and  DFWM,  respectively.  Because  of 
the  specific  characteristics  of  the  X  and  A  states  of  CH,  vibrational  temperatures  do  not 
need  to  be  corrected  for  relaxation,  polarization,  or  saturation  effects.  This  interesting 
aspect  of  high-pressure  CH  spectroscopy  is  discussed  in  detail  in  reference  12.  Therefore 
the  5%  agreement  between  the  rotational  temperature  and  vibrational  temperatures 
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FIG.  4.10.  Boltzmann  plot  of  integrated  Q-branch  transitions  for  2.5  ^  <  20.5  taken  in 

the  YYXX  polarization  configuration.  The  power-broadened  line  widths  of  the 
transitions  used  in  the  analysis  indicate  that  the  average  degree  of  saturation  of  these  data 
is  1.6  /“,.  The  lower  line  (STL)  represents  a  least-squares  linear  regression  analysis  using 
data  corrected  for  the  dependence,  and  the  upper  line  (WZR)  represents  data 
corrected  for  the  Gj  (£4 ,  e, ,  83 ,  Sj ;  dependence.  In  both  linear  regression 

analyses  data  corresponding  to  3.5  20.5  of  resolved  Q2  and  Q2  transitions  are  used 

and  no  coUisional  corrections  are  applied. 
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Term  Energy  (cm 


suggests  that  the  DFWM  signal  does  in  fact  exhibit  a  reduced  (less  than  linear)  collisional 
dependence  when  I  ~  /“ . 

Finally  we  comment  on  the  rotational  temperature  error  of  58  K  (2%).  This  error 
is  more  indicative  of  the  precision  of  these  measurements  rather  than  the  accuracy  because 
several  significant  assumptions  were  made  in  interpreting  of  the  observed  DFWM  signals 
using  the  MCECA  model.  Based  on  these  assumptions  and  the  results  presented  in  this 
chapter,  a  more  realistic  estimate  of  the  accuracy  of  these  measurements  is  140  K  (5%). 
This  conclusion  is  also  supported  by  the  results  of  another  ongoing  project  at  Stanford 
which  involves  the  investigation  of  trace  radical  species  such  as  CH  and  C2  in  an 
atmospheric-pressure  diamond  synthesis  reactor.^^"^^  The  reactor  consists  of  an  rf 
inductively-coupled  argon  plasma  seeded  with  hydrogen  and  methane.  These  studies  have 
focused  on  measuring  CH  and  C2  concentration  profiles  and  CH  vibrational  and  rotational 
temperatures  using  saturated  DFWM.  The  results  of  these  measurements  are  consistent 
with  the  results  of  numerical  simulations  of  the  deposition  environment  over  temperature 
range  of 2500-4000  K.  Furthermore  the  agreement  between  the  CH  vibrational, 
rotational,  and  simulated  temperatures  has  consistently  been  ~5%.^^"^® 

V.  CONCLUSIONS 

We  have  experimentally  evaluated  the  effects  of  polarization  in  DFWM 
spectroscopy  from  the  weak-  to  the  strong-field  limit.  The  effects  of  polarization  were 
rationalized  using  the  ideas  developed  in  Chapter  3  for  the  weak-field  limit.  The  key 
results  are  that  the  saturation  intensity  is  relatively  independent  of  polarization  and  that  the 
DFWM  spatial  effects  represented  by  the  weak-field  geometric  factors  of  Chapter  3  are 
applicable  at  laser  intensities  up  to  twice  the  saturation  intensity.  Furthermore  we  showed 
that  polarization  ratios  measured  with  weak  fields  can  be  used  to  infer  important 
relaxation  information. 

The  error  in  applying  the  weak-field  geometric  factors  to  saturated  DFWM  signals 
was  found  to  be  in  the  range  of  10-30%.  This  degree  of  accuracy  may  be  sufficient  for 
determining  energy  transfer  rates.  In  particular  the  sub-Doppler  resolution  offered  by  PC- 
DFWM  would  be  advantageous  for  resolving  energy  transfer  rate  differences  between 
different  spin-orbit  and  A-doublet  components  that  are  difficult  to  resolve  using  LIF  at 
these  temperatures.  In  determining  relative  population  distributions  for  temperature 
measurements,  however,  10-30%  accuracy  is  often  not  sufficient. 
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We  showed  in  Sec.  IV  that  if  the  absorption  coefficients  and  the  total  geometric 
factors  of  a  set  of  transitions  are  nearly  the  same,  relative  population  distributions  accurate 
to  5%  can  be  extracted  from  saturated  DFWM  signals.  Such  behavior  is  expected  for 
members  of  a  single  branch  once  the  absorption  coefficients  and  the  total  geometric 
factors  reach  a  highJ-limit.  Therefore  small  values  of  7  should  be  avoided  in  the  analysis 
when  accuracy  is  of  primary  importance.  This  result  is  an  important  finding  because  it 
substantiates  a  significant  assumption  in  using  two-level  nondegenerate  models  to 
determine  relative  internal-state  distributions  of  molecules.  We  also  note  that  this  finding 
is  the  same  as  that  of  Altkom  and  Zare^^  for  saturated  LIF. 

We  supported  this  finding  by  presenting  a  rotational  temperature  analysis  of  the 
CH  radical  in  an  atmospheric-pressure  flame.  We  employed  a  novel  polarization 
configuration  in  this  experiment  which  offered  true  zero-background  detection  and  greatly 
discriminated  against  other  intensity  grating  contributions  to  the  DFWM  signal.  The  laser 
intensity  was  maintained  at  approximately  1.7  times  the  experimental  saturation  intensity 
of  the  Qief(8)  transition,  and  a  PC-DFWM  spectrum  of  the  Q  branch  of  the  CH 
a2a(v -0)-X2n(v"=0)  system  was  obtained  in  the  YYXX  polarization  configuration. 
From  a  Boltzmann  analysis  a  rotational  temperature  of 2774  ±  140  K  was  obtained.  This 
rotational  temperature  is  in  good  agreement  with  the  average  LIF  and  DFWM  vibrational 
temperature  of  2762  ±  160  K  previously  reported  by  our  laboratory. 
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APPENDIX 


Calculation  of  the  Polarization  Tensor  Products 


Math,  Math,  and  more  Math 


In  this  appendix  we  establish  the  notation  used  for  expressing  the  polarization  unit 
vectors  and  evaluate  the  polarization  tensor  products  necessary  to  calculate  the  DFWM 
signal  intensity  for  experimental  configurations  not  treated  in  the  text.  In  what  follows,  we 
rely  on  the  definitions  of  symbols  already  introduced. 

An  arbitrary  vector  r  can  be  expressed  in  Cartesian  coordinates  as 

''A (Al) 

and  its  complex  conjugate  r*  as 


r=Ke^+r*e^+r*e^,  (A2) 

where  e^,  e^,  and  are  unit  vectors  and  r,,  and  are  the  standard  Cartesian 
components  of  the  vector  r .  The  conventions  for  relating  the  unit  vectors  and  the  standard 
components  of  the  Cartesian  (real)  basis  to  the  spherical  tensor  (complex)  basis  are  defined 
as  follows: 

=  .  (A3) 

and 

±  K ) ,  =  z; .  (A4) 


It  follows  that 


(A5) 


and 


(r’)J>  =  (-l)*[r.»]‘  (A6) 

where  q  =  0,  +1.  Equations  (A3)  and  (A4)  in  turn  define  the  conventions  for  relating  the 
unit  vectors  and  the  standard  components  of  the  spherical  tensor  basis  to  the  Cartesian 
basis,  i.e., 

+  ^-V)  ’  =  ^0.  (A7) 

and 


-rO)  r  - _ 

'-l  /  >  V  V+1 


(A8) 


Consequently,  we  have  for  the  Cartesian  basis 


♦  *  « 
e=e  ,  e  •  e  ,  =  0 

p  p  ^  p  p  pp  ^ 


(/•*). =[^1*, 


(AlO) 


where  p-x,y,z.  Therefore  adopting  the  standard  definitions,  we  are  left  with  the 
conclusion  of  expressing  an  arbitrary  vector  r  and  its  complex  conjugate  r*  in  the  spherical 
tensor  basis  as 


(All) 


r  =  eo  +  .  (A12) 

Hence  we  can  express  an  orthonormal  set  of  unit  vectors  in  a  generalized 


notation  as 


(A13) 


where  the  standard  components  of  the  unit  vectors  are  defined  as 

(S/)|  =  '/■«(  .  ((6,)ir=e;-c"  ,  e,. 


(A14) 


Special  care  must  be  taken  in  interpreting  Eq.  (A14).  We  take  (e^ ),  to  represent  the  i 
component  of  the  unit  vector  e^,  [(By)J*to  represent  the  complex  conjugate  of  the  / 
component  of  the  unit  vector  Sj,  and  (e*  ),  to  represent  the  /  component  of  the  complex 
conjugate  of  the  unit  vector  i.e.,  e*. 

For  light  in  a  pure  state  of  arbitrary  polarization  propagating  along  the  space-fixed 
Z-axis  (such  a  choice  is  always  possible),  its  polarization  unit  vector  and  its  complex 
conjugate  may  be  represented  by 

Sj  =cos<j>ye^  +  e*^sin<j)^.e_y  =-^  cosily -/e'^^  sinclty  e*_^  --^j^cos(|>y-i-/e*^^sin<|)^  c*,  (A15) 


e’ =cos(|)j  e, +e  *‘’''sin<|)^.c^  =-^|^cos(j)^.-»e  *“''  sin<|)^.Je*, -^|^cos(|)^. +  7c  sm(j)^Je']  (A16) 

where  (t)^  ranges  from  0  to  7t  and  5^  is  the  phase  necessary  to  describe  elliptical 
polarization.  Reference  to  Eqs.  (A13)-(A16)  shows  that 

(8^.),  =COS(J)^  ,  (8*),=COS(l)^.  , 
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(AI7) 


For  completeness  we  define  left  circularly  polarized  light^  ((j)^  =  7t/4,  5^  =  7t/2)  as 

^L=±(e,+iey)  =  e_,,  (A18) 


right  circularly  polarized  light  =  7c/4,  6^  =  -n/2)  as 

^R  =  ±(e,-iey)  =  -el^ , 

and  light  of  arbitraiy  linear  polarization  (5^  =  0)  as 


(A19) 


8^.  =cos(|)^e,+sin(|)^.e^=^ 


e  -e 


(A20) 


The  polarization  tensors  of  Chapter  3,  Eq.  (51)  in  the  text  are  expanded  and 
expressed  in  the  spherical  tensor  basis  as 


F{z„Z„Zj,Z,-K)=  ^(-l)e(8:)i'>  (E,)«g  {z\)%  {z,)% 
6.9.9' 


1  1  K 

'  1 

1  k" 

-q  q+Q  -Q^ 

.-9' 

q'-Q  Q, 

,  (A21) 


where  the  terms  in  parenthesis  are  3-y  symbols.  The  polarization  tensors  in  the  Cartesian 
basis  are  readily  worked  out  using  Eq.  (A4)  and  (A21).  The  resulting  expression  does  not 
have  a  compact  form,  so  the  polarization  tensors  are  presented  for  each  value  of  K,  i.e., 

-^(64,Bj,8y,82;0)  =  — I  [(84)^  (®4)>’  (^7)^'*'  (®4X  (®;)r  ] 

X  [(«; ).  (8, ),  +  (e;  ),  (6, ),  +  (e;  ),  (6, ),]  },  (A22a) 

f'(S4,e,.8,.8j;l)  =  -i|  [(e;),  (s,),-  (si),  (e^),]  x  [(e;),  (e,),,-  (z\\ (8,).] 

+  [(six  (8;)x  -  (e;).  (8^.),]  X  [(e;),  (8,  \  -  (8;),  (8,  ),] 

-  [(el),  (8,),  -  (8;),  (8,),]  X  [(8^),  (8,  ),  -  (8^),  (8,  ),]  |,  (A22b) 


F(84,8,,E^.,E2;2)  =  ||  [2(8;X  (6;).-  «),  (8,),-  (s;)^  (8^.)^] 

x[2(eD.(e,).-(8;),(8,),-(8;),(8,)^]  } 

+  [«)z  (e;)x+  (el)x  (e,x]  x  [(e^),  (e,  )x+  (ej),  (e,).] 

+  [(Oz  (S;  )y  +  (£4);,  (s>)z]  X  [(8^),  (8,  )y  +  (s^)^  (e,  ),] 

+  [(e4)x  (6;  )x  -  (sl)^  (8;)^]  X  [(e^),  (8(  X  -  (82)^  (8,  x] 

+  [(8;)x  (8,)^+  (8:)^  (8^.),]  X  [(e^),  (e,)^+  (e;)^  (e,),]  |.(A22c) 

Equations  (A21)  and  (A22)  can  be  used  to  obtain  the  DFWM  signal  intensity  along  84,  i.e., 

the  DFWM  signal  direction.  It  is  often  possible  to  predict  the  direction  of  84  by  invoking 
symmetry  arguments  ofX^\  These  equations  also  apply  for  polarized  detection  in  which 
case  84  is  the  polarization  axis  of  the  DFWM  detector. 

For  experiments  in  which  it  is  difiicult  to  predict  the  DFWM  signal  direction  by 
symmetry  arguments  or  for  unpolarized  detection,  it  is  more  useful  to  express  the  DFWM 
signal  intensity  in  terms  of  its  components.  Rewriting  the  DFWM  signal  intensity  of 
Chapter  3,  Eq.  (19)  in  Sec  niA  in  terms  of  the  components  of  the  generalized  basis  of 
orthonormal  unit  vectors  defined  by  Eqs.  (A13)  and  (A14)  yields 

=  (A23) 


where  pf  is  the  scalar  amplitude  of  the  j  component  of  P^^^{r,i)  and  is  defined  as 

^  ~  ^Fmi  *8;  •  (A24) 

Equation  (A24)  can  be  rewritten  using  the  notation  of  Chapter  3,  Eqs.  (45)-(5 1)  as 


e,e;e. 


K 
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K 


K  ' 


where  all  of  the  terms  have  their  previous  meanings. 

Equation  (A25)  can  also  be  expressed  in  terms  of  its  standard  components.  In  the 
spherical  tensor  basis  =  e*,  and  the  polarization  tensors  in  Eq.  (A25)  take  on  the  form 

K 

Q=-K 


O 


/Ov 

2  ®e, 


n 


(A26) 


Substitution  of  =  c,  •e_^,  =  into  Eq.  (A21)  gives 


F(e;,e„ej.z,-.K)=  2](-l)«  (e, )?-£ 

e.«' 


[-q  q  +  Q  -Q)[-q'  q'-Q  Q] 


(A27) 


which  can  be  used  in  Eq.  (A25)  to  generate  the  expression  for  the  scalar  amplitude  of  the  q 
component  of  the  nonlinear  polarization  P^^^{r,i).  The  resulting  expression  is  equivalent  to 
Eqs.  (33)  and  (34)  of  reference  2  derived  using  a  density  matrix  approach  in  the  weak-field 
limit.  In  the  Cartesian  basis  8^  =  e*,  and  substitution  of  (s^)^-  =  =  5^^,,  into  Eq. 

(A22)  for  p=y  gives  the  Y  component  of  the  polarization  tensors  as  a  function  of  K, 
namely, 


F(c;,e,,8  .,e2;0)  =  i|  [(s^)^  (s^),  (e,),  +  (e^)^  (e*)^  (e,),  +  (e^^  (e^),  (e,),]  |,  (A28a) 
F(e;,E,,8^,E2;l)  =  ||  [(8^.),  (8^),  (8,)^  -  (8^.),  (8^)^  (8,),] 

+  [(£;).  (6. ).  -  (e,X  (6.  X]  (A28b) 


and 


F(e*,E,,E^.,E2;2)  = 


^^2)z  (.^1)2  (^j)y  (^2)1  i^i)x  (®>)>>  (^2)^;  (^i  )>>]  | 

+||  [(eyX(eDr(e,)y+(e>X(eD>(e.)z] 
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(A28c) 


The  X  and  Z  components  are  readily  evaluated  by  reiterating  the  procedure  for  p=x  and  z, 
respectively. 

For  the  near  collinear  phase  matching  geometries  and  pure  polarization  states 
described  in  this  thesis,  only  one  standard  component  is  nonzero.  This  result  can  be 
verified  by  substituting  the  polarizations  of  the  input  fields  E^,  E2,  and  into  Eq.  (A26) 
and  (A28)  using  the  definitions  of  Eqs.  (A18)  and  (A19)  for  circularly  polarized  light  and 
Eq.  (A20)  for  (J)^  =  0  and  =  7c/2  linearly  polarized  light.  These  substitutions  yield  the 
results  presented  in  Table  3.5.  For  noncollinear  phase  matching  geometries,  mixed 
polarization  states,  or  both,  all  of  the  components  must  be  determined. 

The  above  expressions  refer  to  the  standard  components  of  the  electric  fields  when 
all  the  fields  are  defined  with  respect  to  a  common  reference  frame.  For  example,  in  the 
case  of  noncollinear  phase  matching  geometries,  all  the  fields  must  be  rotated  into  a 
common  frame  of  reference  (we  suggest  the  detection  frame)  before  the  standard 
components  can  be  determined.  Such  frame  rotations  are  readily  carried  out  in  the 
spherical  tensor  basis  using  Wigner  rotation  matrices  or  in  the  Cartesian  basis  using 
direction  cosine  matrices.  For  more  information  on  this  topic,  see  reference  3,  Chapter  3. 
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